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The Problem of the Turbulent Boundary Layer 


By Donatp Cotes, Pasadena, California?) 


NNT EE ee ee RT et eee 


1. Introduction 


Historically, experimental measurements of the magnitude of surface fric- 
| _tion have usually preceded the development of adequate methods of prediction, 
_ and the writer’s study of supersonic boundary layers?) is a case in point. In the 
present survey the turbulent boundary layer problem is examined in some 
_ detail from the experimental point of view, in an effort to identify the features 
which must eventually appear in any successful analysis of boundary layer flow. 
s If experimental measurements could be compared with an exact solution of 
the equations of fluid motion for the particular model geometry involved, any 
_ discrepancies could be assigned to error in the meastrrements or to failure of the 
4 equations to describe the behavior of the fluid in question. In practice, however, 
a comparison of physical and mathematical models for viscous fluid flow over 
- a flat plate can be achieved only in an asymptotic sense, first because of mathe- 
- matical approximations which are not uniformly valid in the entire flow field, 
_ and second because of physical imperfections which are idealized or ignored in 
the usual theoretical treatment. 
E The idea of an asymptotic correspondence between actual and assumed con- 
- ditions is essential to any definition of the boundary layer problem. For laminar 
_ flow this idea, which is explicit in the boundary layer concept of PRANDTL, is 
acquiring an increasingly precise mathematical formulation, and the recent 
; development of methods for improving the theoretical analysis in this respect 
| is a conspicuous sign of current interest in the question [1]). 

On the other hand, the present inadequate state of knowledge of turbulent 
~ shear flows requires the corresponding definition of the turbulent boundary 
layer problem to be almost completely heuristic, so that the validity of a 

particular analysis can only be tested in the light of experimental evidence. It 
is proposed to illustrate this difficulty by a brief discussion of the flat plate 
| Boundary layer, using an analysis which is frankly phenomenological and which 
is therefore less vulnerable to criticism of concept than an analysis which 
assumes one or another mechanism for the various turbulent transport pro- 


cesses. 


1) Guggenheim Aeronautical Laboratory, California Institute of Technology. 
a 2) See Measurements of Turbulent Friction on a Smooth Flat Plate in Supersonic Flow, to be 


- published separately. 
3) Numbers in brackets refer to References, page 201. 
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While the review given here is designed to provide a point of departure for 
the study of compressible boundary layers, it should be remembered that the 
purpose of the discussion is always to facilitate the interpretation of experimen- 
tal data, rather than to suggest or to support a particular theoretical analysis 
of the turbulent boundary layer. 


2. General Considerations 


The remarks of the present paper refer to steady two-dimensional mean flow 
past a smooth flat plate at sufficiently large Reynolds numbers so that the 
boundary layer hypothesis may be assumed to apply. It is further assumed 
that there is no heat transfer between the plate and the fluid, and that the 
pressure, density, and velocity of the external flow are uniform. 

Early experiments, particularly those carried out in towing tanks, were in 
most cases designed to determine a mean friction coefficient from its definition 
in terms of the dynamic pressure and of the drag per unit width of one side 
of a flat surface, 


D(x) = %q C(x) , (1) 


while an alternative method makes use of the equivalence of the drag force and 
the rate of momentum change in the boundary layer or in the wake, 


D(x) = [ 9 u(y —w) ay. (2) 


Finally, the total drag may be equated to the sum of the downstream compo- 
nents of the integrated tangential and normal forces, 


x 


D(x) = ie T(x) dx + Dol(x) , (3) 


0 


where the second term on the right represents the drag of any isolated roughness 
elements which may be present on the plate surface. 
The first two equations above are usually combined in the form 


Colz) = = O(a) , (4) 


in which the boundary layer momentum thickness @ is defined by the definite 
integral 


O(x) = | eg ae (1 - Z) dy. (5) 
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The first and third equations i in turn peowide a gelowiouship between drag and 
local friction, 


_ 1 f tyl4) Dal) 
G(x) == [9 ax + A, (6) 
so that 
Tole) |» €O(2) Ie CADE) 
q anger Ge ano,” (7) 


The notation in these expressions is chosen deliberately in order to empha- 
_ size that x, t,,/¢, Cp, and @ are experimentally measured quantities. The essence 
of most boundary layer research consists of the assumption that there exist. 
corresponding uniquely defined quantities €, c,, cr, and 0, which may be asso- 
ciated with an ideal boundary layer, and of the attempt to determine the latter 
quantities by necessarily imperfect experimental techniques. That is, by defi- 
nition, 
d dé 


Op= 0/{6) = eS a(t) = 2: (8) 
é 
or = orl) = ¢ | o(€) dB = 2 (9) 
0 
1 1 : 
0 = (6) = 5 F op(é) => | on(é) ab -. (10) 


In order to establish a connection between the experimental and ideal 
quantities defined above, it is necessary to introduce an additional hypothesis. 
The assumption usually made is that the effect of conditions near the leading 
edge on the surface friction far downstream will eventually vanish, in the sense 
that the function O(x) for a smooth plate is the ideal relationship 6(&) subject 
to a possible displacement of the streamwise coordinate; in other words, down- 

stream of a particular point x= s, 


Oe = %e)-=-O(x)’. (11) 


Alternatively, a similar but less restrictive assumption may be made in 
terms of the local surface friction; 


og — np) = 0. (12) 


In view of equations (4), (6), and (10), it is clear that the hypothesis (11) 
implies (12), while the converse is not true. Given the hypothesis (12), then 


0(% — %) = O(x) — Oy, (13) 
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where 
S—%» 


Oo gt rea 1 
0) =F" + al WO) dn 5 f en6) a8. (14) 
0 


The definition (14) assumes that the form drag Dy is independent of 
downstream of x = s, as will be the case for example if Dy is associated with 
finite leading edge curvature or with the presence of boundary layer tripping 
devices on the surface of the plate. The experimental evidence suggests that the 
assumption (11) is probably sufficiently general in practice; i.e. that OQ, van- 
ishes. The expression (14) thus serves to connect the drag of a tripping device 
with the position at x = %, of the apparent origin of the turbulent boundary 


layer. 
Finally, equations (8), (9), and (11) imply 
(C) 
Ca(% — Xo) = 2 % : (15) 
dO(x Re 
o,(% — x9) = 2 oo) = Tele (16) 


The remarks above apply to both laminar and turbulent boundary layers 
for arbitrary surface heat transfer and free stream Mach number. The definition 
of momentum thickness in equation (5) above is not rigorously correct for 
turbulent flow. However, since the usual experimental methods for measure- 
ment of mean velocity involve errors of comparable magnitude, no distinction 
will be made here. 

The analysis leading to the equations (13) to (16) is carried out formally 
and in detail both because-of its generality and because the considerations 
mentioned are fundamental to any evaluation of experimental data. The pre- 
sence of the parameter @, in equation (13) immediately raises the question of 
uniqueness, especially for the turbulent boundary layer. The term uniqueness, 
in the sense used here, refers to the validity of a comparison of the ideal function 
O(cs), obtained by eliminating the variable € between the functions 6(€) and 
c,(€) of equations (8) and (10), with its experimental analog O(t,,/¢). The sug- 
gestion that the slope of the experimental curve @(x) may not imply a unique 
value of the ordinate @ is precisely equivalent to the introduction of the 
constant of integration @5. Since 0 increases with &, it is apparent that suffi- 
ciently far downstream the parameters x) and @) must always become small 
compared to « and @ respectively. However, the sense in which an asymptotic 
correspondence between the real and ideal boundary layers is defined by this 
property is not the sense in which the concept is most appropriate to the 
experimental problem, as will be emphasized later. 

The matter of uniqueness can be examined more closely for the laminar 
boundary layer, since the theoretical solution is known within the limitations 


} 


yy 


> 


_ 
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of the boundary layer approximation, and has been verified experimentally 
for low-speed flow by DRYDEN [2], LrepMANN [3], DHAWAN [4], and others. The 
laminar equations are satisfied by a function of a single independent variable 
equivalent to yIVE. The consequent similarity of the various velocity profiles 
along the plate implies that the boundary layer momentum thickness is in- 
versely proportional to the slope of the velocity profile at the wall, so that the 
functional relationship 6(c;) mentioned above is 


cy R, = const = (0-664)?, (17) 


for the laminar boundary layer with constant density, where Ry = u, 6/». 


For the turbulent boundary layer, unfortunately, the question of uniqueness 


_ lies under a darker cloud than can be dissipated by theoretical arguments 


alone. Further consideration of this problem will therefore be deferred until it 


_ becomes convenient below to cite experimental evidence in support of an anal- 
_ ysis, based on considerations of similarity, of the turbulent boundary layer. 


3. Functional-Similarity 


In the early attempts of PRANDTL, VON KARMAN, and TAYLOR to formu- 


late a theoretical treatment of turbulent shear flows, the molecular transport 
_ process of the laminar regime was used as a model for a hypothetical turbulent 
_ transport process. These analyses have in common the concept of a characteris- 


tic length, somewhat analogous to the mean free path in a gas, whose magnitude 
is estimated in terms of some physical property of the mean flow. 

Not only have these theories made possible a.careful extrapolation of empir- 
ical knowledge of turbulent shear flows to conditions outside the contempo- 


rary range of experiment, but they have provided a framework within which 
the effects for example of wall roughness [5] could be fitted. However, it is now 
generally recognized that recent detailed measurements of the structure of 


turbulent flow are making tenancy of the mixing hypothesis progressively 
more uncomfortable [6], [7], although the primitive state of experimental 
research has so far prevented replacing the theory with one more satisfactory. 

The mixing analogy for turbulent flow will not be discussed further here, 
since many detailed summaries are available. Neither will any review be at- 
tempted of the current experimental investigations of LAUFER [8], KLEBANOFF 
[9], TOWNSEND [10], [11], and others, which it is to be hoped will lead first to 
a knowledge of the structure, and finally to an understanding of the mechanism, 
of turbulent shear flow. 

In the absence of an adequate theory describing the turbulent mechanism, 
the present discussion will attempt to organize, rather than to explain, the 


experimental evidence which is available for the low-speed turbulent flat plate 


boundary layer. The general analysis is based on the principle of physical 
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similarity, using certain demonstrable properties of the boundary layer, andi is 
later made specific by an appeal to experimental data. 

The mean velocity distribution u(y) through the turbulent boundary layer 
may be taken to depend on four local parameters 0, Ty, @, and mw; a fifth para- 
meter u, is of course implied as the value of wu when y = 0. The first two quan- 
tities, 6 and t,,, depend on a length coordinate ¢ which is not written explicitly, 
since it will be found that the Reynolds number R = u, é/v for the ideal bound- 
ary layer is adequately determined for a given profile by the assumptions 
already made in this and in the preceding section. 

The function of six variables 


(4, Y, Q, HM, O, Tw) = 9 (18) 


may be conveniently written without loss of generality as 
hem SAE ee 9 5 
w= 62%, 4), (19) 


in terms of a characteristic velocity u, and two characteristic lengths 6 and 
v/u,, where . 


Tw =e uz : (20) 


Suppose that equation (19) can be specialized for the region near the wall, 
on the ground of a certain similarity observed experimentally. Both the laminar 
sublayer profile and the adjacent turbulent mean velocity profile near the 
surface are found to be independent of the boundary layer thickness 6, so that’ 


Ur 


<= (2), y>O. (21) 


Suppose further that in the outer region of the boundary layer similarity is 
observed with respect to a coordinate system attached to the edge of the layer 
and moving with the free stream. It is not surprising, if the shearing stress in 
this region is maintained by a turbulent transport process, that the velocity 
defect is found to be independent of the viscosity of the fluid. In fact, turbulent 
shear flows with a free boundary, such as wakes and jets, apparently have in 
common with the boundary layer the property that near the free boundary the 
characteristic length, when considering the distribution of mean velocity mea- 
sured with respect to the outer fluid, is the width 6 of the shear layer. This 


observation is expressed by the velocity defect Jaw, which in the usual form 
is written 
Uy — 


gol (Slav ai ae (22) 


If now it is assumed that there exists a finite region in which both (21) and 
(22) are valid, it is easily shown ([12]; see also [13], [14]) that (19) in this region 


4 


2 
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must have the form 
U 
Ur 


==, n2* +9(3), (23) 


v 
or more specifically, in view of (21), 


U 1 Y Ur 


Ur aie, In v ¥ v(0) ; (24) 


Note that the existence of equation (23) is a sufficient condition to insure 
the existence of equation (22) since, if (23) is valid throughout the outer part 
of the boundary layer, 
Ur 


= 9(1) — o(5) —~—Int. (25) 


That the condition is also necessary is less commonly recognized. The 
equations (21) and (22), which are plausible on dimensional grounds and are 
supported by experimental evidence, and MILLIKAN’s demonstration of neces- 
sity, already cited, are together a less ambiguous basis for equation (24) than 


_ derivations involving specific assumptions about the mechanism of turbulent 
_ shear flow. In particular, the constant x is introduced during a separation of 


variables in a manner which shows clearly that this parameter should not 


_ depend on the surface shearing stress or on the boundary layer thickness, 


although it may depend on the wall geometry. 
Evaluation of equation (23) at the point y= 6 yields immediately the 


friction law 


uw 1 O Ur 
ne Mgrs tae p(1) , (26) 
which may be applied directly to pipe or channel flow on replacing the boundary 
layer thickness by the pipe radius or channel half width. For a boundary layer, 
however, it is more convenient to replace the boundary layer thickness 6 in 
equation (26) by the momentum thickness 6, with the result that 


ie es Ro + 9(1). (27) 


Ur x Cy — Cy (Ux/U) 


This expression follows immediately from equation (26) and the definitions 


é* fl (1 SF ) dy (28) 
and : 
6= gala) (29) 


on assuming that equation (22) is valid throughout the boundary layer. The 
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parameters c, and c, are then defined by 


Sak One A ys 30 
eer, [FG)43 ee 
0 
and ‘ 
a EE Ai eee = of”) gq» 31 
Due ea eer pews idl) ane (31) 
so that the form parameter 6*/6 is given by 
Om 1 ‘ 
= 3a 
OS eee ee ee 
Cree 


Furthermore, in view of the obvious relationship 


T Ur \2 
= 4% = : 33 

f q 2 es ; ba! 
it is clear that equation (27) is the analog of equation (17) in supplying, for 
the turbulent boundary layer, a connection between local friction coefficient 
and local momentum thickness Reynolds number which may be taken both as 
a definition and as a test for uniqueness. 

In order to complete the discussion of the ideal turbulent boundary layer, 
however, it is necessary to adopt a convention for defining a Reynolds number 
R. The present restricted idea of uniqueness implies a function 

dR 
2 aR = h( Re) (34) 
which may be integrated to yield a relationship 
R= H(R) , (35) 


on satisfying the initial condition RK = 0 when R, = 0. If the surface shearing 
stress is always positive, then R, is a monotonically increasing function of R, 
and the initial condition in question can obviously be satisfied. 

However, the integration can be carried out, and a precise definition of 
R achieved, only if the function / of equation (34) is specified for all positive 
values of Ry. But the function in question is defined experimentally, so that 
the extrapolation to arbitrarily small Reynolds numbers is not itself unique; 
in fact, it is not certain that the extrapolated function will be finite or even 
integrable at Ry = 0. For the present, the dependence of dRy/dR on Ry expressed 
by equation (34) will simply be identified with equation (27), and a discussion 
of the necessarily asymptotic correspondence between these expressions will 
be postponed until the experimental evidence has been reviewed. 
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The transcendental form of equation (27) makes direct integration difficult. 
However, in the absence of pressure gradient and roughness, the local friction 
4 coefficient and the momentum thickness for the ideal turbulent boundary layer 
may be taken to be related by the momentum integral equation, 


6 
ee to. sated d U u 
Sa ag ee a (1—<h) ay, eo) 
0 


E ee ne to equations (8), (29), and (33) above. In the present notation, using 
E22) and (26) in turn in (36), 


dR = e7* 90) 9% (mluz) [és — 7 Gat nes ay ee (37) 


Uy 
It is obvious from equation (27) that R, vanishes, implying R = 0, when 
_ U,/U, = C,/C,. With this initial condition, the definite integral of equation (37) 


a 


_ is the turbulent flat plate friction law for c,(R), 


Cp R= 20, e-* 9) eV Ber 1 (2 a és) = 4 : ( fe 1) «,| 


2 Cy 


4 2 e —% 91) e% (62/1) C 


As a check on this expression, it is readily verified that equation (9), 
R 
op R=2R, = | oR, 
0 


has equation (27) as a definite integral when the local friction coefficient is 
evaluated from (38). Further, it follows that the mean and local friction coeffi- 

_ cients are always finite, and are in fact equal in the limit Rk = 0, for the partic- 
ular friction law implied by functional similarity. Finally, anticipating numer- 
ical values to be obtained later, the last term on the right side of equation (38) 
is of order 40 c, compared to terms of order c, R, and may safely be neglected ; 
Cr Pp =o Cy e7 *9(1) e% (#x/¥r) [1 Ur (- ae *.) i 2 (= . (; mii 2) P (39) 


U, \x 64 x \Uy “ Cy 


The corresponding expression for the mean friction coefficient, from equa- 
tion (27), is 
—H# Hw (U,/UT C Ur 
¢; R=2R,=2,¢ ees rate ce (40) 
In summary, the two equations (39) and (40) above define the ideal local and 
mean turbulent friction coefficients for a flat plate in terms of four empirical 
velocity profile parameters x, p(1), c,, and cy. These parameters may sup- 
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Figure 1 
Direct Measurements of Local Turbulent Friction in Low-Speed Flow. 


posedly be determined by simultaneous local measurements of velocity profile 
and surface friction at a single Reynolds number, and extrapolation of the 
resulting friction formulae to arbitrary Reynolds numbers may then be under- 
taken with the same degree of confidence which is attached to the hypotheses 
of similarity underlying equations (21) and (22). | 


4, Experimental Literature?) 


While the turbulent boundary layer on a flat plate has been the subject of 
a great deal of experimental research, it is not easy to find reliable measurements 
of mean velocity profile and surface shearing stress made under identical experi- 
mental conditions. Before discussing the available measurements, therefore, 
it is convenient to provide a short catalog of several special techniques which 
have been used for determining local surface friction. 

The derivative dO/dx supposedly gives the surface shear directly, provided 
that the flow is uniform and two-dimensional. Since differentiation of experi- 
mental data is necessary, large errors may occur unless careful and complete 
measurements are made of the various velocity profiles. Alternatively, extra- 


polation toward the wall of the turbulent shear — gw’ v’ or of the sublayer 


1) The writer is indebted to Dr. G. B. ScnuBAvER and Mr. P. S. KteBanorr of the National 
Bureau of Standards, Washington; to Dipl.-Phys. W. Tixypmann of the ‘Max-Planck-Institut fiir 
Strémungsforschung’, G6ttingen; and to Dr. F. R. Hama of the State University of Iowa, Iowa 
City, for their courtesy in providing numerical data for some of the measurements cited here. 
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Figure 1 continued 


_ laminar shear yz 0u/0y implies unusual nicety of experimental technique. Direct 
measurement of the shearing force on an isolated element of the surface, 
although it requires a considerable effort in the development of special instru- 
mentation, appears to be the most reliable method available at present. Finally, 
a device based on a relationship between heat transfer and momentum transfer 
in the sublayer has been used successfully, but must be calibrated by comparison 
with a primary standard obtained by one of the more conventional methods of 
local friction measurement. 

The floating element technique is conspicuous among these methods for 
the consistency of the data obtained. Omitting a discussion of a meteorological 
application by SHEPPARD [15], three papers describing direct measurements of 
surface shearing stress on a flat plate in low-speed flow have appeared in the 
boundary layer literature. 

KeEmpPF [16] in 1929 measured the surface friction at several stations on the 
flat bottom of a ship model. The presence of a standing wave system on the 
free surface of the towing tank introduced into the friction coefficient data a 
scatter which is apparently systematic in terms of the Froude number, but 
random in terms of the Reynolds number; no attempt has been made to remove 
this scatter in the data considered here for the first four stations. KEMPr’s 
measurements were made at sufficiently large values of x so that the influence 
of the parameter %) should be negligible. A criticism of the data can, however, 
be based on the possibility of secondary flow in the boundary layer, since the 
ship hull in question, even though flat-bottomed, was relatively narrow in its 
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lateral dimension. A further criticism of the measurements at the largest 
Reynolds numbers has been expressed by FALKNER [17] on the ground that the 
length of run may have been insufficient to guarantee steady flow near the 
aft end of the ship model. 

ScHULTz-GRUNOW [18] later used similar methods to measure the friction 
on the floor of a wind tunnel. Great care was taken to minimize pressure varia- 
tions along and across the small gap provided for clearance around the floating 
element. The tunnel boundary layer was removed at an upstream slot, whose 
downstream lip formed the leading edge of the plate under investigation. 
Transition to turbulent flow occurred immediately downstream of the curved 
leading edge; however, there is some doubt about the magnitude of the para- 
meter x», especially since the Reynolds number was apparently varied by 
changing the velocity at a fixed station along the plate. 

More recently, DHAWAN [4] has obtained further data on the turbulent skin 
friction on a flat plate by means of the floating element technique. DHAWAN 
was able to induce transition by means of leading edge curvature, and to con- 
firm the existence of turbulent flow over most of the working surface by means 
of an axial impact pressure traverse near the plate surface. 

The experimental data of KEmpr, SCHULTZ-GRUNOW, and DHAWAN are col- 
lected in Figure 1. The measurements appear to define a single curve over the 
range of Reynolds numbers from 2 x 10° to 4 x 108, with a scatter so small 
as to be almost unprecedented in boundary layer research. In fact, it may be 
admitted that the present re-examination of the turbulent boundary layer 
problem was to a large extent inspired by the remarkable consistency of these 
local friction data. 

The existence of a unique relationship between the ordinate and the slope 
of the momentum thickness distribution, O(«), has already been proposed as a 
criterion for uniqueness in the turbulent flat plate boundary layer. The imme- 
diate question is one of demonstrating that such a relationship exists for the 
experimental data, without necessarily identifying the resulting function for 
example with equation (27) above. In view of the wide variations which occur 
from one experiment to another in various factors which influence boundary 
layer transition, it is not obvious a priori that a state is ever reached in which 
the dependence of the turbulent boundary layer on its early history is no longer 
measurable in terms of the local mean properties of the flow. 

With the restricted definition of uniqueness considered here, it is sufficient 
as a beginning to consider the momentum thickness distribution only. The 
measurements which are most useful for this purpose were made by WIEG- 
HARDT in 1943, as part of an extensive research program carried out since 
1939 in a boundary layer tunnel at the ‘Institut fiir Str6mungsforschung’ at 
Géttingen. While this program has placed considerable emphasis on the deter- 
mination of the drag of surface irregularities, including uniform roughness, and 


SE NTE ee ee ee ee 


Vol. V, 1954 The Problem of the Turbulent Boundary Layer 193 


on the investigation of turbulent boundary layer flow in a pressure gradient, 


the case of a flat plate boundary layer has also been investigated in some detail. 

Continuity in these German experiments is provided by WIEGHARDT’S velo- 
city profile measurements in uniform flow, carried out at free stream velocities 
of 17-8 and 33-0 m/s respectively. In addition, ScHutz-GRUNOW’s experiments 


at 19-4 m/s were made in the same channel. Of these data, WIEGHARDT’S mea- 
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Figure 2 
Velocity Profile Friction Measurements in Low-Speed Flow. 


_ surements at 17-8 m/s, by virtue of early transition and almost immediate coin- 
_ cidence for the actual and apparent origins, are in very good agreement with 
_ the properties of the ideal turbulent boundary layer. The local friction coeffi- 


cients for this series, obtained by differentiating the experimental momentum 
thickness distribution, are plotted against Re in Figure 2. 

Some recent measurements in a flat plate boundary layer at the National 
Bureau of Standards, reported by KLEBANOFF and DIEHL [19], were under- 
taken in part to test the hypothesis of uniqueness for thick turbulent boundary 
layers produced for example by initial roughness. Figure 2 shows the relation- 
ship between 2 dO/dx and Ro observed in experiments with natural transition 
at a free stream velocity of 108 ft/s, and with initial sand roughness at 35, 55, 
and 108 ft/s. Except in the last instance, the most downstream data are in good 
agreement with WIEGHARDT’S curve ; the agreement is in fact better than might 
be expected considering the uncertainty in the local friction coefficients. Exam- 
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ination of the experimental velocity profiles shows that the present criterion 
for uniqueness is equivalent to the test proposed by KLEBANOFF and DIgHL, of 
similarity in the velocity distribution, and serves equally well to identify the 
profiles which should be studied in detail. 
Finally, Hama [20] has published some velocity profile measurements in 
turbulent boundary layers, and his experimental values for the distribution of 
@(x) have been used to obtain the local friction coefficients plotted in Figure 2. 
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Figure 3 
The Law of the Wall. 


Several velocity profiles which should be representative of the fully devel- 
oped turbulent boundary layer are shown in Figure 3, in the coordinate system 
of equation (21), in order to exhibit the occurence of a logarithmic relationship 
in the profile near the surface. According to equation (24), the value of the 
parameter x is determined by the slope of the straight line in the figure, and 
y(0) by the intercept at y u,/y = 1. 

In establishing a value for the boundary layer thickness 6, needed in deter- 
mining the parameter c, from equation (30), it has been found convenient to 
plot (1 — u/u,)** against y and to fit a straight line to the outer part of the 
profile in this coordinate system. The thickness 6 is defined as the intercept on 
the y-axis. This convention is neither new nor significant; it is suggested for 
example by the mixing length theory of PRANDTL when the shearing stress 
varies linearly with y and the mixing length is constant, but it is justified at 


present only by its convenience as a well-defined method for specifying the 
boundary Jayer-thickness. 
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Figure 4 shows the profiles of Figure 3 in the coordinate system required by 
the velocity defect law, equation (22). Again fitting a straight line with a slope 
determined by x to the logarithmic region, the intercept at y = 6 = 6* Uy|Cy Uy 
is p(1) — (0), according to equations (24), (26), and (30). 

This analysis of the turbulent velocity profiles has been carried out for all 
the data available to the writer, using a value of x of 0-40 and determining 
values for (0), (1), c,, and cg, the latter from equation (31). For WIEGHARDT’S 
measurements at 33-0 m/s a modified procedure was followed, the local friction 
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Figure 4 
The Velocity Defect Law. 


coefficients being obtained from the expression (27) using the measured values 
of Ro and preliminary estimates for c,, cz, and g(1). The collected values for 
several of the profile parameters are plotted against Re in Figure 5, and show 
a satisfactory degree of consistency. 

It is necessary finally to choose particular values for the parameters x, (1), 
c,, and Cy, in order to complete the numerical analysis of the turbulent bound- 
ary layer. Recognizing the asymptotic nature of the treatment by emphasizing 
the data at the larger Reynolds numbers, then approximately 


x= 0-400, 9(0)=5-10, o(1)=7-90, c= 405, c= 29.0. 


These parameters have been used to evaluate R, Ro, Cr, and 6*/@ as func- 
tions of c, from equations (39), (40), and (32), and the resulting data are collected 
in Table I. The values tabulated in the first and second columns have been 
plotted in Figure 1, and those from the first and fourth columns in Figure 2. 
The values tabulated in the fourth and fifth columns are compared with the 
experimental values of the shape parameter 6*/@ in Figure 6, the agreement 
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Figure 5 
Velocity Profile Parameters. 
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Figure 6 
Boundary Layer Shape Parameter. 


here of course being expected since these experimental data were used to 
determine Cy. 

It may be pointed out that the curves in Figures 1 and 2 are in no sense fitted 
to the measurements in question, since the ideal curves are based on detailed 
information about the velocity profiles, while the measurements are not. 
Furthermore, although experimental values for uw, were used to eliminate the 
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Table I 


Suggested Properties of the Ideal Turbulent Boundary Layer in a Fluid of Constant Density 


BA te. 
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8-0 1.017 1-032 _,| 5-25 1.828 0.0632 
7.8 1-143 1-005 “on a tere 1-809 0.0624 
76 1-291 9.78 6-31 1-790 0-0616 
74 1.463 9.52 | 6-96 } x10 1-772 0-0608 
7:2 1-666 9.25 7-70 1-753 0-0600 
7-0 1-906 8-98 8-55 1-735 0.0592 
6-8 2.193 8-70 9.54 77 0-0583 
6-6 2-536 t «104 | 8.43 1-069 1-699 0-0575 
6-4 2-950 8-16 1-204 1-681 0-0566 
6-2 3-45 7-89 1-362 1-663 0-0557 
6-0 4.07 7-62 1-551 1.645 0-0548 
5-8 4.84 7.35 1-776 1-628 0-0539 
5-6 5-79 7.07 2-047 1-610 0-0529 
5-4 7-01 6-80 2383 TS ais 1-593 0-0520 
5-2 8-49 | 6-53 2771 ie 1-575 0-0510 
5.0 1-050 ] 6-26 3-29 1-558 0-0500 
4.8 1-310 5-99 3-92 1-540 0-0490 
4-6 | 19-8 | 1-653 | 5-72 4.73 1-523 0-0480 
4-4 Sh Mia eg eer ae 1-506 0-0469 
4.2 2.748 t 5.19 743 1-488 0.0458 
4.0 3-64 4.92 8-95 1-471 0.0447 
3-8 4.91 4.66 1-143 1-453 0-0436 
3-6 6-78 4.39 1-487 1-436 0.0424 
3.4 9.72 4.13 2.004 1-419 0-0412 
3-2 1-398 3-87 2-700 } «103 1-401 0-0400 
3-0 2-107 3-60 3-80 1-383 0-0387 
2.8 330 $ x108 | 3.35 5-53 1-366 0-0374 
2-6 5-43 | 3-09 8-38 1-348 0-0361 
2-4 9.43 2-835 1-337 1-330 0.0346 
2.2 1-755 l 2.581 2.265 | oe 1-311 0-0332 
2:0 3-55 | x10? | 2.331 4.14 1-293 0-0316 
1-8 8-00 2.083 8-33 1.274 0-0300 
a ei ey 504 Baan ea asa oy |\at 268 
1.2 2.583 x 109 1-354 1.748 | 408 1-213 0-0245 
1.0 | 1.518101 | 1.117 8-48 ay 1-191 0.0224 


influence of Reynolds number in the analysis of the profiles, it was not neces- 
sary to know the magnitude of the Reynolds numbers in question. In fact, the 
discrepancies in Figure 1 are perhaps connected with the assumption that the 
real and apparent origins coincide for the experimental boundary layers; com- 
pare Figure 2, where local quantities only are involved. 

Finally, WrecHarpt1’s velocity profile data at a free stream velocity of 
33-0 m/s are shown in detail in Figure 7 in order to illustrate the consistency 
of these particular measurements in terms of the law of the wall, equation (21), 
and the velocity defect law, equation (22). 
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Figure 7a 


Gottingen Measurements at 33-0 m/s. 
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Figure 7b 


Gottingen Measurements at 33-0 m/s. 
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5. Discussion 


Most of the ideas which appear in the present study are not new. For 
example, equation (23) was postulated on empirical grounds, and equation (39) 
derived as a consequence, by voN KArMAN in 1932 [21], and ScHuLTz-GRUNOW 


_ in 1939 carried an analysis of his experimental data almost to the point where 


the numerical results of the present paper could have been obtained. More 


_ recent papers by Rortta [22], BAINEs [23], and LANDWEBER [24] have discussed 


the representation of the turbulent boundary layer in terms of the law of the 


wall and the velocity defect law, but have not emphasized the physical inter- 


pretation of these formulae. The most important single element in the discussion, 
and one deserving a place among major contributions to the boundary layer 


_ literature, is the demonstration by MILLIKAN [12] that the two similarity laws 
_ together imply a logarithmic region in the velocity profile. 


In reviewing the asymptotic nature of the turbulent boundary layer prob- 
lem against the background of the similarity laws, it is obvious that extrapola- 
tion of a local friction law based on functional similarity to indefinitely small 
Reynolds numbers is permissible only in a restricted sense. For example, the 
parameter c, is evaluated from equation (30) ; however, the integral in the latter 
expression neglects the dependence of the sublayer flow on Reynolds number, 
and so cannot be uniformly valid for all the experimental data considered here. 
The contribution of the sublayer flow to the integral defining the displacement 
thickness increases with decreasing Reynolds number, and the approximation 
involved in identifying the measured quantity u, 6*/u, 6 with c, for small 


_ values of R, is apparent in Figure 5. 


Although this fault is easily corrected at the cost of slightly more compli- 
cated formulae [24], it is evident that the two similarity laws cited here, and 


- therefore equation (27), must inevitably fail at sufficiently small values of the 


Reynolds number. The region of common validity is marked by a logarithmic 


velocity profile, which is known to obtain only for values of y w,/y greater than 


say 50. However, Figure 4 shows that the logarithmic law can describe the fully 
developed velocity profile only within the inner 20 percent of the boundary 
layer in any case, so that the smallest value of 6 u,/y which satisfies both 
criteria is about 250; the corresponding value of u, 0/y is about 700. Consulting 


the experiments, the approximate value of the abscissa in Figure 5, below 
_ which the parameters g(0) and (1) appear to be no longer sensibly independent 


of Reynolds number, is “, O/y = 1500. 
It is clear, from the evidence considered here, that the velocity defect law 


| is of primary importance in determining the asymptotic properties of the ideal 
- turbulent boundary layer for Reynolds numbers of technical interest. The prob- 


able physical significance of this fact may be discussed in terms of the charac- 
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teristic scales which were introduced early in the present organization of the 
boundary layer problem. 

Turbulent fluid motion as observed experimentally is characterized at one 
extreme by the presence of large eddies, which are responsible for accepting 
energy from a prime mover, and at the other extreme by the presence of small 
eddies, in which viscous dissipation accompanies an approach to isotropy; the 
corresponding characteristic scales are distinct provided that the Reynolds 
number is sufficiently large. Although study of fluctuating phenomena in a tur- 
bulent boundary layer in terms of spectra and scale has not yet reached a state 
where it is possible always to distinguish cause and effect, it may be supposed 
that the characteristic lengths 6 and »/w, which appear in the two velocity laws 
are in fact useful measures of a macroscale and microscale for the turbulent flow. 

Measurements of the space correlation across a turbulent channel flow by 
LAUFER [25], across a free jet by CoRRSIN and UBEROI [26], and across a bound- 
ary layer by KLEBANOFF [9], together with the observation that a velocity 
defect law of the kind considered here holds at least approximately for such 
flows, are convincing evidence that large eddies are important in the general 
turbulent shear flow whether the boundaries are free or solid or both. Since the 
velocity defect is also the momentum defect for a fluid of constant density, it 
is to be expected that the defect law and the equations of mean motion to- 
gether must determine the shearing stress not only at the wall but throughout 
the fully turbulent shear flow. Furthermore, both the velocity defect and the 
turbulent shearing stress have as a characteristic length the thickness 6, which 
represents the large eddy structure and thus in some sense the history of the 
flow; the same characteristic length is found for the phenomenon of intermit- 
tent turbulence at the free edge of the shear flow [9], [11]. These observations 
perhaps explain some deficiencies in mixing hypotheses patterned after a mole- 
cular model, since the present assertion is first that the actual transport process 
is connected with the recent history of the motion as well as with purely local 
conditions, and second that momentum transport is accomplished by the 
largest eddies in the flow. 

The physical significance of the law of the wall is not yet clearly recognized, 
except as a consequence of the constraint applied to the turbulent motion by the 
condition of vanishing velocity at the fixed boundary!). The law of the wall 
appears to be related primarily to local conditions and thus to the small-scale 
turbulent structure, since the dimensionless velocity profile is found experi- 
mentally to be insensitive to variations in ambient pressure distribution [27] 
or to changes in the nature of the other boundaries of the flow. In this connec- 
tion it is important to note that, although the law of the wall can be studied 
onlyin the presence of a finite shearing stress at the surface, the law is inherently 
incapable of providing any estimate of magnitude for the shear. 

1) See The Law of the Wall in Turbulent Shear Flow, to be published separately. 
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The present study suggests strongly that the assumptions which lead to 
the logarithmic velocity distribution are in effect equivalent to the assumption 
that there is a region in the flow where the mean velocity is dependent on the 
small eddy structure and thus on the dissipation, while in the same region the 
momentum defect and shearing stress are dependent on the large eddy struc- 
ture. Thus a search for a relationship between shear and velocity is in this 
sense a search for a relationship between the turbulent structures at large and 
small wave numbers, and so might be undertaken using methods which have 
been found useful in the analysis of isotropic turbulence. 

The widespread current revival of interest in the theoretical and experimen- 


tal aspects of the boundary layer problem can hardly fail to lead to important 


discoveries, both in a knowledge of the structure of turbulent shear flows and in 
the formulation of an adequate theory. In the meantime, the writer hopes that 
the discussion given here, of the more rudimentary properties of one such flow, 


_ may serve to allow the identification of a proper environment for further experi- 


mental research, and also may serve to provide a point of departure for study 
of the effects of pressure gradient, compressibility, and roughness. 
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Summary 


Existing measurements of low-speed turbulent surface friction on a flat plate, 


in the absence of pressure gradient and roughness, are shown to be consistent 
with a simple analysis based on functional similarity in the velocity profile. In 
particular, the fully developed turbulent boundary layer is found to be unique 
within the accuracy of the experimental data, with uniqueness defined as the 
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existence of a definite correspondence between local friction coefficient and mo- 
mentum thickness Reynolds number. The relationships known as the law of the 
wall and the velocity defect law are found to describe the turbulent velocity pro- 
files accurately for a considerable range of Reynolds numbers, and an effort is 
made to clarify the physical significance of these formulae. Finally, the proper 
definition of a length Reynolds number is discussed in terms of the asymptotic 
local properties of the ideal boundary layer, and numerical values for ideal mean 
and local friction coefficients are tabulated against Reynolds numbers based on 


- momentum thickness and on distance from the leading edge. 


Zusammenfassung 


Es wird gezeigt, dass vorhandene Messungen der turbulenten Wandschub- 
spannung an der glatten ebenen Platte in inkompressibler Strémung ohne Druck- 
gradient durch eine einfache Berechnung in Ubereinstimmung gebracht werden 
kénnen. Die Rechnung beruht auf einer funktionellen Ahnlichkeit der Geschwin- 
digkeitsverteilung. Es wird im besonderen gefunden, dass die vollentwickelte 
turbulente Grenzschicht innerhalb der Messgenauigkeit einem eindeutigen Zu- 
sammenhang zwischen dem Ortlichen Reibungskoeffizienten und der Reynolds- 
schen Zahl, bezogen auf die Impulsdicke, folgt. Die Beziehungen, die als Wand- 
gesetz und Mittengesetz bekannt sind, beschreiben die Geschwindigkeitsverteilung 
genau innerhalb eines erheblichen Bereiches Reynoldsscher Zahlen, und es wird 
versucht, den physikalischen Inhalt dieser Gesetzmassigkeiten zu vertiefen. Ab- 
schliessend wird eine zweckmassige Definition der auf Plattenlange bezogenen 
Reynoldsschen Zahl diskutiert, die auf dem asymptotischen 6rtlichen Zustand 
der idealen Grenzschicht beruht. Rechenwerte der idealen, mittleren und 6rtlichen 
Reibungskoeffizienten, bezogen auf beide obigen Definitionen der Reynoldsschen 
Zahl, werden tabelliert. 


(Received: August 17, 1953.) 


Untersuchungen tber Elektronenstromungen 


Von JOHANNES MULLER, Berlin 


Gewisse Eigenschaften im Verhalten von Elektronenstré6mungen auf klas- 
sischer Grundlage werden aus Matrizendarstellungen gewonnen. 

Die Bedingungen, wie sie durch die Erhaltungssatze fiir Energie und Ladung 
bei Elektronenstrémungen auftreten, werden fiir den Fall kleiner Abweichun- 


-gen vom stationaren Zustand ermittelt. Die Anordnungen sollen sich also linear 


verhalten, das heisst, die Wirkung einer Anregung soll ohne gegenseitige St6- 
rung mit der Wirkung einer anderen Anregung tiberlagert werden kénnen. Die 
Beeinflussung der Elektronen erfolgt nur iiber das elektromagnetische Feld. 
Wechselwirkungen zwischen Elektronen, die zu einer erheblichen Anderung 
ihrer Geschwindigkeit fithren — das sind Zusammenstésse im Sinne der kine- 
tischen Gastheorie —, werden vernachlissigt. Die Gesetzmassigkeiten werden 
fiir Strémungen mit einer und mehreren Teilchengeschwindigkeiten sowie mit 
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Geschwindigkeitsspektrum aufgestellt. Fiir ein Exponentialspektrum, wie es 
bei der Emission der Elektronen von einer Gliihkathode entsteht, ergeben sich 
zwanglos Ausdriicke, wie sie unmittelbar fiir die Theorie der Schwankungen 
von Bedeutung sind, was mit dem energetischen Charakter des Rauschens 
zusammenhangt. Durch ein hohes Magnetfeld wird der raumliche Verlauf des 
Weges der Elektronen vorgeschrieben. Dadurch ist direkt nur das elektrische 
Feld in Richtung des Weges fiir die Bewegungsvorgange verantwortlich. 


Aufstellung einer Leistungshilanz 


Gegeben sei eine Strémung als Strahl mit einer Teilchengeschwindigkeit an 
jedem Ort, deren Elektronen die Wegstrecke von O iiber s, s bis s,, zurtick- 
legen (Figur 1). Die statischen Laufzeiten betragen dabei 0, t(s), t(s) und t(s,,), 
der Strahlgleichstrom 7. Den statischen elektrischen Feldern ist langs der Weg- 
strecke ein schwaches elektrisches Wechselfeld iiberlagert. Wegen des starken 
magnetischen Fiihrungsfeldes haben Wechselmagnetfelder keinen Einfluss. Zu- 
nachst werde nur ein Wechselfeld €(s) langs der Wegstrecke ds am Orte 5 
angenommen. Dadurch erhalten die Elektronen eine infinitesimale Volt- 
geschwindigkeitsmodulation (= %-Modulation) von &(s) ds exp jm ¢t zur Zeit ¢ 


Ws) Lls) 04, C4) 


Figur 1 


Richtungen der positiven Halbwellen von Wechsel-Feldstarke &(s), -Voltgeschwindigkeit @ und 
-Konvektionsstrom £ sowie die Bezeichnungen der Wegstrecken und Laufzeiten. Bewegungsrichtung 
der Elektronen von O biss,,. 


an 5. Zur Zeit ¢ + t(s) — t(s) werden diese Elektronen mit dem unverandert 
gebliebenen Zuwachs an ¥-Modulation am Punkt s ankommen (Erhaltung 
der Energie). Die statischen Laufzeiten sind eine eindeutige Funktion der An- 
fangsenergie oder der statischen Anfangsvoltgeschwindigkeit U, der Elektronen, 
wenn zunichst keine Kreisbewegung mit Zyklotron-Frequenz senkrecht zu dem 
magnetischen Fiihrungsfeld vorhanden ist. Fiir die Erhdhung der Laufzeit der 
zur Zeit ¢ an s ankommenden Elektronen ergibt sich dann: 

oes E(s) ds ef@t-eles)—a9)) 
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Mit dem Gleichstrom 7 multipliziert, ergibt dies den Uberschuss der Elek- 
_ tronenraumladung zwischen s und s zur Zeit ¢. Aus der zeitlichen Anderung 
_ dieser Ladung erhalt man den infinitesimalen Konvektionswechselstrom an s 
zur Zeit ¢ (Erhaltung der Ladung). Das Vorzeichen hangt mit den gewahlten 
Richtungen (Figur 1) zusammen. 


dQ(s) =e sii (s) | &(s) ds j o e) tJ [x(s) —2(8)] 
; ‘ 1) 
; —j [x(s) —7(8)] ( 
SOA dee ae 


é@ 


Da es sich im folgenden stets um periodische Vorgange handelt, wird die 
Zeitfunktion weggelassen. Wenn auf der gesamten Wegstrecke ein Wechselfeld 
besteht, so erhalt man durch Integration die durch Feldmodulation entstehende 
Konvektionsstrommodulation (= £-Modulation), da bei kleinen Abweichungen 
die Einzelwirkungen sich addieren, aus: 


—j w [t(s)—1(5)] 
9-4/4 i ills cS (2) 


U, 


Se Rn Ce a Ne 


ll le it ol 


OP Aw 2 ® 


Die Riickwirkungsfelder infolge der Abstrahlung elektromagnetischer Ener- 
gie durch die veranderlichen Bewegungen der Elektronen in ihrer Gesamtheit 
(die Strahlung des einzelnen Elektrons wird vernachlassigt) sind in der angenom- 
menen Wechselfeldverteilung €(s) mit enthalten. Wird (1) auf die Anfangs- 
¥-Modulation angewendet [€(s) ds wird durch B(0) ersetzt, t(s) = 0], so ergibt 
sich dieser Anteil der £(s)-Modulation zu: 


dei ers) 


aU, (3) 


LNs) = 1 B(0) — 


Damit erhalt man fiir die gesamte 2-Modulation an s, wenn noch die Wirkung 
der Anfangs-£-Modulation hinzugenommen wird: 


s) 


-y de~i oles) -718)] de-i” ae 
2) =1 | (5) ds 2 + 1B) “+ eI.) 


é 


Wenn die Feldverteilung erzwungen wird, wie bei der Ableitung der Lei- 
stungsbilanz, so kann die Differentiation nach U, vor das Integral gesetzt 
werden. Aus (4): 


Ris) =1 a; / &(s) ds e771 eles)— ENE B(0) e “Fax8) | 4 (0) 0) en fet) (5) 
0 
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Der Ausdruck in der eckigen Klammer ist nichts anderes als die 8-Modulation 


an S. 


B(s) = [ &) ds e— 1 [x(s)—1(5)] el B(0) e771 ers) : (6) 

; ‘ 
Wenn das Wechselfeld nicht erzwungen wird, so gilt die folgende Gleichung (7) _ 
Br dB(s) —j wr(s) | 
Lis) = 1 ay + 2(0) (7) 


nicht bei endlicher Stromstarke, da eine Riickwirkung auf die Feldverteilung- 
wegen der sich verandernden Raumladungsverteilung beim Ubergang zu einer 


um dU, erhéhten statischen Anfangsvoltgeschwindigkeit fiir die gesamte Str6- 
mung eintritt. Dieser Einfluss geht gegen Null, wenn der Gleichstrom der 
Strémung gegen Null geht. Allgemein gilt: 


Qs) = i SF) 


+ (0) e707) | (8) 
Durch Umformung kann aus (6) nach dem Wechselfeld €(s) aufgeloést werden: 


?: a[ B(s) ef O*)] 
j @t(s) 
é Sopsgarr rs. (9) 


Die durch das elektromagnetische Wechselfeld der Strémung zugefiihrte Lei- 
stung betragt: 


E(s) = 


sm 


N= / €(s) ds 2*(s) 2). (10) 
0 
€(s) aus (9), L(s) aus (7) eingesetzt und partiell integriert, ergibt: 
N = B(5;p) L# (5p) — B(0) L*(0) — Fw / B(s) Be(s) i 479 as. (1) 
é é 


Diese Gleichung besagt, dass die durch das Wechselfeld zugefiihrte Wirklei- 
stung [Realteil von §t] und die am Anfang der Strémung zugefiihrte Wirklei- 
stung [Realteil von B(0) L*(0)] zusammen die abgegebene Wirkleistung am 
Ende der Strémung [Realteil von B(s,,) 2*(s,,)] ergibt. Da die Voltgeschwindig- 
keit stets ein Mass fiir die Energie der Elektronen ist, gilt die Beziehung auch 
bei relativistischer Massenanderung. Das dritte Glied auf der rechten Seite von 
(11) stellt eine reine Blindleistung dar. Werden die Wechselgeschwindigkeiten 
eingefiihrt, so ergibt sich nichtrelativistisch dieser Anteil zu: 


7 w | mo v(s) v*(s) dn. 
0 


| 
| 
| 


: 
| 
: 


:| 


; 
; 
; 


: 
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_ Dies ist die Blindleistung der Energie der Wechselbewegung und entspricht 


vorzeichenmassig einer Induktivitat. du sind die Anzahl der Elektronen mit 
Wechselgeschwindigkeiten p(s) und Ruhemasse my. 

Es soll nunmehr der Einfluss einer Bewegung senkrecht zum magnetischen 
Fuhrungsfeld untersucht werden. Zu diesem Zweck wird das Quadrat der Be- 
wegungsgrosse senkrecht'zu dem Fiihrungsfeld 


Prt: ahi 2) = 1p ex) 


nach der Zeit differenziert. e, ist der Einheitsvektor in Richtung der magneti- 


_ schen Induktion. Unter Benutzung des Bewegungsgesetzes 


apn) PN EB 


dt m 


(bei Elektronen ist die Elementarladung e positiv einzusetzen) ergibt sich, 
wenn die ersten rdumlichen Ableitungen der Induktion stetig sind: 

dB 
aq | 


dp? p? 
T+ ——2e(p, E,)-2-4. 


U) ° 


 Dabei ist dB/dq, die Anderung der Induktion in Richtung des augenblicklich 


bestehenden Impulses # , . Wenn die Intensitat des Magnetfeldes gegen oo geht, 
verbleibt senkrecht zum Magnetfeld nur eine Kreisbewegung mit Zyklotron- 
Frequenz. Damit wird der zeitliche Mittelwert von (6, E,) Null und von 
4B/dq,, wegen der Quellenfreiheit der Induktion, —dB/2 ds. Es folgt: 


Pirie indy CB. oa PA, 
B 


: Por 
1 at aE Re jee 


0 


= const — 


Dies bedeutet, dass die Bewegungsgrésse #, keine Modulation durch die Wech- 
selfelder erfahrt. #% /B ist proportional dem Drall der Kreisbewegung, der mit 
B + oo gegen Null geht. Es kann jetzt angegeben werden, wie die Bildung von 


_dt/dU, in (1) bei Vorhandensein einer Querbewegung vorzunehmen ist. Die 


Bewegungsgrésse senkrecht zum Fiihrungsfeld #) , am Anfang wird dabei unver- 
andert gelassen. Wenn die statischen elektrischen und magnetischen Felder in 
ihrer Richtung nicht zusammenfallen, erfordert eine Realisierung derartiger 
Strémungen allerdings erhebliche Magnetfeldstarken. Die grundsatzliche Még- 


-_lichkeit solcher gefiihrter Str6mungen ergibt die Berechtigung, das Verhalten von 


Str6mungen mit Steuerungen in Fortschreitungsrichtung in allgemeiner Weise 
zu behandeln. Die Ergebnisse kénnen sofort auf Falle angewendet werden, bei 
denen kein oder nur ein schwaches Magnetfeld besteht, wenn die Wirkungen 
von Querbewegungen (bedingt durch seitliche thermische Anfangsenergien 
oder durch Quersteuerung langs der Stromung) vernachlassigt werden konnen. 

Die Wirkleistungsbilanz gilt auch fiir Elektronenwege mit Umkehrstelle. 
Wahrend bei Wegen ohne Umkehrstelle die Integrationsgrenzen fest sind, ver- 
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andert sich der Integrationsbereich durch die Ortsanderung der Umkehrstelle. 


* 


Ein Teil des Weges, der die Umkehrstelle enthalt, werde fiir sich betrachtet. _ 


Dieser Teil soll einmal so gross sein, dass am Anfang und Ende dieser Teilstrecke 
die Wechselgréssen klein gegeniiber den statischen Werten sind, andrerseits soll 
er so klein sein, dass trotz der zeitlich veranderlichen Vorgange die Verhalt- 
nisse in diesem Gebiet als statisch angesehen werden kénnen, so dass die Wirk- 
leistung am Anfang und Ende dieser kurzen Wegstrecke die gleiche ist. Die 
Wirkleistungsbilanz wird damit durch derartige Umkehrstellen nicht verancert. 


I. Str6mungen ohne Anfangssteuerung 


Die Wirkleistungsbilanz kann auf beliebig viele Einzelstromungen erweitert 
werden. In einem elektromagnetischen Hohlraum!) ohne elektromagnetische 
Verluste sollen sich m Einzelstrémungen befinden, deren Wege an den Begren- 
zungen des Hohlraums beginnen und enden. Damit wird den Strémungen ins- 
gesamt weder elektromagnetische Wirkleistung entnommen noch zugefihrt, und 
die Wirkleistung am Anfang der Wege ist gleich der Wirkleistung am Ende der 
Wege. Seien 2, und B, bzw. 2” und B™ die Wechselkonvektionsstréme und 
die Wechselvoltgeschwindigkeiten am Anfang bzw. am Ende des yw-ten Weges, 
so gilt bei m Einzelstromungen: 


Sy Qo Qe 1 Ql) Qylu)e ne B, ee oe sm B* : (12) 
wal pw=1 

Auf beiden Seiten steht die doppelte Wirkleistung. Bei kleinen Abweichungen 
besteht ein linearer Zusammenhang zwischen den End- und Anfangsmodula- 
tionen, darstellbar durch 2 m Gleichungen. Die Anfangsmodulationen sind bei 
Str6mungen ohne Anfangssteuerung voneinander unabhiangig. Die Endmodu- 


lationen des w-ten Weges als lineare Funktion von den Anfangsmodulationen 
sind angeschrieben?) : 


gH) —Q 24+ 2 B,+---+2 


- ~ 
HQ, HB, 


Re — BY Ati == Y 


Ah 
HQ 


HQ Lin “8 Lag, Bin : | 
Bite + Big Qin + Bu os, Bin 
(4 eh 2S ott) 


1) Der Hohlraum dient dabei nur dazu, den Energieaustausch auf definierte Stellen durch 
Vermeidung von Abstrahlung zu lokalisieren, und ist bei der Anwendung auf Anordnungen, bei 
denen dies ohne weiteres berechtigt ist, nicht wesentlich. 

) Derartige Koeffizienten sind bei ebener Strémung von F. B. LLEWELLYN, Bell Syst. Tech. J- 
14, 652 (1935); F. B. LleweLiyn und L. C. Perrerson, Proc. I. R. E. 32, 144 (1944), und bei ebenen 
Strémungen mit Ionen von W.O. ScHumMANN, Naturwissenschaften 31, 140 (1943), berechnet 
worden. H. W. Kénic, Arch. elektr. Ubertr., Nov. 1952, 445, findet bei ebener Stroémung mit einer 
Tragergeschwindigkeit eine einfache Beziehung zwischen den Koeffizienten. 

Anmerkung bei der Korrektur: Durch die Verwendung der Voltgeschwindigkeit anstatt der 
Geschwindigkeit werden die Beziehungen formal vollstandig analog denjenigen bei einem verlust- 
freien System von m Leitungen bzw. einem verlustlosen (2 m)-Pol. 


(13) 


HQ, 


a 


To te ee Se ee ST ee SE ee Re Ne eee ee . | 
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oder in Matrix-Form: 


: | 
und transponiert | : 
(2) BY...) = (2, B,...) A. (14b) 
Die (2 m)-reihige quadratische Matrix Y lautet: 
Lig. Lie, Zio, Lig 
Bs = Bi, &, By 9, B, a 
»)f a x, Q, se B, ae Q, sue B, alee . (15) 
Boo Bo. Boo Bo, 
(14b) wird in der “cee Weise erweitert: 
: oh || RO) * BV, * 
=(, BY, ee | QT’ | QQ = (2, G, ee 3 | Duy . (16) 


ae 


| | 
| Pe 


Auf der linken Seite steht dann die doppelte Wirkleistung der Endmodulatio- 
nen, die mittels (12) durch die doppelte Wirkleistung der Anfangsmodulationen 
ersetzt wird. (14a) werde links mit der Vertauschungsmatrix © multipliziert. 
Dadurch wird die 2. mit der 1., die 4. mit der 3. Zeile usf. vertauscht. Multipli- 
kation der %{-Matrix rechts mit © bewirkt die entsprechende Spaltenvertau- 
schung. € € ergibt keine Vertauschung. Aus (14a) entsteht: 


ay yey ey sh 
CACC B,)=C a al =C€AC] & | ; 
| ih sas | 


t 


Fiir € AC wird  gesetzt. W geht aus W durch Vertauschung der 2. Zeile und 


Spalte mit der 1. Zeile und Spalte usf. hervor. Es gilt W = W und Y-4 = w-1. 
Die Reihenfolge von Vertauschung und Transponierung bzw. Bildung der in- 


versen Matrix ist gleichgiiltig. Die Determinanten von Y& und Q sind gleich. Der 
konjugiert-komplexe Wert von (17) gebildet und in (16) eingesetzt, ergibt: 


B, ||* | B, || 
(2, B,...) WM W* || 2, |) = (2, Bi...) Qa |. (18) 
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Die Gleichung ist fiir beliebige Anfangsmodulationen erfillt, wenn: 


-O/ Y* = E und transponiert WU W/* = E. 


(19) 


W’* ist damit die zu Y inverse Matrix Y-1, und wenn die Elemente von Y die 
Bedingung (19) erfiillen, stellt {&{ eine Transformation dar, bei der die Wirk-_ 
leistung unverdndert bleibt. E ist die Einheitsmatrix. Da die Multiplikation in — 
(19) vertauschbar ist, ergeben sich sowohl durch Zeilen-Spalten- wie durch — 
Spalten-Zeilen-Multiplikation Beziehungen zwischen den Elementen. Die glei- 
chen Ergebnisse erhalt man durch Zeilen-Zeilen- und Spalten-Spalten-Multiph- 


kation von & mit W*. Spalten-Spalten-Multiplikation liefert Gesetzmassigkei- 
ten, die unmittelbar mit der iibertragenen Leistung in Verbindung gebracht 
werden kénnen. Multiplikation gleicher Zeilen bzw. Spalten muss 1, verschie- 


dener 0 ergeben. 


Aus der Spaltenmultiplikation von Q mit w* folgt: 


Spalte 2u4—1% x Spalte2 yu W*: 
Saal AD ee oe ee 

a ees ae Is Xu 
Spalte 2 u Wx Spalte 2 u— 1 U*: | 


* * aS 
2/2, Bia, + Be, My, 9 | 


(== Ties... ME) 


liefert je m reelle Gleichungen fiir je # Unbekannte. 


Spalte 2 4— 1% x Spalte 20 %*: | 


y2,, Bt +B, OF =0 
y=1 Xu Vd Xo 


Spalte 2 wu Wx Spalte2o0—1 wl: | 


XZ, BE +B, Lt =0 
pel Bu 


VOR UN LAN 
Bo & p Bo 


(6, Oeyh, Zac Mh) 


(uo. Vy 2,525, 9m 


(21) 


(22) 


(23) 


liefert je (m® — m)/2 verschiedene komplexe Gleichungen fiir je m2 — m Unbe- 


kannte. 


Spalte 24 — 1% x Spalte 20 —1%*: | 


a * ths 
2 ng Bi, + Bie My 


» 
be Bo 


ye (u = 0) 
7 10 (+o) | 


» (pw, o = 1, 2, ...., m) 


(24) 


| 
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Biicfert m*” komplexe Gleichungen fiir 2m? Unbekannte. Es ergeben sich ins- 
gesamt 2 m® + m verschiedene Gleichungen fiir 4 m? Unbekannte. Da bei m 
_ Stromungen 4 m? Koeffizienten mit 8 m2 Bestimmungsstiicken (Betrag, Phase) 
vorliegen, sind durch die Wirkleistungsbilanz die Halfte derselben dadurch fest- 
4 gelegt. Es zeigt sich, dass man keineswegs beliebige Bestimmungsstiicke vor- 
_ schreiben und die andere Halfte daraus ermitteln kann. Multipliziert man (20) 
_ bis (24) der Reihe nach mit 2, 2+, B, B*, L, L*, B, BF, L, Bs, so stellen 
_ die linken Seiten von (20) und (21) die deppae Wirkleistnte ie reiner 2- bzw. 
B-Modulation dar, wahrend der Realteil von (22) bis (24) die wechselseitige 
P 

é 


EN ee ek ee ae 


_ Wirkleistung der Endmodulationen ergeben. Nur wenn die 2- und %-Anfangs- 
_ modulationen in der gleichen Strémung wirken, ist die Wirkleistung der End- 


- modulationen nicht Null. Durch Zeilen-Zeilen-Multiplikation von % mit * 
_ ergeben sich weitere Gleichungen, die aber keine neuen Beziehungen fiir die 
_ Elemente liefern, da simtliche Méglichkeiten der Wirkung beliebiger Anfangs- 
- modulationen durch die obigen Gleichungen bereits erfasst sind. Es handelt 
_ sich dabei um Umformungen, die aus den Gleichungen (20) bis (24) durch Ver- 
_ tauschung der gesteuerten mit den steuernden Modulationen hervorgehen. Es 


Biolet fir ,a0= 1, 2,...,m 


: es Lig ig + Lug Lig = 9: (25) 
» Big Big + Big Big =, (26) 
m 1 (=o), 
ee a (27) 
i alate » 10 (u+#a) 


Durch Bildung der Determinante von (19) ergibt sich, dass | {| komplex 
ist und den Betrag 1 hat. Dies bedeutet, dass bei periodischen Vorgangen keine 
Anfangsmodulationen eingestellt werden kénnen, bei denen die Endmodula- 
tionen Null sind. 

Bei den bisherigen Ausfithrungen findet eine Zu- oder Abstrahlung elektro- 
magnetischer Wirkleistung in bezug auf die Gesamtstroémung nicht statt. Um 
das Verhalten der Stroémungen gegeniiber dusseren anregenden elektromagne- 
tischen Feldern und insbesondere die Vierpoleigenschaften zu gewinnen, werden 
zwei Klemmenpaare 1 und 2 eingefiihrt. An diesen Klemmenpaaren werden 
Hilfsstr6mungen m+ 1 (an 1) und m+ 2 (an 2) angebracht. Wenn die An- 
schliisse punktférmig sind, so ist die Laufzeit dieser Hilfsstro6mungen Null. Ist 
dies nicht der Fall, so sollen die Ladungstrager der Hilfsstrémungen die grésste 
Geschwindigkeit, und zwar Lichtgeschwindigkeit c, besitzen. Dabei ist die Lauf- 
zeit t. In beiden Fallen sind dadurch viele Koeffizienten der Hilfsstrémungen 


UNA JouANNES MULLER 5 ‘ ZAMP 


Null. Samtliche Koeffizienten, die B,,, bzw. By»+2 als steuernde Modulationen 
aufweisen, sind bis auf 


eoet, —jwt, 1 —— —JOTm 192 
iT m+ ee, t 
Brstg , Binsrg 


Null. Ausserdem sind die Koeffizienten, die &,,,, bzw. Lm 2 als gesteuerte 
Modulationen enthalten, bis auf 


== git = g71?tm+2 


+1 


m+1 g 
m+1 , m+2 
Sn. Ln +2 


Null. Im folgenden wird mit dem punktférmigen Anschluss (t = 0) gearbeitet. 
Um die Verbindung mit dem ausseren System zu erhalten, wird 


Bri = ee BVmig= Ug und Lyi = 5,7) » dens = sary 


gesetzt. Durch die Wahl der Richtung der Einstr6mung wird erreicht, dass bei 
positivem Realteil von Uf §* Wirkleistung von aussen zugefiihrt wird. 2 und Y* 
wird mit den Koeffizienten der beiden Hilfsstrémungen gerandert. Multiplika- 
tion gleicher Zeilen bzw. Spalten muss 1, verschiedener 0 ergeben. Die Be- 
ziehungen fiir die Koeffizienten der m Stro6mungen bleiben dabei ungedndert. 


Spalte 2y—1 2 2m+1 2m+2 2m+3 2m+4 
ar | Zeile 
€ | 
ae HM “gy = ics 0) | 2 fh = if 
HQ HR, Bag 0 Mx 0 2p 
1 he 
0 0 1 0 0 0 (2m 4- 1 (28) 
Ase 1a ae 1 = 1 Fie 0) 2m + 2 
0 0 0 0 1 0 |lam+3 
i ; 1) Ee Up, 10) —U,. 1 || 2m 4 
xy Sy AVE =u || 
ee * 
g 
aN Big 0 = Be. 0 — Brig 
: g 0 Q 0 P 
HY anu HX, Lig, 
x = (29) 
iain ing 2 + jo Pilg Re we gy 
O 0) 0 1 O 0 
Us, Us 0 Ua. il — Up, 
0 0 0) @) 0 1 
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Zeile 2m + 2% x Zeile 2m +3 U*: 


_ Der Ausdruck fiir die wechselseitigen epee Ry. + RE ist ein Mass fiir 
_ die in der Anordnung umgesetzte wechselseitige Wirkleistung. Ist Ry. + R* 
| Null, so bedeutet dies, dass die Wirkleistungen, die bei einer Anregung an 1 
bzw. 2 allein in der Anordnung auftreten, einfach superponiert werden kénnen 
bei einer gleichzeitigen Anregung an 1 und 2. Wird Index 2 gleich Index 1 
-gesetzt, so erhdlt man aus (30) den doppelten Eingangswirkwiderstand an 1: 


me aS UK. aad Ri se Re, a ee ig Us, oe Yh Ue (31) 


a on “ere 


_ oder auch aus Zeile 2 m + 2% x Zeile 2m + 1 W*. 
'  Zeile 24-1 x Zeile 2m+1%*: 


; HS, 
> Zeile 2 uM x Zeile2m+1U*: 
= oe Mie T ech Ut D3 (33) 


Die Koeffizienten der Endmodulationen gegentiber den Einstromungen am 
- Klemmenpaar sind durch die Koeffizienten der m Strémungen und den 2 m 
_induzierten Leerlaufspannungen am _Klemmenpaar vollstandig bestimmt. 


Spalte 2m + 3 UW x Spalte 2 m + 2 %*: 


2) 
r 


Uy. + Uw, = Ris a i ve, = aa a a _ Big ai Big iy, : (34) 


a 
3 Spalte 2¥— 1% x Spalte 2m + 2 U*: 
= 32g Bry ett Bie Vig, ; (35) 
| Spalte 2 y Wx Spalte 2m + 2 U*: 
Yay -3 Las, Bie, 5 Busy, ag (36) 


Die Formeln miissen so beschaffen sein, dass daraus erstens die unabhangige 
Serienschaltung bei entkoppelten Einzelstrémungen hervorgeht. Zweitens muss 
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» in den Formeln der Fall enthalten sein, dass die Str6mungen, wenn sie sich nur 
in ihrem Gleichstromanteil unterscheiden, das Verhalten einer einzigen Str6- 
mung zeigen. Die erste Forderung wird in den quadratischen Formen fir 
Rio + R* dadurch erfiillt, dass keine Glieder auftreten, deren beide Faktoren— 
verschiedenen Strémungen angeh6éren. Bei den Beziehungen (32) bis (36) sind 
bei entkoppelten Einzelstrémungen alle Koeffizienten fiir » + Null, so dass 
ein Ausdruck verbleibt, der nur die Koeffizienten einer Stré6mung enthalt. 
Im zweiten Fall ist die Wirkung einer 2-Modulation in jeder Teilstr6mung die 
gleiche, und die Wirkung einer %-Modulation der einzelnen Teilstromungen 
unterscheidet sich entsprechend dem Gleichstromanteil, abgesehen von Y,,, , 
Bug, und Yj... Hierbei gilt zum Beispiel: : 

4 4 
fury a es Big = [Be — ot") re | 
(37) 
firv=u: 2, = [2 — e197] “ a a | 
Bu 


Gleichung (30) geht unter Beachtung, dass 3’ 7,/: = 1, iiber in: 
v=1 


Re + Rh = D7 Us, UW, 2 + Uh, * Us, = Us, Us, + Hh, UE. 68) 
yp=1 “F - = ” 


— 
2B 4 ea 2 g 


Aus Gleichung (32) wird mit (37): 


m . . . 4 

i, vhs 7 i ter: 2 os 

Lig +: » Uy, mys [Lo are a nas “tr Uy. as [[£e ree od — + é ne 
Ct valet He 


u 


Tare gle’ (39) 
te Dy Ue ae 
p=1 we 
} 
= [Ui 29 + Ut, Sy] *. 
ie os Z 


In analoger Weise ergeben auch die Gleichungen (33) bis (36) das Verhalten 
wie bei einer einzigen Strémung. ; 

£- und BY-Modulation sind in bezug auf die Ergebnisse der Wirkleistungs- 
bilanz gleichwertig. Eine Vertauschung von 2 und Q in (20) bis (27), (30) bis 
(36) bei den Anfangsmodulationen allein, den Endmodulationen allein und den 
Anfangs- und Endmodulationen gleichzeitig fiihrt insgesamt wieder zu den 
gleichen Beziehungen. Der Ubergang von leer laufenden Klemmenpaaren zu 
kurzgeschlossenen Klemmenpaaren kann bei Einhaltung der Richtung von 3 
und WU durch Vertauschung von § mit Mf in den Formeln vorgenommen werden. 
Die wechselseitigen Impedanzen werden durch komplexe Steilheiten und die 
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war 


Wirkwiderstande durch Wirkleitwerte ersetzt. Aus (30) bis (33) wird: 


q 

4 ~ =f Ce a ae S as S* = S ~ Ok ~ N* 

j Siy, 7 By, = Sr ah 2 Sig, Sey is Sig, J2Q. , (40) 

>: 

+ Sty, - Sty, = 6,+ Sf, = >, ie is as ie taal (41) 
pal 7) yp v xy 

: 

Sn, i x, Sig, Lue, a Sig, Lag, : i) 

: m 

Bi = a, Sig, Bug + Sig, Bus, (43) 

2 on 

ei. : 

: Sind weitere Klemmenpaare vorhanden, so treten o und xin %, U, SG an Stelle 

_ von 1 und 2. Liegt Riickwirkungsfreiheit vom x-ten zum g-ten Klemmenpaar 

S vor, so ist R,,, bzw. G,,, Null. 

Bei nicht rein periodischen zeitlichen Vorgangen kénnen die aus der Lei- 

E stungsbilanz gewonnenen Beziehungen aus der Y-Matrix nicht unmittelbar 

_ benutzt werden. Wenn die Anfangsmodulationen nach einer rein reellen, zeitlich 


_ansteigenden Exponentialfunktion verlaufen, so ist die am Anfang der Str6mun- 
_ gen zugefiihrte Leistung stets grésser als die Leistung am Ende der Strémungen. 
_ Trotzdem k6nnen die fiir rein periodische Vorgange erhaltenen Beziehungen 
fiir nicht rein periodische Zeitfunktionen verwertet werden. Bei einer rein 
periodischen Zeitfunktion exp 7m? sind die Koeffizienten Funktion von 7 w. 
Bei nicht rein periodischen Zeitfunktionen mit komplexer Frequenz wird in den 
_ Formeln die reelle Frequenz durch die komplexe Frequenz ersetzt. 

Ladung und Masse der Elektronen treten direkt in den Koeffizientenbezie- 
hungen nicht auf. Die Beziehungen gelten deshalb nicht nur fiir Elektronen, 
sondern fiir beliebige Ladungstrager. 

Es ist keineswegs erforderlich, dass die einzelnen Stro6mungen an den Be- 
- grenzungen des Hohlraums beginnen und enden. Anfangs- und Endpunkt 
_ k6nnen auch innerhalb des Hohlraums liegen. Zur Aufrechterhaltung der Kon- 
 tinuitat des Wechselkonvektionsstromes ist dann an diesen Punkten eine Wech- 
- selladung entsprechend + 7m Q = % (— Endpunkt, + Anfangspunkt) wirksam, 
_ die auf die Giiltigkeit der Wirkleistungsbilanz keinen Einfluss hat. 

Der Weg einer Strémung (Index v) werde in zwei Teile geteilt. Die Lauf- 
- zeit langs des ersten Wegstiicks sei t,,. Zwischen den Koeffizienten, die durch 
~ 2-Modulation auf dem gesamten Weg und den beiden Teilstiicken vorgenommen 
_ wird, besteht dann folgende Beziehung: 


Yo, = Do, + De, 7°: (44a) 


a dies ehh. winneaiiie 


se ee 


9) ist eine gesteuerte Modulation, die der v-ten Strémung selbst nicht angehort. 
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Ye, gibt zum Beispiel die Wirkung an, wenn eine 2-Modulation nur langs des” 
zweiten Teilstiickes besteht. ; 

Die Ladungstrager sollen den v-ten Weg einmal in der einen Richtung, 
@-Modulalion und in umgekehrter Richtung @-Modulation durchlaufen. Sei 1, 
die maximale Laufzeit und 1,(s) die ortsabhangige Laufzeit, gerechnet vom 
Anfang der &-Modulation, so ist die erzwungene Verteilung des Wechselkon- 
vektionsstromes durch die &- bzw. &-Modulation, & exp —7 @t,(s) bzw. 
a exp —j w(t, — 1,(s)]. Angenommen, die Str6mungen in dem verlustlosen 
elektromagnetischen Hohlraum haben alle verschwindende Raumladung, so 
sind die magnetischen Feldstaérken und Stréme in den Begrenzungen des Hohl- 
raumes, herriihrend von der Stromstarke eines infinitesimalen Stromelements, 
mit der Stromstarke desselben in Phase. Sei h(s) der Proportionalitatsfaktor 
zwischen Stromstarke eines Stromelements und der interessierenden Strom- 
starke an einer beliebigen Stelle des Hohlraums, so bestimmen sich die Koeffi- 
zienten Jy, und Jg zu: ‘ 

Sm Sm 
€ = | Hs) CRAY ss Nee tes a h(s) e~7? [»—5)] ds . 
Gre: 6 
Hieraus folgt fiir die Beziehung zwischen den Koeffizienten: 
Se = ae 2: 1°™. (44b) 


Bei den elektrischen Feldern und Spannungen an beliebigen Stellen des Hohl-- 
lt : : 
raums gilt eis te piety (44) 


Bei der Ableitung tritt an die Stelle der erzwungenen Verteilung des Wechsel- 
konvektionsstromes die der Wechselladung. Derartige Umformungen werden 
im Abschnitt IT bei anfangsgesteuerten Str6mungen bendtigt. 

Fiir den einfachsten Fall einer Strémung mit einer Tragergeschwindigkeit 
soll ermittelt werden, welche Bestimmungsstiicke der Koeffizienten frei gewahlt 
werden kénnen. Es wird laufend durchnumeriert. Ungerader Index bedeutet 
%-Modulation, gerader Index B-Modulation. Es wird gewahlt: 


—_ a 19u- — —_ 1Pi2- 
dig Qi = A775 Lig. = Ay, = Ay, 0! P| 
= p=. 1 Par pee Se t Poo 
Bi, =ioy Go) Big. = Ugg = Apne *™ . 
Zeilenmultiplikation von %{ mit Y* ergibt: 
| 
Mi, Aye | a A a3 
QT — , = ’ 
M91 Ag9 ais ait 


M41 O35 + Ay2 OF, = 1; a4, 0% + yy ay, = 0; ay, 0% + aah, =0. (45) 


Ri Eh en =" 4 b fet ene ht ss a 


ae J ee | 
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Durch Umformung aus (45) gewinnt man: 


B69 ET eee 
ay ass ats ay, Me) 
_ und hieraus die Phasenbeziehungen: 
at 
Pu Var ete Meg Pats. (47) 


-Dabei sind die Phasen um z unbestimmt. Das Vorzeichen wird deshalb zu den 


a-Werten genommen, die nicht als die stets positiven Betrage anzusehen sind. 
Aus der ersten Gleichung von (45) folgt: 


Ay Ag9 + Aya Ay, = 1. (48) 


Damit k6nnen ein Phasenwert und drei a-Werte frei gewahlt werden. Drei Phasen- 
werte und ein a-Wert sind dann festgelegt. Die Beziehung a,,/a*, = —a»,/a%, 


die nicht unmittelbar aus (45) folgt, kann direkt durch Multiplikation der 


Spalte 1 von Y& mit Spalte 2 von Y* erhalten werden. 

Die geringste Anzahl der Koeffizienten bei m Str6mungen betragt 2 m, und 
zwar muss mindestens ein Koeffizient in jeder Zeile und Spalte von Y stehen, 
damit die Determinante nicht Null werden kann. Die Verteilung der Koeffi- 
zienten ist dabei nicht willkiirlich, sondern muss so erfolgen, dass bei der Bil- 
dung der a-Matrix das gleiche Bild fiir diesen Spezialfall der kleinsten Koeffi- 
zientenzahl entsteht, damit bei der Produktbildung von. YY %* eine Einheits- 
matrix zustande kommt. Sei bei m= 2 nur Qy3, Qo4, G32 uNd ay, nicht Null, so 
ergeben sich die Koeffizientenbedingungen zu: 


330% =1; agai=1. (49) 


Auch hierbei kénnen die Koeffizienten nicht willktirlich gewahlt werden. 


Differentialbeziehungen 


Wenn ein Geschwindigkeitsspektrum vorliegt, geht die Gleichstromstarke 
in jeder Strémung gegen Null. Damit kann die Gleichung (8) angewendet 
werden, deren Glied mit Az erweitert wird zu: 

= Ai ad A (50) 
vy OR gE Xe 


v 


Bu 


Zwischen der Wechselvoltgeschwindigkeit und dem Wechselkonvektionsstrom 
an einer beliebigen Stelle der y-ten Strémung, herrithrend von einer Modulation 
%X, besteht damit dieser Zusammenhang. Bei X kann es sich um die Anfangs- 
modulationen der anderen Strémungen oder um die Einstr6mung bzw. erzwun- 
gene Spannung an den Klemmenpaaren handeln, das heisst um solche pri- 
maren Modulationen, die eine Beeinflussung der »-ten Strémung nur indirekt 


218 JOHANNES MULLER : ZAMP 


auf dem Umweg iiber das elektrische Wechselfeld bewirken. U, (Abkiirzung 
fiir U,,) sei die statische Anfangsvoltgeschwindigkeit der Strémung y. Eine 


ahnliche Beziehung kann mittels (4) erhalten werden. Zu diesem Zweck wird _ 


| 
: 


die durch Feldmodulation der Gesamtstr6mung mit Geschwindigkeitsspektrum _ 
entstehende Konvektionswechselstromdichte /(s, x, y) und der durch - und 


%-Modulation in der »-ten Strémung erzwungene Konvektionswechselstrom 
&,(s) in die Gleichungen fiir die elektromagnetischen Felder eingefihrt: 


rot E=—JwomwH, rot H=jwoen€E+ Us, x, y) + Vis). (51) 
Dabei ist: 
Yeas : jw [e(s, x y)—7(8, x, 9)] 
* ° > dei lus, ¥ ¥)—7(5, % ¥ 
is, 2,9)= Ge | ge aU, | Els, x, y) ds ~—_G, , (52) 
Veryan Me 
defor) ber 
(Ss) ae Ai, B, (0) 7a &, (0) é : a (53) 


dU, 
[di/dU,] dU, ist die statische Stromdichte im Voltgeschwindigkeitsintervall U, 
und U, + dU, eine Funktion der Koordinaten x, y auf der emittierenden Ober- 
flache. Der Faktor dFj/dF tragt einem nichtparallelen magnetischen Fiihrungs- 
feld Rechnung. Aus der Gleichung (53) ist ersichtlich, dass die erzwungene Ver- 
teilung des Konvektionswechselstromes, herriithrend von einer 3-Modulation, in 
der y-ten Str6mung aus der erzwungenen Verteilung des Konvektionswechsel- 
stromes durch £-Modulation in der »-ten Str6mung mittels Differentiation nach 
U, und Multiplikation mit Az, entsteht. Zwischen den gesteuerten Modulationen 
Yy, und Ye, ausgelést durch B,- und L,-Modulation, existiert die Beziehung : 
Aa 
Ve, = 4, aa (54) 
Bei Y) handelt es sich um Modulationen, die durch reine Feldmodulationen 
zustande kommen, Wechselvoltgeschwindigkeiten, Wechselkonvektionsstréme 
der anderen Str6mungen, Leerlaufspannungen und Kurzschlu8stréme an den 
Klemmen sowie das elektrische und magnetische Wechselfeld. Die Bildung von 
(54) ist nur dann streng zuladssig, wenn U, in den durch die Konvektionsstréme 
erweiterten elektromagnetischen Gleichungen, abgesehen von den erzwungenen 
Verteilungen, nicht vorkommt. U, ist jedoch in /(s, x, y) enthalten, aber nur 
mit einer Str6mung mit infinitesimalem Gleichstromanteil verkntipft, so dass 
eine Riickwirkung zu vernachlassigen ist. 
Die Koeffizienten, die sich auf die gleiche Stromung beziehen, setzen sich, 
abgesehen von dem durch reine Feldmodulation entstehenden Koeffizienten 
Bro» aus zwei Teilen zusammen, die sich um eine Grdssenordnung unterscheiden. 


B, . (55) 


5) , 
PR ds 


a 
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_ Die ersten Glieder entstehen durch direkte Modulation. Aus (7) und (6) folgt: 


Deh ea aa aig aman -jor 
"S, v aU, ? a, SER, 7S, = Mi (56) 


Die zweiten Glieder entstehen auf dem Umweg iiber das elektrische Wechsel- 
feld und kénnen aber nicht aus dem Grenziibergang mu = v der Koeffizienten 


Logg) Bro : Bro gewonnen werden. Die letzteren lassen sich mittels der Diffe- 
rentialbeziehungen durch &%, a, in folgender Weise darstellen: 
AB, on a B,o 
Ley, = Ag, aU, = Mi, Ai, aU, qu, ; (57) 
AB, 
Gro a= Ae, aay : (58) 
AB, 
De i, aU, (59) 


Demgegeniiber muss fiir w = v auf die Bestimmungsgleichungen (2) und (6) 
zurtickgegangen werden. 


Ss A é SAG A 2 
. defol-H@) — f FEgl8) def ole)-H18)] 
Pe, i Eq (5) ds i. ange | 7 s (60) 
; de-io[rls) 718) 
v . e—1%[t,(s)—7,(5 
3, — 44, Je E,,(s) ds“ (61) 


0 
i j w[r,(s) —t4(8)] eae) j wt, (s) —7,(3)] 
pee —j o[t,(s)—t,($)] _ A, a + eee —7 w[T,(s)—7,(8 
B,, — | Ey,(s) dse Ai, | —zp— ase . (62) 
0 


LU, 
0 


Durch Addition von (61) und (62) folgt: 


0 . 5 d j j t,(s) —T,,(& 
Bg + Big = Ate ar | Co,l8) ds oF 9 (63) 
0 
dB,o 
= Mi, —. (64) 


Durch Vergleich von (64) mit (58) und (59) fiir 4 = v ist der Unterschied bei 
gleichem Index zu erkennen. Sind noch Einzelstromungen mit diskreter Ge- 
schwindigkeit und endlicher Gleichstromstarke vorhanden, so andern sich zwar 
wertmassig die Koeffizienten, die Differentialbeziehungen bleiben erhalten, 
kénnen aber nur auf die Einzelstrémungen selbst angewendet werden, wenn 


ihre Gleichstromstarke gegen Null geht. 
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Integralbeziehungen 


Unter Benutzung der Differentialbeziehungen werden nunmehr die durch | 
Matrizenmultiplikation gewonnenen Summenausdriicke umgeformt. Bei — 
Summe (22) werde die Umformung vorgenommen, die sich bei den anderen — 
Summen (20) und (21), (23) bis (27) in ahnlicher Weise wiederholt. Aus der 
Summe (22) werden vier Glieder fiir » = ~ und » =o herausgenommen, da — 
zwei Glieder ihren Wert sprunghaft andern. Es folgt: 


0=)'2, BY +B, Lt Wtpo 
ey fu & Ln Se ( M (65) 
tig DBE BOE ty Bie Om ee | 


Unter Beachtung von (55), (56) und (58) sind die beiden mittleren Einzelglieder 

von erster Ordnung in der Gleichstromstarke klein gegeniiber den beiden ande- 

ren Einzelgliedern und werden vernachlassigt. Bei Verwendung von (58) in der 

Summe ergibt sich (67). Wenn die Unterteilung des Spektrums in diskrete 

Einzelstr6mungen geniigend fein ist, kann mit y= «4 und »= o summiert werden. 

Der wahre Koeffizient &,,. hat zwar seinen Wert gegentiber 2,, sprunghaft 
g ae 


lt 
gedndert. Der Wert von 2, bleibt jedoch beim Ubergang zu » = u von gleicher 
7 


Ordnung, und wegen der in der Grenze infinitesimalen Gleichstromstarke der 
-ten und o-ten Str6mung ist der Einfluss der Erweiterung zu vernachlassigen. 
Es wird davon erst beim Ubergang zur Integraldarstellung Gebrauch gemacht. 
Die Ergebnisse der Umformungen von (22) bis (27) werden der Reihe nach — 
zusammengestellt. Es werden sémtliche Beziehungen (55) bis (64) bendtigt. 


-A+o, M=o, 
A,c=9, A,,=09, 
Aus (22) bis (24) ee as 
Ai, Ai, aU, d ‘Tia ‘ C= Oe 
Ai aA, o 0 Ai Au, 
a TO , a aU, =a) 
(66) 
U+o, hae; 
O= Ai, Ai, ee 
= At Aid, dU, dU 0 =D" 
Aus (25) bis ( 
ve B,C? 0= 5) 
Ly lh 
giedpereas 1s 
— a aU, , 0) = At, - ° | 
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Samtliche Gleichungen werden erfiillt, wenn vier Ausdriicke Null sind: 


ee ae ty one Se 
0=A,, iS) At, dU, [Bo, Br, oF Big ¢ 1OTm 4 Bee, ef O%or (isk sO) , 


v=1 
(67) 
“ : d j WT -—jJ wt 
[oe = 2) 4, aU, [Bre, Bee, ge a Bro, en Oty (68) 
Reea SA pe mel 4 2 Re Aiz [ 1 Gq (s) ds 4 ei?™u®), (69) 
oz gees "By x, uf dU, *« aU, 
0 
= d 
= = * 
- j of y(Sm) ~7 y(3)] (7) 
; ad F de~1?|7y'Sm) — 7 y(8 
2 j ot ,,(sm) 
+2Re Az | ag [ 2,9) Orn | 4% aU, | 
0 


(67) und (68) sind von nullter Ordnung, (69) und (70) von zweiter Ordnung in 
der Gleichstromstarke. 
Es ist sinnvoll, Str6mungsimpedanzen einzufihren. In der yu-ten Strémung 


- werde eine £-Modulation vorgenommen. Dann ist am Ende dieser Stroémung 
_ der Konvektionswechselstrom &,, exp —]7 w t,. Die wechselseitige Strémungs- 
impedanz r, ,, sei das negative Verhaltnis der Wechselvoltgeschwindigkeit am 


Ende der o-ten Strémung zum Wechselkonvektionsstrom am Ende der y-ten 


Strémung. et es ee Pte (71) 


o, fh oa, 
Die Stromungsimpedanz der y-ten Strémung selbst betragt: 


t, p= =-B,. e°. (72) 
: 8, | 
Durch das negative Vorzeichen wird erreicht, dass die Strémungsimpedanzen 
direkt in die Impedanzen iibergehen, mit denen die Anordnung fiir aussere 
Anregungen wirkt. Seien w und o laufzeitlose Stré6mungen an Klemmenpaaren 


o und x, so wird: 
Tyo Ro, x oe ee Roo: (73) 


t,,¢ und r, ,, in (67) und (68) eingefiihrt, ergibt: 


yap 2. 27° B ead 
a AT 
“ Xo 


(vy + wu, 0) 
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Fiir verschwindende Raumladung sind die Summen, da sie den gesamten 
Gleichstrom der Strémungen enthalten, Null zu setzen. Fiir zwei bzw. eine 
Stromung gilt dann: 

toga ese (75) 
beziehungsweise . 
+r% = 0. (76) 


a} Vad 


Angewendet auf zwei bzw. ein Klemmenpaar: 


Re RE =O, (77) 
Ror Re = Oe : (78) 


‘Das sind die bekannten allgemeinen Impedanzbedingungen fiir eine Anordnung 
ohne innere Verluste. (69), (70) zusammen mit der entsprechend (79) umge- 
formten Beziehung (76) 

Sm 


+ rk, =2ReB,, ef =2 Re | €y,(9) ds ei?) —0,. (79) 


Thu 
7 


0 

stellen nichts anderes als die notwendigen Bedingungen fiir die Verteilung des 
elektrischen Wechselfeldes, ausgelost durch eine 2 ,-Modulation, langs der u-ten 
Strémung bei verschwindender Raumladung unter Benutzung der Differential- 
beziehungen dar, damit die Wirkleistungsbilanz erfiillt ist. 

Fiir einen einzigen Stromfaden mit einem Geschwindigkeitsspektrum, des- 
sen obere und untere Grenzen der statischen Anfangsvoltgeschwindigkeiten U, 
und U,, sind, ergibt sich aus (74) durch Anwendung einer partiellen Integration: 


dt . ; 
Bp es ht 1 OT 7@T,)%* 
di, . , 
=| eo [MR —_ g@ Ty] * 
DB, 
et dU B JOT, B jor, ) x 
ad “adU2 v va e VQ e | ? 
¥ v <o Cu 
Uy 
di, di 
% | ee * ey Se * 
Hea "i ay fab Bus, Br, eer Bus, Bre, | 
i 81 
* ai, Fe are (81) 
ey dU2 fai) HQ, Bi, > 
Uu 


Bei einem Stromfaden bedeutet der Index » nur eine Abhangigkeit von den An- 
fangsvoltgeschwindigkeiten. Bei mehreren Stromfaden bedeutet Index » auch 
eine Abhangigkeit vom Ort des Emissionspunktes, so dass dann eine weitere 


Bi a i ects 


aN ee ee OP RT TL eee Oa! ee ek 
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Integration oder Summation in (80) und (81) erforderlich ist. Unter einem 
Stromfaden wird eine Anordnung verstanden, bei der die Felder in einem 
Querschnitt senkrecht zur Hauptbewegungsrichtung als konstant angesehen 
werden kénnen. Ein Stromfaden kann deshalb zu einer ebenen Strémung aus- 
arten. Bei endlicher Gesamtgleichstromdichte im Stromfaden richtet sich der 


- Gesamtgleichstrom des Stromfadens nach der zuldssigen Querschnittsflache. 


Die weiteren Untersuchungen sollen sich auf das konstante und exponen- 
tielle Spektrum erstrecken. Fiir das konstante Spektrum gilt 

dt, dt 

dU, au, 


= Const. 


Gleichzeitig hat diese Verteilung auch die Eigenschaft, dass der Quotient aus 
der statischen Ladungsdichte dg langs des Stromfadens und dem zugehérigen 
Geschwindigkeitsintervall du ebenfalls konstant ist: 


dq e di 


Dihigee m.- AU= 


_ Ein Stromfaden mit beliebigem Spektrum werde einer reinen Feldmodulation 


unterworfen. Die »-te Strémung mit dem Gleichstromanteil Az, und den 
Anfangsvoltgeschwindigkeiten zwischen U, und U,+ AU, erfahrt an einem 


_ beliebigen Punkte der Strémung eine Wechselvoltgeschwindigkeit 


ne. 


Bi(s} und B(s) + AU 


Der Gleichstrom Ai, verteilt sich jetzt auf ein Spannungsintervall 


Av, [1+ Se]. 


_ Andrerseits kommt aber nach (50) ein Wechselkonvektionsstrom zum Gleich- 
stromanteil hinzu, so dass fiir jeden Punkt s der Strémung der Quotient 


[di,/dU,], unverandert bleibt: 


di, adie 
tap — AU, [i + 4B,(s)/dU,) ~~ aU, ° 


Ai, [1 + dB,(s)/dU,) _ di (82) 


Dies gilt fiir jedes beliebige Spektrum. Speziell beim konstanten Spektrum 
tritt damit durch die Feldmodulation keine Verainderung der spektralen Ver- 
teilung ein. Wird (50) integriert, so ist der gesamte Wechselkonvektionsstrom 
@() am Ende des Stromfadens mit einem konstanten Spektrum nur durch die 
Wechselvoltgeschwindigkeiten $8) und ¥) an der oberen und unteren Ge- 
schwindigkeitsgrenze bestimmt. 


Ga) a0. [p00 — Bw]. (83) 
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Es wird eine wirksame Voltgeschwindigkeit als mittlere Voltgeschwindigkeit 
der Endmodulationen eingefiihrt : 


1 


ms 0 u) } 
Bi) = = [BO + Bw]. (84)5 


Analog sei fiir die Anfangsmodulationen, die durch eine feldmodulierte Vor- 
strecke entstanden sind, gesetzt: 


di 1 
= 7 [8 - BJ, B= 7[B— BJ. (85) 
Durch Aufldsen nach 8°), VB, B, und B, ergibt sich: 
0 dU, c 6) dU, ( ). 
eel rar ar Case Me, Re oY ode (86) 
dU, dU, 
B, = B+ aq B, = Bi — 4a, 1; (87) 


Mit den Wechselvoltgeschwindigkeiten an den Grenzen des Spektrums sind 
auch Wechselkonvektionsstr6me verbunden: 


1 a 1 a 
ig) : tt) Shee Dee Pe ee Le 
QO = > Qi + ris Bo; gw — — ga ae Bl): (88) 
il a 1 at 
= 7 it ag Bi; Lu= 2 g aU, B,. (89) 


Es soll nun untersucht werden, ob mit diesen Ersatzmodulationen der 
Stromfaden mit konstantem Spektrum in seinem Verhalten auf eine Einzel- 
stro6mung mit diskreter Geschwindigkeit zuriickgefithrt werden kann. Zu diesem 
Zweck wird die Wirkleistungsbilanz mit diesen Ersatzmodulationen gepriift. 
Entsprechend dem Ausdruck fiir die Wirkleistung N = Re % &* ist die infini- 
tesimale Wirkleistung am Ende der Str6mungen durch eine reine Feldmodula- 
tion (Index ©) in der Hauptstrecke: 


1 di 1 dt 
aN = cea a [Wg dBVe + VE dBe| = eR d| Beg BE] (90) 
und fiir das gesamte Spektrum: 
N' 1 di (0) (0) % (w) (0) x 


Analog ist die zugefiihrte Wirkleistung am Anfang der Strémungen in der 


Hauptstrecke: 1 


N,= ay a0 [B, Vy bs B,, Bi ; So. 


Die zusammengesetzten Wechselvoltgeschwindigkeiten an der oberen und un- 
teren Grenze des Spektrums am Ende der Strémungen betragen unter Beriick- 


- 


ee eS Le ol eee oe eee CWerente bik band 2 — 
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“sichtigung der Laufzeiten t, und t,, 
BP) = B+ Before; yw) — Be) + B efor, (93) 
Bam Ende der erases wird dann die Wirkleistung: 


IN a (c ) glo) 4 QW B*] 
1 di ‘ . “ . 1 ; 
=>s pee (04) 
at. iG (°) ei to _ Be Be ef tu] 


abgegeben. Wegen der zusdtzlichen Str6me an den Geschwindigkeitsgrenzen 
durch. die Anfangsmodulationen muss das elektrische Wechselfeld noch eine 
zusatzliche Wirkleistung N” aufbringen: 


sm 
di j WT o(S j wTy(3 
=Re | G(s) ds eat rl = Bp ahem 
0 


Br civro | &(s) ds e 1 [to-tol8)] (95) 


0 


Sm 

— Rx crore | Es) ds seid 

% : 
0 


Die beiden Integrale in (95) auf der rechten Seite stellen nichts anderes 


als Bi?) bzw. BY dar. Damit ist N” gleich dem letzten Glied von (94), und 


wegen (91), (92), (94) und (95) ist N‘) = N’+ N,+ N". Die durch das Wech- 
selfeld der Hauptstrecke zugefiihrte Wirkleistung N’+ N” ist gleich der Dif- 
ferenz der von den Ladungstragern abgegebenen und zugefiihrten Wirkleistung. 
Eine Wiederholung des Beweises wurde vorgenommen, um die speziellen Aus- 
driicke beim konstanten Spektrum zu gewinnen. Nunmehr werden die Ersatz- 
modulationen aus (87) und (86) in (92) und (94) eingefiihrt: 


1 
N, = > [B. SF + 2, BF]; NO = > [BP ah + BP a]. — (96) 
Das sind die Wirkleistungen fiir eine aap Aus (93) erhalt man: 
09 = OP + 8p [eter + cio] + B, re [etere— een], (07 


; ie af au —Jwr ~jor 
BY = BY + B, > [eter + e }+ 2, 2- Cs A a (98) 


Lasst man bei konstantem Gesamtstrom die Breite des Spektrums klein werden 


A 
mat 5: Tatts: wo Ar<l1, 


_ZAMP V/15 
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so entsteht: 
o de~i® 
QP = + Qe fe Bt 


a * BO — BY +B, cie™. (99) 


Das sind entsprechend (4) und (6) die Zusammenhange fiir eine Strémung. Ist 


w At > 1, so kann (97) und (98) nicht durch Wahl einer Ersatzlaufzeit und : 
einer geeigneten Abhangigkeit der Laufzeit von der Anfangsvoltgeschwindig- _ 


keit unabhangig von den Anfangsmodulationen auf eine Strémung ohne Spek- 
trum zuriickgefiihrt werden. Das gleiche gilt fiir die wechselseitige Leistung 
zwischen Feld- und Anfangsmodulationen N”, die ebenfalls die Laufzeiten ent- 
halt. PO werden die durch das Wechselfeld entstehenden Modulationen 
2@ und BY) lineare Funktionen von den Konvektionswechselstrémen 2%, und 
Sy Sein, thd Q,, entstehen durch Verschiebung und Anderung der Breite des 
Spektrums am Anfang der Str6mungen und enthalten nach (89) 2, und Q,. 
Insgesamt lassen sich also die Ersatzendmodulationen durch die Ersatzanfangs- 
modulationen nach (97) und (98) wie folgt darstellen: 


Ce he hg Br SP Be SB he (109) 


In dieser Beziehung besteht kein Unterschied zwischen einer Str6mung mit 
konstantem Spektrum und einer Stro6mung mit diskreter Geschwindigkeit. 
Diese linearen Zusammenhange zwischen den Anfangs- und Endmodulationen 
bilden mit der Invarianz der am Anfang und Ende der Str6mung auftretenden 
Wirkleistungen, deren Wert durch die Ersatzmodulationen wie bei einer Ein- 
zelstromung wiedergegeben wird (96), die Voraussetzungen dafiir, dass die Er- 
gebnisse der Matrizendarstellung (20) bis (27), (30) bis (36), (40) bis (43), (45) 
bis (48) verwendet werden kénnen. Damit kénnen diese Koeffizientenbezie- 
hungen auf Strémungen mit konstanten Spektren bzw. eine Kombination mit 
Str6mungen diskreter Geschwindigkeiten erweitert werden. Der Unterschied 
zwischen dem Verhalten der beiden Str6mungsarten besteht in der verschie- 
denen Ortsabhangigkeit, die in den Koeffizienten selbst nicht direkt erscheint. 

Es werden die Umformungen angegeben, die fiir eine gesteuerte Modulation 
Y beziiglich ihrer Abhangigkeit von den an den Grenzen des Spektrums tat- 
sdchlich auftretenden steuernden Wechselkonvektionsstrémen 2, und @,, bzw. 
ihren Ersatzmodulationen bestehen. 


D= Vo, 2+ Do, Su= La > [o, + Yo, ]+ Sc ro dem Baul) Oe 


Seth 


&, und &,, sind aus (89) durch 2, und B, ersetzt worden. Die auf die Ersatzmo- 
dulationen bezogenen Koeffizienten lauten: 


Vo,= + (9s, + Deals Vg. = i [Do, — Va, (102) 
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Boder 


1 a iy EHO 
Yo, ee es Wi Yo, Yo, = Ve, ae SE 


a 


(103) 


_Geht die Breite des Spektrums bei unverandertem Gesamtstrom gegen Null, so 
wird Yo, = Yeo, = Ye,, das heisst, die Wirkung einer 2-Modulation unmittelbar 
-rechts und links des auf eine Geschwindigkeit zusammengezogenen Spektrums 
ist die gleiche. Unter Benutzung von (86) und (103) in (80) und (81) (die Inte- 
_grale entfallen wegen d*7,/dU2 = 0) ergibt sich: 


Dorit Fx oo [Brg eietuls 4 By xed %e [29 eiete|* (104) 
acto 


we Sy ACRE RE ne abana de : 


Mu a 

Pps BT iD, gt BU | (105) 

Auf Klemmenpaare 1 und 2 angewendet oa =, = 0,2, = =n Lo = — Ses 
B,= UU, B, = u,) folgt: 

‘ Ris aig Re 3 L,. Bi Si Bi or. a Uy. UW, F U,.. Be (106) 


: in Ubereinstimmung mit den Gleichungen (34) und (30) fiir eine Strémung. 
_ Umgekehrt stellen (104) und (105) Beziehungen fiir die Str6mungsimpedanzen 
_ bei einer Strémung mit diskreter Geschwindigkeit dar. 

i Fir das exponentielle Spektrum nimmt di,/dU, und seine erste Ableitung 
die folgenden Werte an: 


di, et — _[ a, ou | & svete 
wae he = lar. = ; Abaar Sram | 
(107) 
oo, ees ae 
lee oh t es ae 


T Temperatur des Emissionsgebietes des Stromfadens mit dem Gesamtgleich- 
strom 7. Die wechselseitigen Stromungsimpedanzen aus (80) und (81), der Stré- 
-mungswirkwiderstand 7,,,, [wird aus (108) und (109) erhalten fiir o = w], und 
die Beziehungen (69) und (70) ergeben sich damit, wenn in allen Gleichungen 
mit 2k T multipliziert und eine Umstellung vorgenommen wird, zu: 


Upg= 


2e dt, B,, Bro, = 247 B Bia, +2kT[x,, + 1% ,] ei [To] , (108) 


Ug. 
Uy=0 


2 edi, B,. Bl, =2¢i By Bi +2T [t,o +t ,], (109) 
0 


[2edi, Bs Bi =2ei Bio Bi AAT, « (110) 
EB 
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[2edi, Bio Bh, =2e7 Bue Bio +4kT Hy» (111) 
0 


[edi B_, Bi, = 267 Big By —4hT Re Ai; 


d de! ®!*) 
iri Es, (8) 4s — a 


Lu 
0 


[ edi, 2,. Mt, = 2etL,, Ot, —4 bT Re Ad 


ins (113) 
t ; —7 @[t (Sm) —T,(3)] 
ee j @Ty(Sm) del sate 
Sf caus (eg enk ane ee aU, 


Die ersten vier Gleichungen sind von nullter Ordnung, die beiden letzten sind 
in allen ihren Gliedern von zweiter Ordnung in Ai, da die Koeffizienten Beg 
und 2, von erster Ordnung in Ai, sind [(59), (58)}. 2 e di, ist das infinitesimale 
Schwankungsquadrat!) pro Hertz Bandbreite 2, 2* der v-ten Str6mung mit 
der infinitesimalen Gleichstromstarke di, an der Emissionsstelle. 2 e7 ist das 
Schwankungsquadrat des Emissionsstromes des gesamten Stromfadens. (108) 
bis (113) ergeben somit zwanglos Beziehungen, die auf Rauschvorgange in 
linear sich verhaltenden Anordnungen bezogen werden kénnen. Die linke Seite 
von (111) stellt das Schwankungsquadrat der Voltgeschwindigkeit der La- 
dungstrager am Ende der w-ten Strémung dar. Es setzt sich zusammen aus 
einem Schwankungsquadrat der Voltgeschwindigkeit am Ende der yu-ten Stré- 
mung, das entsteht, wenn der gesamte Schwankungsstrom in der Strémung 
mit der niedrigsten Geschwindigkeit (U, = 0) wirkt und einem Schwankungs- 
quadrat einer Spannung, das man erhalt, wenn man dem Str6mungswirkwider- 
stand 7, ,, entsprechend der Rauschformel von Nyqgurst?) eine Rauschtempe- 
ratur 7 zuordnet. Bezieht sich w auf ein Klemmenpaar 1 (%,, 7,,,,, gehen iiber 
in U,, R,,), so entsteht aus (111): 


co 


[2edi, My, Ut, = Zed Uh, Uh, +42T Ry. (114) 


0 


Beim Ubergang vom Widerstandersatzbild zum Leitwertersatzbild tritt fiir 


1) Alle vorkommenden Schwankungs- oder Rauschquadrate oder -produkte beziehen sich 
stets auf 1 Hertz Bandbreite. 
*) H. Nyoutsr, Phys. Rev. 32, 110 (1928). 


en 


- 
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YU, und R,,, 3, und S,,: 


(ee edi, Sig Sy = 207 Sh +4 kT Sy. (115) 
0 


_ Die Zerlegung in zwei Rauschanteile ist von FREEMAN?) fiir ebene Stromungen 


bei verschwindender Raumladung berechnet worden. Sie kann hier auf den 


_ Fall mit Raumladung im elektromagnetischen Feld und fiir die Randwerte der 


_ Strémungen selbst erweitert werden. Im Abschnitt II iiber anfangsgesteuerte 


Strémungen wird gezeigt, welche Beziehungen bei Anwesenheit eines Potential- 


_ minimums bestehen. Die Schwankungsstréme in verschiedenen Strémungen 
_ einer emittierenden Oberflache stehen in keinem phasenmassigen Zusammen- 
: hang”). Dagegen besteht eine phasenmidssige Zuordnung bei den durch die 
_ Emissionsschwankungen ausgelésten Schwankungen. Fiir die Schwankungen 
_ der Voltgeschwindigkeiten an den Enden der Strémungen regelt dies (109) und 
_ fiir die Schwankungen der KurzschluBstréme an zwei Klemmenpaaren 1 und 2 
& (116) (in (109) wird B,., Bq, Ty,0 Vo, « ersetzt durch 3,; 32, Siz, Sai]: 


f 26d, Sig Shy = 205 Sig Th + 2kT (Set Sh]. (116) 
0 


‘Die Integrale in (109) und (116) sind bestimmend fiir den wechselseitigen 
Einfluss bei der Bildung wirksamer Schwankungsquadrate. Sind die Integrale 


- in (109) und (116) Null, so kann gerechnet werden, als ob die Schwankungs- 


quadrate an mw oder 1 bzw. o oder 2 unabhangig voneinander sind. 
Die linke Seite von (113) stellt das induzierte Schwankungsquadrat des 
Konvektionsstromes am Ende der y-ten Stro6mung dar, das proportional Az? 


ist. Da das Schwankungsstromquadrat der w-ten Strémung proportional der 
 Gleichstromstarke Az, ist, so tritt das induzierte gegentiber dem direkten 


Schwankungsstromquadrat in jeder Strémung bei verschwindendem 41, zu- 
riick. An Klemmenpaaren liefert (113) und (112) keinen Beitrag. Die Glieder 
mit T werden Null wegen t,, = 0, und jeder einzelne Koeffizient wird in den 
anderen Gliedern Null. 

Aus (108) und (109) sowie (110) und (111) erhaélt man die Vertauschungs- 
relationen durch Subtraktion, wobei die Koeffizienten durch die Str6mungs- 
impedanzen ersetzt werden. 


fore) for) 
. . . * . 
; [2edi,t, 12, —2eity tha = 2edi,t, ot, —2ei toy, (117) 
0 


0 
1) J. J. Freeman, J. appl. Phys. 23, 1223 (1952). 
2) W. Scuorrxy, Ann. Phys. 57, 541 (1918). 
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| 20d, Tha 2 ete ten [208 if, — 2e%t, i Ss (118) 


und fiir uw = u: 


[red t,t z.~ [20a ty thy. (119) 


Die linken Seiten gehen aus den rechten Seiten durch Vertauschung von 
y mit uw, » mit o, w mit w und uw mit o sowie Bildung des konjugiert komplexen 
Wertes in (117) hervor. 

Wird von einem Stromfaden auf ein Emissionsgebiet F, iibergegangen, so 
bedeutet Index v in den Ausgangsgleichungen, zum Beispiel (81), sowohl eine 
Abhangigkeit von U, als auch von x, y, dem Ort des Emissionspunktes. Wird in 
(107) die Stromdichte i(x, y, U,) eingefiihrt, so folgt aus (81), wenn mit 2 mul- 
tipliziert und umgestellt wird, wobei auf das Leitwertersatzbild an zwei Klem- 
menpaaren durch Austausch von %, und %, gegen die KurzschluBstréme J, 
und %, bezogen wird: 


; Riya dF, 31,(%, ¥, Us) h(x, y, Ue) 
(120) 


2 
= [BED ar Sil 9.0) 30.9.0 +2 [Sat SR]. 


Hieraus folgt, dass die durch die Emissionsschwankungen entstehenden Schwan- 
kungsprodukte und -quadrate der KuzschluBstr6me an den Klemmenpaaren 
nur dann durch (120) dargestellt werden kénnen, wenn die Temperatur der 
emittierenden Oberflache T(x, y) konstant ist. 

In Figur 2 ist ein verlustloser elektromagnetischer Hohlraum H angegeben, 
der eine Elektronenstro6mung mit der emittierenden Flache F, und einer ein- 
heitlichen.Emissionstemperatur einschliesst. Drei weitere Hohlraume, zweian 
der Begrenzung von H und einer innerhalb des Raumes, sind angeordnet. Die 
Elektronen sollen ungehindert hindurchgehen kénnen, und die Laufzeiten der 
Elektronen innerhalb dieser drei Hohlréume sollen verschwindend sein. Wenn 
die Elektronen durch H, und H, nur in einer Richtung entweder von aussen 
nach innen oder umgekehrt gehen, ist der Wirkleitwert fiir eine Anregung in- 
nerhalb von H, und H, Null. Die wechselseitigen Leitwerte zwischen H, und 
H, sind Null, wenn die Elektronen durch H, und H, in der gleichen Richtung, 
bezogen auf die gleiche Flachennormale, gehen. Fiir H; sind die Leitwerte im 
allgemeinen nicht Null. Die bei einer inneren Anregung von H, geschwindig- 
keitsmoduliert hindurchgegangenen Elektronen ergeben eine Riickwirkung auf 


RE 


«* we 
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_ die neu eintretenden Elektronen. Ebenso sind die wechselseitigen Leitwerte 


Si, + 3, von H; mit H, oder H, aus dem gleichen Grund im allgemeinen 
_ nicht Null. In den Fallen, in denen die Leitwerte Null sind, fallt in (120) das 
Glied 2 (S,, + GS) weg. Die Gleichung (120) besagt dann, dass die durch die 
Emissionsschwankungen an den inneren schmalen Randern von H, und H, 
_ induzierten wechselseitigen Schwankungsprodukte und -quadrate der Kurz- 
-schluBstréme genau so gross sind, wie wenn in jedem Stromfaden die Schwan- 
kungsstréme der einzelnen Strémungen dieses Stromfadens zusammengefasst 


4 


; 
: 


am Anfang der langsamsten Stré- 
mung wirken wiirden. Lasst man 


_ die seitlichen Abmessungen von 
_ H, und Hy gegen Null gehen, so 


werden die Kurzschlussstréme 


_ genau so gross wie die auf den 


 treffenden 


Querschnitt von H, und H, auf- 
Konvektionsstréme. 
Vergréssert man die seitlichen 
Abmessungen von H, und Hg, so 
dass sie nicht mehr als klein ge- 


_ gentiber der Vakuumwellenlange 


der hochfrequenten Vorgange 
sind, so besteht ein Unterschied 
zwischen KurzschluBstr6men und 
Konvektionsstr6men, da bei nor- 
malen Mediumkonstanten Trans- 


_ formationen in H, und H, auf- 


treten. Besitzen H, und H, ein 


Figur 2 


Elektromagnetischer Hohlraum H mit emittierender 
Flache F,, und drei weiteren Hohlraumen Hj, H5, H3 


mit den Grundflachen Fy, Fs, F3. 


Medium mit der besonderen Eigenschaft einer Permeabilitat Null, so treten 
keine linearen Transformationsverzerrungen auf, und die KurzschluBstrome 
werden auch bei endlicher Flache mit den Konvektionsstrémen identisch. 
Werden die Hohlraume H, und H, entfernt und Flachen / und fF, auf der 
inneren Oberflaiche von H gewahlt, so kann in (120) 3, bzw. 3, durch die auf- 
treffenden Konvektionswechselstréme £, bzw. 2, an /, bzw. J ersetzt werden, 
wobei die Leitwerte Null gesetzt werden kénnen. Fiir die wechselseitigen 
Schwankungsprodukte der auftreffenden Konvektionsstréme an /, und F, gilt 


dann: 


Fe Ueg=00 


ch J 2edile. 9, U) dF Lag(e, 9, Ue) M(x ¥, U) 


Ue=0 


(121) 


Fe 
= [ 2ei(x, 9) dF Lael, ¥, 0) M(x, 9, 0) 
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Das Schwankungsquadrat der Konvektionsstréme an F, betragt: 


Fe Ue=00 
2edi(x, y, U,) dF, &1,(*, y, U,) Li (%, y, Ue) 


U0 


(122) 
Fe 
= [ 2ei(x, y) AE Lig(x, y, 0) L(x, ¥, 0) 


Bei senkrechter Richtung des magnetischen Fiihrungsfeldes zur emittieren- 
den Oberfliche liefert die normale Gliithemission zwanglos ein Exponential- 
spektrum. Bei schiefer Richtung wird der Sattigungsstrom auf cos y (p= Winkel 
zwischen Flachennormale und Magnetfeld) geschwacht. Es ergibt sich das 
gleiche Spektrum, wenn angenommen wird, dass der Bruchteil 1 — cos p der 
Elektronen, die wieder auf die Kathode zuriickfallen, keine Vorzugsrichtung 
der Emission bewirkt. Man bezieht zweckmassigerweise das Verteilungsgesetz 
auf die Richtung des magnetischen Fiihrungsfeldes. Dabei findet man ebenfalls 
ein exponentielles Verteilungsgesetz von der Form 


: P Sie e dU, cos 
di =idF po, apy, exp RE (mo + RT/ce) (RT)? ’ 


wobei die Emissionsenergie W = e U,+ W, ist. (W) Emissionsenergie, wenn 
nur eine Bewegungsgrésse py, quer zum Fiihrungsfeld besteht.) 

Fir eine Str6mung mit diskreter Geschwindigkeit werden noch die Verein- 
fachungen angegeben, die bei verschwindender Raumladung bestehen. Wird | 
bei den Ausdriicken fiir das Leitwertersatzbild (40) und (41) die Differential- 
beziehung (54) benutzt, so ergibt sich fiir die wechselseitigen Leitwerte und den 
Wirkleitwert an zwei bzw. einem Klemmenpaar: 


zr [Sin Be]. 023) 


cae eh a 
Sie + SH = At aU, ne Seal 2 Sy, = Ai 

Der Koeffizient By, aus (43) folgt mit (56) zu: 
Bu, er ane im Ra ° (124) 


Bei den Untersuchungen in diesem Abschnitt I stehen die Anfangsmodula- 
tionen in keinem gegenseitigen Zusammenhang. Im folgenden Abschnitt II 
werden Vorgange behandelt, bei denen die Anfangsmodulationen teilweise oder 
vollstandig abhangig sind von den Endmodulationen anderer Strémungen. Dies 
hat Bedeutung bei dem wichtigen Fall des Vorhandenseins von Potentialminima 
in den Strémungen. 


(Eingegangen: 26. Juni 1953.) Abschnitt IT folgt. 
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Der Quotienten-Differenzen-Algorithmus 


Von HEInz RUTISHAUSER, Ziirich?) 


1. Einleitung 


Im Anschluss an eine praktische Anwendung des BO-Algorithmus (Biortho- 


~ gonalisierungs-Algorithmus von C. Lanczos [4], [5]2) machte mich Herr Prof. 


__E. STIEFEL, ETH., auf das Problem aufmerksam, die h6heren Eigenwerte direkt 
aus den sogenannten Schwarzschen Konstanten zu bestimmen, das heisst ohne 
_ den Umweg iiber die Orthogonalisierung. Auf diese Anregung hin entwickelte 
der Verfasser einen Algorithmus, der die gestellte Aufgabe lést. 

Allerdings gab bereits A. C. AITKEN [1] eine Methode an, welche haupt- 
sachlich zur Auflésung algebraischer Gleichungen gedacht war, aber auch die 
Bestimmung hoherer Eigenwerte aus Schwarzschen Konstanten gestattet. 
Ferner stammt von C. LANczos ein Algorithmus?) zur Bestimmung des charak- 
teristischen Polynoms einer Matrix aus Schwarzschen Konstanten. Uberdies 
entwickelte J. HADAMARD in seiner Dissertation [2] eine Methode zur Bestim- 
- mung der Pole einer durch ihre Potenzreihe gegebenen Funktion. Er hat damit, 
wie §2 zeigen wird, auch das eingangs erwahnte Eigenwertproblem gelést. 
Wenn hier das schon geléste Problem nochmals aufgegriffen wird, so geschieht 
dies deshalb, weil der entwickelte Algorithmus eine Reihe von weiteren An- 
wendungen gestattet und insbesondere auch wertvolle Beziehungen zur Ketten- 
bruchtheorie vermittelt*). 

Die vorliegende Arbeit befasst sich mit den theoretischen Grundlagen, wah- 
rend tiber die zahlreichen Anwendungen (Aufldésung algebraischer Gleichungen, 
Eigenwertprobleme, Umwandlung einer Potenzreihe in einen Kettenbruch und 
umgekehrt, Summation schlecht konvergenter Reihen usw.) spater berichtet 
werden soll. 

Es sei schliesslich noch bemerkt, dass an der GaMM.-Tagung 1953 (21. bis 
24. April in Aachen) auszugsweise tiber diese Arbeit berichtet wurde [7]. 


2. Problemstellung 


Sei A eine m-reihige Matrix, A* ihre Transponierte, ferner seien %, Vg Zwei 
Vektoren, von denen x, in bezug auf die Matrix A und yp in bezug auf die 


1) Institut fiir angewandte Mathematik der ETH. 

2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 251. 

3) Es handelt sich nicht um den BO-Algorithmus, vgl. vielmehr Kapitel VI bei [4] oder 
S. 173—179 bei [5]. 

4) Herrn Prof. Stiefel verdanke ich auch die Anregung zur Vereinfachung einiger Beweise mit 


Hilfe der Kettenbruchtheorie. 
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Matrix A* in allgemeiner Lage sei, das heisst, im Koordinatensystem der Eigen- 

und Hauptvektoren von A (bzw. A*) verschwinde keine Komponente von x9 

(bzw. yy). Ich bilde nun aus % und 4p die unendlich vielen Vektoren x, = A” Xo 

und y, = (A*)” y) sowie die innern Produkte s,,, = (%,, y), die bekanntlich nur 

von der Summe der Indizes ws, vy abhdngen und als Schwarzsche Konstanten der 

Matrix A in bezug auf die Anfangsvektoren x) und y, bezeichnet werden. 
Die mit diesen Schwarzschen Konstanten s,, gebildete Funktion 


co 


Sy 
te) = D) yar (1) 
0 
ist nun rational, und ihre Pole sind die Eigenwerte von A, denn es ist 
2) =(B 44%). Mt B=(zk— Al, 


Damit ist das gestellte Problem, namlich die Bestimmung der Eigenwerte 
einer Matrix A aus Schwarzschen Konstanten, zuriickgefiihrt auf die Bestim- 
mung der Pole einer rationalen Funktion 


f(s) =P 4+ 340+. 


Ist jedoch A eine unendliche Matrix oder ein Integraloperator, so ist die 
Funktion f(z) zwar nicht mehr rational, aber ihre Pole sind immer noch Eigen- 
werte von A. In gewissen Fallen kann man beweisen, dass /(z) fiir z + 0 mero- 
morph ist, beispielsweise dann, wenn die Schwarzschen Konstanten einem 
selbstadjungierten und volldefiniten Eigenwertproblem entstammen, fiir wel- 
ches der Entwicklungssatz gilt. 

Die dargelegten Griinde rechtfertigen es wohl ohne weiteres, nicht das Pro- 
blem der Eigenwertbestimmung, sondern wie J. HADAMARD [2] das etwas all- 
gemeinere Problem der Bestimmung der Pole einer durch eine Reihe (1) gege- 
benen Funktion /(z) zu behandeln, wobei die s,, beliebig komplex sein diirfen. 
Dariiber hinaus wird in gewissen Fallen auch noch die Partialbruchzerlegung 
der Funktion /(z) gesucht sein, falls eine solche existiert. 


3. Der Quotienten-Differenzen-Algorithmus (QD-Algorithmus ) 


Fiir den Fall, dass die erzeugende Funktion 


oo Sy 

/ (2) = a getl 
ausserhalb eines gewissen Kreises | z| = R genau einen einfachen Pol A, besitzt, 
gab bereits D. BERNOULLI die Formel 

Ay = lim 5 (2a) 


v—>0O0 Sy 
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a 


an. Zur Bestimmung weiterer Pole dienen die Formeln (21) bei AITKEN pals 
beispielsweise erhalt man unter der Voraussetzung |A,| >|A_| > |A;| fiir den 
Pol A, die Beziehung 


oa Sy4a Sy42 

| Spe Syis 

j 

‘ : s 

; 44m (2b) 
4 v—>00 Sy Sy44 

Sp44 Sp42 

S 


v 


und entsprechende Formeln kénnen fiir die weiteren Pole aufgestellt werden. 
‘Nun gilt aber die Identitat 


Sy Sy41 
Sy4a Spa! __ ( Spee “+ 
ies coc: | ? 
Sy Syd Sy 


_ der wir entnehmen, dass man A, auch durch folgenden Prozess aus den Schwarz- 
schen Konstanten s, = s{)1) ermitteln kann: Man bildet die Quotienten 

— gf) = st» /s®, dann die Differenzen dieser ¢”: d\”) = qi’*?) — gq”, und schliess- 
lich die Produkte 


Dann ist A, nach (2b) der Grenzwert der Quotienten gy) = s+» /s¥) (fiir »> 0), 
genau so wie A, nach (2a) der Grenzwert von q\) ist. 

Fiir die praktische Rechnung wird man die Zahlen s°’, gq”), .... unter- 
einanderschreiben (vgl. das Beispiel auf Seite 236) ; man spricht dann von einer 
« Quotientenkolonne » gi”) bzw. von einer « Differenzenkolonne » usw. 

Man wird nun vermuten, dass sich durch konsequente Fortsetzung dieses 
Prozesses auch die weiteren Pole A, A,, ... als Grenzwerte entsprechender Quo- 
tienten g¥), g!), ... ergeben (vorausgesetzt, dass |A,| >|A,|>|A3|>|Aa|>---). 

Dies ist in der Tat der Fall, nur muss man die weiteren Differenzenkolonnen 
a) (o = 2,3, ...) vor der Multiplikation mit s!’*”) modifizieren; es ergibt sich 
dann die folgende allgemeine Vorschrift zur Aufstellung des QD-Schemas zu 
den gegebenen Reihenkoeffizienten (bzw. Schwarzschen Konstanten) s, = s{”: 

Ausgehend von den s‘”) berechne man nacheinander (fiir ‘o = 1, 2, 3, ...) 


s@) = sth qo) (= aA Sy Spa 
v 


Sy44 Spy 2 


git = s@tB js) Quotientenreihe der s”) , 

d” = q°*!)_ q” Differenzenreihe der q”), B) 
e”) = d+ e?*" modifizierte Differenzen?) , 

RR Phe ae e”) neue s-Reihe . 


1) Wir bezeichnen die Schwarzschen Konstanten s, fortan auch gelegentlich mit sf?) 
2) Fiiro = 1 stimmen die el?) mit den al”) iiberein, also ist immer el”) = 0 zu setzen. 
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Falls man die s”)-Werte (« > 1) nicht benétigt, kann man auch mit der 
folgenden komprimierten Vorschrift arbeiten (ausgehend von ée”) = 0 und 


qn = Sy4a/Sy)* 


a) = a 4 go — of, 
(+3) qa 2 ee 
bis Loti el ( ) (4) 
Vo41 Pas G5 (P) et 
te 


Beispiel zum QD-Algorithmus 
Die s¥) seien die sogenannten Fibonaccischen Zahlen ; es ist dann 


z 


ze—z-—1° 
a ee oe eee ee wey 
i 
1 1 1 
2 -1 
2 —0,5 -1 +0,5 0 0 
Us —0,5 
3 0,166 667 0,5 — 0,166 667 0 0 
1,666 667 — 0,666 667 
5 — 0,066 667 — 0,333 333 + 0,066 667 0 0 
1,6 — 0,6 
8 + 0,025 +0;2 — 0,025 0 0 
1,625 — 0,625 
13 — 0,009615 — 0,125 
1,615 384 
21 J | 
A = 1,618034 Ay & — 0,618 034 


An diesem einfachen Beispiel tritt bereits eine Erscheinung zutage, die sich 
als allgemeingiiltig erweisen wird: 


Wenn f(z) eine rationale Funktion ist, deren Nenner den Grad n besitzt, so 
bricht das QD-Schema nach der n-ten modifizierten Differenzenkolonne ab, ganz 
ahnlich, wie das gewohnliche Differenzenschema eines Polynoms nach der n-ten 
Differenz abbricht. 


Diese Eigenschaft erlaubt uns, dasQD-Schema in einem solchen Fall auch 
von oben nach unten fortschreitend aufzubauen, wenn die in der obersten 
Schragzeile stehenden Werte 


0) (0). (0 0 0 
Soi 1 Gee aaa aa 
bekannt sind. Dieses Vorgehen, welches einem wohlbekannten Prozess beim 
gewohnlichen Differenzenschema vdllig analog ist, wird vor allem fiir die An- 
wendungen von grosser Bedeutung sein. 


TN) Sele 


5 
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4. Die zugeordneten Polynome p") (z) 


Gewissermassen als Vorlaufer des BO-Algorithmus hat C. Lanczos das 
bereits erwahnte Verfahren zur Bestimmung des charakteristischen Polynoms 
einer Matrix aus Schwarzschen Konstanten angegeben. Im Verlaufe dieses 
Prozesses, der tibrigens mit dem QD-Algorithmus nahe verwandt ist), wird 


_ eine Kette von Polynomen ~,)=1, ?,(A), f(A), ... aufgebaut, deren letztes 


- Glied #,,(A) das Minimalpolynom der Matrix A ist. 


Entsprechend bilden wir hier das zweidimensionale Schema der Polynome 


_ pl)(z) vom Grade o (P-Schema) : 


(0) 
0 


(0) 
: (0) 
1 
pt 
(1 (0 
5 ) PS ) 
2 (1 
po” Oe 
(2) (1) 
1 : 3 


das wir nach der folgenden Regel aufbauen: 
py =1 fir aller, 
PM (2) = 2 PPP) — af) O92) (v= 0,1,...5 6=1,2,...). 


Die Regel (5) kann auch schematisch dargestellt werden, so dass man sofort 
ersieht, in welcher Weise ein Element im P-Schema aus den links davon ste- 
henden Elementen p+) und p” , gebildet wird: 


(5) 


Oe 


ty see 
ao—i1 . 


Wir wollen nun fiir diese Polynome #'”)(z) die beiden folgenden Relationen be- 
weisen : 


o-1 


pag) 0 = 0. 
ey tl? +e ae Poy 9 ft oo > 0. 


(v+1 Xs) eee AE Aas) ts 6 
fez— piv) — —el) pl fir o 20 (6) 


(7) 


1) Vergleiche hierzu § 8. 
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Beweis von (6) 


Wir verwenden als Induktionsvoraussetzung die Formel (6) mit o — 1 an 
Stelle von o, also die Beziehung 


(ued Gast) 
pet = dk th pity =.0, 


welche fiir o = 1 wegen p\ = 1 und ef) = 0 sicher fiir alle y erfiillt ist. 
Wir bilden nun die Grésse 


ic p° (v+1) — p” 4 el”) pee J 
die wegen (5) zu 


z[por?) — +h] — g?*) ) poh) =e q\”) it el?) per) 
wird. Unter Verwendung von 


gery) ES el”) 22 qi” a 
[siehe (4)] ergibt sich 


o=2(p%" iene Pe (oy) — q? [per (v+1) — pi”) i}t et Sate : 
oder unter Verwendung der Induktionsvoraussetzung: 
ee 4 ert) pb Mee 1a qe (>) ee pas 1) a el? font) per), 
was sich schliesslich unter Verwendung von 
g'”) el”) =o ee 1) 


[siehe (4)] zu 
6 <3 elt) {—z z pre) Bs "hie pith of porn 


reduziert. Dieser Ausdruck verschwindet aber wegen (5); was zu beweisen war. 


Beweis von (7) 
Unter Verwendung der mit z multiplizierten Gleichung (6) ergibt (7): 
pi”, — 2 perry) pe ze”) porn “ qo, MW as el) p”) ae el) q”) p”, 
=P ehe. Oe hela ee be eee ee Pele 


Wegen (5) verschwinden aber die eingeklammerten Ausdriicke ; was zu beweisen 
war. 


eae LS, Viele ae 


a 
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Als Beispiel sei hier das zum QD-Schema auf Seite 236 gehorige P-Schema 
angegeben : 
i pt pe? 
1 
z—1 
ui z2@—zg—] 
z—2 
a z2—z—1 
z—1,5 
1 z2—zg—] 
z — 1,666667 
ik 


Man erkennt aus diesem Beispiel sofort eine weitere grundlegende Eigenschaft 
des QD-Algorithmus, die wir spater (§ 7) beweisen werden: 


Ist co & 
f(z) = oF et 
0 


eine rationale Funktion mit dem Nenner N(z) = 2” + ---, so steht im zugeordneten 
P-Schema in der Kolonne p\? fiir alle v dasselbe Polynom, némlich N(z). 


5. Beziehungen zur Kettenbruchtheorie 


Wie die Formel (7) zeigt, sind die Polynome #¥)(z) fiir festes » gerade die 
Naherungsnenner des folgenden J-Kettenbruchs?) : 


ee iN ate ie ON See 
v 


fag [e-@-P  -P -@ 


(8) 


Dieser ist wiederum der sogenannte « gerade Teil » des nachstehenden S-Ketten- 
bruchs, in welchen er durch Extension transformiert werden kann?) : 


le ig Ss 
rts had | 2 ier = ee. aa (9) 


Bildet man nun durch Kontraktion den « ungeraden Teil» des Kettenbruchs (9), 
so resultiert offenbar: 


as a) | a) qe) | 


fe) = Se {1+ —" 


2 
pf) p-@-—  fe—a— 


1) Vgl. H. Watt [8], §§ 50 und 51. 
2) Vgl. H. Watt [8], § 96, insbesondere Formeln (96.6), (96. 12) und (96.13). 
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oder wegen (4): 


if 1 
sig Dione eelipi eee 
| = qv*h) |z eS qu ee he el? +) 
ars ge ay (2) 
es geri — e+) oo BAN a 


Es ist also fiir alle v 


je) = te + Heal 


und damit offenbar fiir |z| > lim sup |s,,1/s, |: 


oo 


fle) = SP peat (10) 


Das bedeutet insbesondere, dass die Formeln (8) und (9) fiir » = 0 Ketten- 
bruchentwicklungen der urspriinglich gegebenen Funktion /(z) darstellen; der 
QD-Algorithmus leistet somit auch die Kettenbruchentwicklung einer durch 


foe) 
eine Potenzreihe )’ s,/z*+1 gegebenen Funktion. 
0 


Beispielsweise ergibt sich aus dem QD-Schema fiir die Fibonaccischen Zah- 
len (Seite 236) der endliche S-Kettenbruch 
Yo A Hg A 
are MES te 


fiir die Funktion 2/(z? — z—1). 

Auf Grund dieses Zusammenhanges lassen sich sofort gewisse Ergebnisse 
aus der Kettenbruchtheorie auf den QD-Algorithmus iibertragen. Bae 
weise sind die im J-Kettenbruch (8) als Teilzahler auftretenden Gréssen e!) gq’? 
mit den Reihenkoeffizienten der durch den Kettenbruch dargestellten Funktion 
(10) in folgender Weise verkniipft?) : 


(y) (v) 
el) gi) Ho+1 Het | 
ao to (H' (2 . 
Hierbei ist H® =1 und 
0 
Ss, S, +1 Syyo-1 (11) 
S S Ss 
(yy v+1 y+2 v+o 
Vs be) 
Syi¢ a Syic vy+2o0-2 


1) Vgl. § 42, insbesondere Formeln (42.3) und (42.7) bei H. Watt [8]. 


Vee Th Nee Se ee 


4 Damit wegen (11): 
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_ Ferner ist ja nach (4) 


o-1 
v J! v v 
Sy = Cony Soon = {ee GI }s0), oder si) = sh) TT (el) gi) . 
wal 


und gq?) = —* se (12) 


Aber auch die Naherungsnenner #'”)(z) des Kettenbruchs (8) kann man explizite 
durch die Reihenkoeffizienten ausdriicken. Weil namlich der o-te Naherungs- 


__bruch von /,(z) (dessen Nenner 9) ist), bis zum Glied 1/z®% einschliesslich mit 


f,(z) ibereinstimmen muss (vgl. [8], §51), miissen im Produkt /,(z) p2)(z) die 
Glieder mit 


samtliche verschwinden. Daraus folgt aber wegen 


f,(2) = Sy + Sy41 ae Sy42 ae Pike 


z Z2 23 


und weil £(z) mit 2° + --- beginnt: 


S, Sy4d SW eal I 
(v) 1 Sy41 Sy+2 nee Spe & (13) 
£2) = 
ty 5 HORE Sa OI oS Gaceteeo eterna Fe 
g ge 
Syio Syto41 y+2a0-1 


6. Schwierigkeiten bei der Bildung des QD-Schemas 

Die Quotientenbildung im QD-Schema wird natiirlich gestért, sobald einer 
der Werte s'”) verschwindet; man muss aber zwei verschiedene Falle unter- 
scheiden: 

a) Es sei s?) =-0; aber, s+) + 0 

Dann wird q'”) = oo, so dass das QD-Schema in einem von diesem Element 
aus nach rechts 6ffnenden Winkelraum nicht definiert ist. Man kann leicht fest- 
stellen (vgl. [8], §50), dass dann auch der J-Kettenbruch (8) fiir die Funktion 
/,(2) nicht existiert. 


b) Es sei s”) = 0 fiir alle » => us 


In diesem Fall kann man das QD-Schema sinnvoll nach rechts und unten 
fortsetzen, indem man fiir die sich zunachst als 0/0 ergebenden Quotienten gq’ 


ZAMP V/16 
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(v = w) konsequent den Wert 0 einsetzt. Es erscheinen dann im QD-Schema 
rechts und unterhalb des Elementes s\“) (welches im nachstehenden Beispiel 
eingerahmt ist) nur noch Nullen. Fiir den Kettenbruch (8) der Funktion 
},(z) (v = pw) bedeutet dies, dass die Naherungsbriiche R,, R,,1 usw. sich nur 
noch formal durch Faktoren 2 im Zahler und Nenner unterscheiden und nach 
dem Kiirzen alle gleich sind. /,(z) ist dann also rational und damit auch f(z). 

Das identische Verschwinden der e?) und g) rechts und unterhalb einer 
gewissen Stelle wirkt sich entsprechend auch auf die Polynome #”) aus. In der | 
Tat muss nach (6) f)(z) von v unabhangig sein, wenn alle e”) verschwinden; 
ausserdem wird p”) = z p+), wenn q”) = 0 ist. 


ead FE: z2+2 
Beispiel: QD-Schema fiir f(z) = > @+1) 
sv) =s, g”) el?) si?) qs) aly) el?) si”) qi”) av) el”) s(?) 
if 
-1 
-1 —2 2 
—3 3 
eo +2 6 —3 —1 —6 
—1 == 0 0 
6B 0 fo | 0 0 0 0 Lo Pe 
-—1 ol 0 0 
+3 0 0) (0) 0 0 
Zugehiriges P-Schema 
po p) py py? p”) 
il 
z+1 
1 : 22-3 
z+3 2® (g+1) 
1 2 (z+1) 23 (z+1) 
. z+1 z2 (¢+1) 
1 2 (+1) 


7. Grundlegende Eigenschaften des QD-Algorithmus 


Wir wollen hier die bereits in §3 auf Grund eines numerischen Beispiels ver- 
mutete Eigenschaft des QD-Algorithmus beweisen. Zundchst folgt aus § 6 
unmittelbar : 


Satz 1. Verschwinden die n-ten modifizierten Differenzen e) im QD-Schema 
> fiir alle v, so ist 
4 Ss 
f(z) = a par 


eine rationale Funktion, deren Nenner hichstens den Grad n hat. 
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Es gilt aber auch umgekehrt: 


Satz 2. Wird mit den Koeffizienten S, einer rationalen Funktion 


lee) 


i) = 2) sr 


deren Nenner N(z) den Grad n hat, das QD-Schema aufgestellt), so verschwinden 
die n-ten modifizierten Differenzen e\”) fiir alle v. 

Ferner hat das zugeordnete P-Schema die Eigenschaft, dass die Polynome 
_ pl)(z) fiir alle v gleich dem Nenner N(z) der Funktion f(z) sind. 


Die Voraussetzung von Satz 2 ist insbesondere dann erfiillt, wenn die s, 
Schwarzsche Konstanten einer endlichen Matrix A in bezug auf beliebige 
Anfangsvektoren sind. 


Bewets von Satz 2 


Ist /(z) rational und regular an der Stelle z = 0, so sind auch die /,(z) rational, 
und zwar alle mit demselben Nenner N(z). Nun ist aber der J-Kettenbruch 
einer rationalen Funktion, falls er iberhaupt existiert?), immer endlich, und die 
Folge der Naherungsnenner desselben endet mit dem Nenner N/(z) der ratio- 
nalen Funktion. Im Fall der Funktion /,(z) heisst dies, dass die Folge der Poly- 
nome #!”)(z) fiir jedes feste v wie folgt beschaffen ist: 


PP =1; PMR), APE), -. 2712), BQ) = NE). 


Damit ist der zweite Teil des Satzes bereits bewiesen, und der erste Teil 
ergibt sich anschliessend aus (6) mit o = n, weil ja die Polynome p” , = 2"~1+4 --- 
sicher nicht identisch verschwinden. 

Ist hingegen z = 0 ein r-facher Pol von /(z), so ist f,(z) eine rationale Funk- 
tion, die an der Stelle z = 0 regular ist und deren Nenner den Grad m — ¢ hat. 
Somit folgt aus dem ersten Teil des Beweises, dass el”), = 0 fiir » 27 (und 
ebenso s'?),_,= 0). Nun kann man offenbar das QD-Schema in der in §6 
unter b) angegebenen Weise erganzen und findet, dass bereits fir » 27 —1 
2) 4 ~0 ist, ferner e?) ,=0 fiir y =>7—2 usw, bis schliesslich el”) — O fiir 
alle y > 0; was zu beweisen war. 


Eine weitere wichtige Eigenschaft des QD-Schemas, die sich insbesondere 
auch zur Bestimmung der Nullstellen eines Polynoms eignet (vgl. die angektin- 


1) Falls das QD-Schema existiert. Ist z = 0 ein Pol von f(z), so gilt Satz 2 nur dann, wenn das 
QD-Schema gemass der in § 6 unter b gemachten Bemerkung bis zur Kolonne a) erganzt wird. 
2) Nichtexistenz dieses Kettenbruchs fiihrt nach [8], §§ 41 und 42, sofort auf H) = 0 fir ein 


gewisses o < ” und damit wegen (12) aus s(?) = 0. Es muss damit Fall § 6a vorliegen, denn § 6b 
wiirde die Existenz nicht storen. 
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digte Arbeit iiber die Anwendungen des QD-Algorithmus) ergibt sich fast un- 
mittelbar aus der erwahnten Arbeit von AITKEN: 


Satz 3. Ist die Funktion f(z) rational und haben thre Pole verschiedene Abso- 
lutbetrage: | A, | > |4,| > |As| > +++ so gilt fiir das QD-Schema der Koeffi- 
zienten von f(z): 

Mink 2° =O Nea ere (14) 


o 


lim g” =A, 


yp—> Oo y—> Coo 


Beweis. Nach AITKEN gilt unter den gemachten Voraussetzungen: 


Ht) 
Aphe ncn lim, es 
Schreibt man dies fiir m =o und m =o —1 an und dividiert, so folgt wegen 
(12) sofort die erste der Beziehungen (14). Fiir den Beweis der zweiten Behaup- 
tung betrachten wir den Quotienten 


ee F G+ 
oS SEY 
el”) gi? ) 


der fiir y +00 gegen 4,,,/A, konvergiert. Also strebt e wie (A,.3/A,)” gegen 0; 
was zu beweisen war. 


Die Aussage von Satz 3 kann wie folgt verallgemeinert werden: 
Satz 4. Es kénne die Funktion 


in folgender Weise zerlegt werden: f(z) = R(z) + q(z). Dabei sei R(z) eine 
fiir |z| << R rationale Funktion mit dem Nenner N(z) = 2" + --- und (2) 
eine nicht notwendigerweise rationale, aber fiir |z| = R (@ < 1) regulare 
Funktion. 

Dann gilt, sofern das QD-Schema existiert und insbesondere die Polynome 


p) (2) (o.= 0,1, 2,..., 2) fiir unendlich viele Indexwerte y; definiert sind: 
lim p)(z) = N(z) . (15) 


I 


Der Beweis dieser Behauptung kann ohne weiteres der erwahnten Arbeit 
von HADAMARD [2], §§ 14 bis 21, entnommen werden; allerdings wird dort eine 
Funktion /(z) betrachtet, die nach aufsteigenden Potenzen von z entwickelt ist. 
Die bei HADAMARD auftretenden Polynome 


B= 14+ AD 24 AP Go + AW a, 
PB, = 1+ AM e+ AP) 224 0. 4 A of | 


Ee ee 
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deren Koeffizienten durch die Gleichungen (18), (18’) bei [2] definiert sind, ent- 


_ sprechen den in dieser Arbeit verwendeten Polynomen Dy” (z), pO”(z) usw.; 
_ wobei die Indizes #, g, 7 so bestimmt sind, dass 


[41 | =|42| = sent Ase 0 |Apaa|=|Ao8 |= = 
=|Aq|= @ > |Agua|= SBN aetna 


HaADAMARD beweist dann, dass diese Polynome §,, und §;,, fiir m > oo gegen 
Polynome $$ bzw. ‘f’ konvergieren, deren Nullstellen genau die ersten p bzw. 


- die ersten q Pole von f(z) sind, und zwar ist die Konvergenz etwa so gut wie die 


einer geometrischen Reihe mit dem Quotienten 0'/o bzw. 0”/o’. 


Dem Satz 4 kann man folgende ances entnehmen, die wir hier ohne 
Beweis angeben: 


Folgerung 7. Fir die Giiltigkeit der Aussage (14) (Satz 3) fiir einen einzelnen 
Indexwert o=k ist bereits die Bedingung |A,_,| > |Ax| > |Agay| hin- 
reichend. 

Folgerung 2. Baut man mit den Werten gq’) und e) des QD-Schemas 
analog zu (5) die Polynome #)(z) nach folgender Vorschrift auf: 


B)(z)=1 fir alle v, 
; | (16) 
Be) = 2 BLM) — a B42) (C= m+ 1 s.0), | 
so gilt unter der Voraussetzung 
| Am | > | Ama | 2 | Am+e | aa 2 |4,| > | Ansa | Pe 
die Beziehung 
Lim Bl(e) = (e — Ama) (2 — Amas) «+» (2 — An) (17) 


8. Beziehungen zum BO-Algorithmus von C. Lanczos 


Wie bereits in §4 erwahnt, stammt von Lanczos ein Verfahren zur Bestim- 
mung des charakteristischen Polynoms einer Matrix A aus Schwarzschen Kon- 
stanten s, (bei LANczos mit c, bezeichnet). Es werden dort nach einem bestimm- 
ten Verfahren aus den s, die Polynome fy = 1; (2), f2(2), ... gebildet, deren 
letztes, ,,(z), das Minimalpolynom von A ist. 

Aus dem Gleichungssystem (48) bei [5], dem die Koeffizienten des Polynoms 


bolt) = 2° + gat + mt + Mo 


geniigen miissen, geht hervor, dass p,(z) ebenfalls Naherungsnenner des J-Ket- 
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tenbruchs von 


sein muss, mit andern Worten: 


Die von Lanczos konstruierten Polynome #,(z) stimmen mit den zuge- 
ordneten Polynomen #')(z) des QD-Schemas der c, iiberein. 


Genauere Betrachtungen zeigen ferner, dass das algorithmische Schema 
von Lanczos!) mit den hier verwendeten Bezeichnungen wie folgt ausgefullt 
werden kann: 

0 0,5 t iA) Z P58} 


= (0) (1) (0) (1) (0) (1) 
i Sy Sy Sy Sy S3 Sy 
0 0 (0 (0 (0 (0) 
a= |-a | -| - | | =| -4 
So i 0 0 0 0 0 
Ss i 0 0 
Pee Ale et ers Ve eee (18) 
So 0 
[| q q rs] 
2e|8s|2s|25 
BS. 8S. Sc. 5 San 
ISS Si So tS 
— — a a 
es a oe eee aes ee ee es 
OTs Oo $5 Oe &, os 
| i i M4 MM 


Lanczos hat ferner gezeigt, dass dieselben Polynome #')(z) auch auftreten, 
wenn man die Vektoren %, V9, die zur Erzeugung der Schwarzschen Konstanten 
s, der Matrix A dienten, als Anfangsvektoren fiir den BO-Algorithmus ver- 
wendet?). Da diese Polynome ferner mit den Elementen «,, 6, der Kodiagonal- 
matrix durch die Relationen 


Py (2) — (2 — %) py = 0, 
Peril) — (2 = %41) Pz) + B, PP) = 90 (6 = 1,2,...) 


verkniipft sind, folgt durch Vergleich mit (7) sofort eine Beziehung zwischen 
der mit dem BO-Algorithmus erzeugten Kodiagonalform der Matrix A und 
den mit denselben Anfangsvektoren x, y mit dem QD-Algorithmus erhaltenen 
Werten ¢) und el): 


Oy OT Oa Opt Col iB Go Ce AOR Lies.) ee 


1) Tabelle 1, Seite 263, bei [4] oder Seite 176 bei [5]. 
2) Diese Polynome sind in der Arbeit [6], Seite 41, mit Pg (A) bezeichnet. 


NO ee Pe ee eee 


ee De 


> Sa eee 
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Diese Formeln vereinfachen auch die Schreibweise des J-Kettenbruchs (8) der 


- Funktion f(z) zu 


Nace ets a as 


rea tes a.er det, 
Wie man leicht einsieht, haben nicht nur die Gréssen 


gq +e, und g® © 
eine bestimmte Bedeutung fiir den BO-Algorithmus, sondern auch 


gag ad, pa ge (¢ = 1,2). »=0,1,...). (20) 


Es sind dies namlich die Kodiagonalelemente, die man erhalt, wenn man den 


BO-Algorithmus mit den Anfangsvektoren A*x und (A*)-*y1) an Stelle von 
x, y beginnt. 


Demnach enthalt das OD-Schema soviel Information, wie man durch viel- 
fache Anwendung des BO-Algorithmus erhalt; wir werden in den Anwen- 
dungen darauf zuriickkommen. 


9. Beziehungen zum cg-Algorithmus?) 


Die Beziehung, die nach [3], § 18, zwischen dem cg-Algorithmus und den 
Kettenbriichen besteht, lasst auch eine Querverbindung zwischen cg-Algorith- 
mus und QD-Algorithmus vermuten: 

Sei A eine reell symmetrische Matrix; mit dieser bilden wir einerseits die 
Schwarzschen Konstanten s, fiir die Anfangsvektoren x und y = x sowie an- 
schliessend das QD-Schema. 

Andererseits fiihren wir aber mit demselben Anfangsvektor 7) = x den 
cg-Algorithmus durch (Formelgruppe 5.1 bei [3]), allerdings ohne das Glei- 
chungssystem aufzuldsen. Die dadurch erzeugten Residuenvektoren 7, geniigen 
nach den Formeln (5. 8) und (5. 9) bei [3] den Relationen 


a 1 
RS ass Ser aT | 
Vig) ag ee eget ( tafe 1.) Pee FP  S(j ait 2 (21) 
ta? a; t+1 a; aly a Be 7 i-1 pe aty Pie) ey 
(7,,7;) = 0 fir + +7. 


1) Die Wahl von x beeinflusst dann zwar das Biorthogonalsystem, nicht aber die Kodiagonal- _ 


_ matrix. 


2) cg-Algorithmus = Synonym fiir: Verfahren der konjugierten Gradienten von STIEFEL 
HeEsTENES (siehe [3)]). 
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Dieselben Vektoren 7; liefert aber auch der BO-Algorithmus in der Form 
(2'), (3’) bei [6], wenn man fiir die dort auftretenden s;,,,-Skalare — 1/a,,_, ein- 
setzt ; man erhalt dann (wegen A = A*): 


Ye an 
Ao 

Axi,=— Kou hy Xe Von Xg-1 (o = 2, 3, +++), 4) 
o—1 


(%3, %;) = (0) fiir t=]. 


Vergleich von (21) und (22) liefert die folgenden Korrespondenzen zwischen den 
beim cg- und BO-Algorithmus auftretenden Grdssen 


Mes ft S 1 be-1 
He Gr a= Gan Oot ‘ips (Res 
b b (23) 
o—1 o—1 
=— oder =<, 
Yo Ag1 Ba a5-1 


Ein weiterer Vergleich mit (19) ergibt nun endlich die gewiinschte Beziehung 
zwischen den q, e des QD-Algorithmus und den Gewichtskoeffizienten a, b des 
cg-Algorithmus: 1 b 

gal = 5 OO ee SE ag See (24) 


Ag_1 Ag_1 


10. Die Partialbruchzerlegung der Funktion f(z) 


Nachdem wir Methoden zur Bestimmung der Pole einer Funktion /(z) ent- 
wickelt haben, wollen wir nun versuchen, auch die zugehGrigen Residuen zu 
ermitteln. Zur Lésung dieser Aufgabe gehen wir von der Formel (13) aus, 
welche dasselbe sagt wie 


Satz 5. Es ist 
RO 
p\(2) f,(2) = Polynom + + ---, (25) 


ot 
C4 


das heisst der o-te Naherungsbruch des Kettenbruchs (8) weicht in erster An- 
naherung um den Betrag s!”) ,/227*1 von f,(z) ab. 


Beweis, In der Tat ergibt sich ja das Verschwinden der Koeffizienten von 


aus §5, wahrend der Koeffizient von 1/z°** aus folgender Betrachtung hervorgeht : 


Man erhalt den Koeffizienten von 1/z“*! im Produkt p''(z) f,(z) | 
dadurch, dass man in der Formel (13) die Potenzen z* durch die } (26) 
Werte s, 4... ersetzt. 
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Daraus folgt — mit 4 < o — sofort wieder das Verschwinden der Koeffizien- 
ten von eh os 1 


Ve ee EG. 


aber fir 4 = o wird (13) mit der Substitution (26) offenbar zu H”) ,/H), was 
wegen (12) gleich s’) , ist; was zu beweisen war. 
Fir die Bestimmung der Residuen von f(z) benétigen wir nun die Formel 


ais fir o+o0, 


sai f 16) 1) 2%) (27) 


| 3, fick, o =o; 


in welcher C irgendeinen Weg bedeutet, der alle Singularitaten von f(z) um- 
schliesst. Diese Formel folgt sofort aus (25), denn fiir @ < o fehlt in der Ent- 


wicklung von He) BP (2) BO) 


das Glied mit 1/z; aus Symmetriegriinden folgt dasselbe fiir @ > o, und fiir 
o = g ist der Koeffizient von 1/z tatsachlich s?) ,. 


Nehmen wir nun an, /(z) sei rational und habe  einfache Pole: 


n 


a (28) 


1 


so erhalt man durch Einsetzen in (27) 


n 0 fir o+o, 
y u (0) (0) = 
= C, Ps (A,) i (A,) ia 0) 


s! 


Leet La OC. 


Diese Gleichung kann in Matrizenform geschrieben werden; definiert man 
namlich die Matrizen 


Cy 3; 
Aa) Sg U 
C= reg as 
0) 0 
e, si) 
t (29) 
PP (A1) PY (Aa) PY) (An) 
POA) (Ap) P(A) 
P= 
() (Ay)! 4(Aa) 1) (An) | 
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so lautet sie PC P* = S. Daraus ergibt sich endlich die Formel fiir die gesuchten 


Residuen c,: ° 
GAS Pus £, (30) 


oder ausfiihrlich geschrieben: 
OW Se ftir: year, 


n—-1 1 ni 
Sa PA OU) =) (30’) 
o=0 So+1 a fey Wir 
Da aus (28) sofort folgt, dass fiir alle » 
3 = Os oA, (31) 
t=1 


ist, kann man die Formeln fiir die Bestimmung der c, auch gebrauchen, um 
aus den Werten s(0), s(1), ..., s(2 m — 1) eine Funktion 


s(x). = os EuAC 
it 


zu konstruieren, die an den Stellen vy = 0, 1, 2, ..., 2 — 1 die vorgeschriebenen 
Werte s(v) = s, annimmt (Interpolation durch Exponentialsummen). 
Beispiel. Gesucht eine Funktion s(x) = c, Aj + cy 23, die die Werte 


s(0) = —1, s(1) = 2, s(2) = 16 undys(3) = 68 


annimmt. Wir stellen das QD-Schema auf: 


x a a st d 
—1 
=—2Z 
yo 10 20 
8 —3 
16 —3,75 —60 
4,25 
68 


Hier ist 
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und damit 

1 

z 0 

C-1= P*S-1pa— 
il 

OMee 

also 


¢,=4, C=—5 oder s(x) =4-3%—5-22, 
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Summary 


The quotient-difference (= QD) algorithm developed by the author may be 
considered as an extension of BERNOULLI’s method for solving algebraic equa- 
tions. Whereas BERNOULLI’s method gives the dominant root as the limit of a 
sequence of quotients qi” = a [s formed from a certain numerical sequence 


si), the QD-algorithm gives (under certain conditions) all the roots A, as the limits 


of similiar quotient sequences q”) = ide i Close relationship exists between 
this method and the theory of continued fractions. In fact the QD-algorithm 
permits developing a function given in the form of a power series into a con- 


tinued fraction in a remarkably simple manner. 
In this paper only the theoretical aspects of the method are discussed. Prac- 


’ tical applications will be discussed later. 


(Eingegangen: 5. August 1953.) 
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723 A/B Reflex Klystron Performance 


By Giacinto Zito, Turin?) 


1. Foreword 


During a systematic series of measurements carried out on a 723 A/B 
Reflex Klystron it was possible to outline the results in a number of diagrams 
which proved to be of a considerable practical value inasmuch as they make 
it possible to follow, at a given load, the behaviour of the tube during the 
simultaneous variations of the applied voltages to the anode and the reflector. 
The practical value of these diagrams is evident in that they enable one to obtain 
very pure amplitude and frequency modulation and to reach conditions of 
very stable operation. 


2. Klystron Characteristic Curves and Methods Used in Obtaining 
Them 


Consider the circuit shown in Figure 1. Through varying simultaneously 
the anode (resonator) voltage V, and the reflector voltage V, so as to keep 


Crystal 


Meter 


Frequency 
Meter 


Terminator 
VF 


Attenuator 
SE a 


i 
Directional Coupler Thermistor 
Mount 


Regulated 
Power 
Supply 


Figure 1 
Circuit for the measurement of static klystron performance. 


constant the frequency generated by the klystron tube, a nearly linear relation 
between V, and V, is obtained. In other words the ratio between the variations 


1) Istituto Elettrotecnico Nazionale « Galileo Ferraris». 


NN 


aS 
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; of the anode voltage and the reflector voltage along a constant frequency 
_ curve is nearly constant. 


a 


The approximate linearity of the curves is to be interpreted as a conse- 


quence of the linear variation of the potential in the drift space which would 


be obtained exactly in the case of parallel flat electrodes and in the absence 
of space charge. 


280 
270 
260 
250 
240 


230 


Frequency Deviation Mc 


220 Per cent Power Output 


SWR 1 
210 100% Power Output(mw) 30 


Frequency (Mc) 9369 


200 


| iter 
-—140 ¥ - 130 


- 150 


- 160 


Figure 2 
Performance diagram, first mode. 


Slight deviations from linearity may be due to the different geometry of the 
electrodes, to the effect of the space charge, also to the fact that in general 
along a curve of constant frequency the dissipated power is variable. It follows 
that a variation of the operating temperature in the tube may modify its 
geometry with consequent frequency variations. 

A number of these curves at different frequencies are shown in the diagrams 
of Figures 2, 3, and 4, corresponding to three operating modes of the klystron 
tube. 

Simultaneous power measurements were carried out by means of a thermi- 
stor and a direct reading bridge (see Figure 1). On the diagrams of Figures 2, 3, 
and 4 a series of constant power curves have been sketched. In the diagrams 
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Figure 3 Figure 4 
Performance diagram, second mode. Performance diagram, third mode. 


the curve A/ = 0 refers to the frequency of 9369 Mc. The frequency difference 
between two adjacent curves is 5 Mc. 

The output power (30mW) at the point V, = 300 V on the curve Aj = 0 
(see Figure 2), has been taken as the reference value of 100° for the constant 
power curves. 

Since the thermistor mount does not produce appreciable reflections across 
the entire band, the load is constant and equal to the characteristic impe- 
dance of the waveguide. 

The constant power curves are similar to a series of parabolas, with power 
values increasing from the outside parabola to the inside ones. The point 
marked 100% in the diagram of Figure 2 would be the vertex of the parabola 
100°% power which would develop wholly above the 300 V abscissa, but which 
is not drawn in the diagram. 

At the ends of the constant frequency curves the power diminishes quickly 
to zero. 
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3. Examples of Utilization 


3.1. Amplitude Modulation 


If amplitude modulation is produced by feeding only one of the tube elec- 
trodes, as is usually done, undesired frequency modulation cannot be avoided. 
For this reason, it is often necessary to resort to some method of pulse 


_ modulation. In some cases, however, especially in the measurement field, the 
_ large frequency spectrum inherent to pulse modulation may cause certain 
_ difficulties which limit the utility of this system. In these cases, a considerable 


advantage can be obtained if the amplitude modulation is produced by modu- 
lating anode and reflector electrodes at the same time, utilizing a region on 
the diagrams where the power varies considerably along one of the constant 
frequency curves. The diagram of Figure 2 presents the wanted requirement 
in the region between the frequency curves Af = 0 and Af = +5 and the power 
curves P = 50% and P= 100%. 

In order that the point of operation should stay on a constant frequency 
curve during the modulation process, it is necessary to vary simultaneously the 
anode voltage and the reflector voltage so that, at every instant, the ratio 
between their peak variations AV, and AV, is constant and equal to the slope 
of a curve V, versus V,, as shown in Figure 2. The amplitude and phase relation 
between voltages V, and V,.can be easily obtained by means of two transformers 
having the secondary windings in series with the anode and reflector circuits. 
It may be necessary to limit the maximum value of the depth of modulation to 
about 50%, in order to avoid distortion. However, this condition does not 
limit the practical usefulness of the system for measurement purposes. 


3. 2. Frequency Modulation 


The frequency modulation is generally obtained by the superposition of a 
sinusoidal or saw-tooth voltage on the bias voltage of the reflector. The fre- 
quency variation will then depend on the amplitude of the modulating voltage. 
However, a certain amount of amplitude modulation will always accompany 
the desired frequency modulation. 

From Figures 2, 3, and 4, it appears possible to obtain frequency modula- 
tion free from amplitude modulation, if the nearly straight parts of the con- 
stant power curves are employed where they cross a great number of constant 
frequency curves of uniform spacing. In view of the shape of the curves, the 
further from the vertices the more favourable are the regions for such a system. 

The mode of Figure 4 is more suitable for frequency modulation, both be- 
cause the constant frequency curves are closer to each other, and because the 
nearly straight curves are in a region of not too high anode voltages. 
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In order to maintain the operating point on the straight part of the constant 


power curves during the modulation process, it is necessary, as in the previous | 
case, to keep the ratio between the anode-voltage and reflector-voltage varia- _ 


tions constant. 


By choosing the correct point of operation, simply varying the above men- _ 


tioned ratio, it is furthermore possible to pass from amplitude mudulation to 
frequency modulation. In some cases this possibility may be of practical value. 


3. 3. Stability 


In each mode of operation, for lower frequencies, the branches of the 
curves P = const are parallel to the constant frequency curves for considerable 
lengths. This fact is accentuated especially in the modes of Figures 3 and 4. 

This coincidence makes it possible for the tube to operate with high stability 
if anode voltage and reflector voltage are regulated so that at any instant 
their ratio is constant. 

Values of output power which can be as low as 20% of the maximum cor- 
respond to these regions of good stability. However, the power may be high 
enough for laboratory purposes. 


4. Experimental Tests 
4.1. Circuit 


By means of the circuit shown in Figure 5, a number of experiments have 
been conducted. By refering to curves obtained in these experiments, the cor- 
rectness of the previous deductions has been verified. 

The 50 c/s power-line voltage simultaneously modulates anode and reflector 
electrodes of the klystron tube by means of two transformers: 7; and 7, (see 
Figure 5). A variable ratio transformer 7, inserted on the primary 73 allows 
a continuous variation of the ratio between the secondary voltages of 7; and 
T,. A phasing network compensates the phase shift in the transformers. Inas- 
much as in obtaining amplitude modulation, frequency modulation and good 
stability, a decrease of reflector voltage must correspond to an increase of 
anode voltage, the connections of the two transformers are such that the two 
secondary voltages are 180° out of phase. 

The klystron-output power, properly attenuated, is detected by crystal-1, 
and the rectified voltage is applied through a switch to the vertical plates of a 
cathode-ray oscilloscope having a D.C. amplifier. 

By means of a directional coupler, inserted before the attenuator, a portion 
of the power is taken out and applied to a frequency meter. The crystal-2 
voltage, through the above mentioned switch, can in turn be applied, after due 
amplification, to the vertical plates of the oscilloscope. 


ASAT AS ten 
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The Aaa -sweep voltage is derived from the power line through a 
- Phasing network. 


4. 2. Measurements carried out 


4.2.1. In order to obtain amplitude modulation, a curve of frequency a 
_ little higher than the curve Af = 0 of Figure 2, has been utilized in the region 
_ between the power curves P = 50% and P = 100%. 
By adjusting amplitude and phase of the sweep voltage it was possible to 
_ obtain nearly linear amplitude modulation. 
| The maximum frequency deviation was measured with a cavity wave-meter 
_ (Figure 5). It turned out to be about - 1-3 Mc, corresponding to a percentage 
- variation of 1-3 x 10-4. 
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Meter 
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Figure 5 
Circuit for measuring dynamic performance. 


4.2.2. As previously stated, in the case of frequency modulation better 
results are obtained the further the point of operation is moved from the lower 
vertices of the constant power curves, towards the regions where they are 
straighter. It is therefore possible to obtain better performance by moving the 
point of operation into the region above 300 V. The circuit shown in Figure 5 
was used to investigate how the tube works in the region above 300 V. 

By means of the transformer inserted on the anode circuit, a 50-cycle 
voltage was superposed to the 300 V anode voltage. The 50-cycle voltage had 


ZAMP V/17 
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a peak value variable up to a maximum of 150 V. For various values of the 
reflector voltage the crystal-1 response was observed at the oscilloscope. Thus 
it was possible to ascertain the behaviour of the tube in a broad region which 
had not been explored with static measurements not to damage the tube. 

This investigation has verified the expected regularity of the curves even 
above the 300 V ordinate. 

By utilizing the mode of Figure 4, an anode voltage of 350 V and the - 
left branch of the curve P = 45% (not included in the diagram), favourable 
conditions were found to obtain a sufficiently high and constant power together 
with very large frequency swing and fair linearity. 

In this region it was possible to obtain a 40 Mc sweep with a constant output 
power P = 45% (about 13 mW). By accepting power variations of about 10%, 
during the modulation period, and a slight loss in modulation linearity, it was 
possible to increase the bandwidth up to 60 Mc by increasing the modulation 
voltages applied to the two electrodes. 

These large frequency swings can be very useful for measurement purposes, 
for instance to study the response curve of wide-band filters. Of course, with 
smaller frequency swings it is possible to obtain a better linearity by utilizing 
regions where the curve distribution is more uniform. 

It would be possible to further extend the useful sweep width by means of 
a fast acting automatic volume control system. The error voltage obtained from 
the control circuit could be applied to the reflector electrode in such a manner 
that the output amplitude remains almost constant during the modulation 
cycle. 

4.2.3. Finally, the expected operating point with optimum stability was 
observed. To this end, the mode of Figure 4 was used at the point V, = 270, 
V,, = 59. In this point a branch of the curve P = 35% is parallel to the frequency 
curve Af = 15. By means of transformers T; and 7), two sinusoidal voltages 
of peak values 10 V and 2:5 V were superimposed to anode voltage and reflector 
voltage respectively. The change of frequency with respect to the previous 
value was 0:5 Mc while the output power remained practically unchanged. 


4.3. Validity Limits of the Diagrams 


If the static anode voltage is changed to a higher value, the tube dissipates 
more energy. It is accordingly heated up and this, as pointed out at paragraph 2, 
changes the cavity dimensions and hence the frequency. However, if the voltage 
is modulated at audio frequency there is not enough time for the cavity to heat 
up and cool down during the audio cycle. Accordingly there is only an electronic 
frequency change. 

Measurements on the 2 K 22-klystron type show a remarkable thermal fre- 
quency sensitivity with anode voltage variations. The thermal anode sensitivity 
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_is of the same order or still larger than the reflector sensitivity, and the fre- 
_ quency shift is in the opposite direction than for fast voltage variations. 

The behaviour of the 723 A/B tube is much more normal and more predict- 
able from electronic theory, but a small thermal effect will surely be present too. 

Therefore it should be kept in mind that from the static measurement of 
_ the performance curves one cannot exactly predict dynamic behaviour under 
_ audio frequency modulation. For this it would be necessary to measure directly 
dynamic performance. 
| Of course these facts do not invalidate the measurements, because for slow 
_ changes they are perfectly right and for quick modulation they show at least 
the general features of the klystron characteristics. It is expected that the 
_ dynamic curves would be quite similar to the static ones except that their in- 
clination would be changed. 
. It is important to notice that, simply by means of the circuit of Figure 5, 

it is possible to find out, by cut-and-try in various cases, the best conditions of 

operation. Then, the diagram would only give a first indication, and can be 
sketched getting only few points on the curves. The exact condition of operation 
can be found in an easier way by reading the instruments inserted in the supply 
and modulation circuits. 


5. Conclusion 


By means of diagrams drawn out for the 723 A/B Reflex Klystron tube 
conditions of operation were found which proved to be of practical value. 

By simply relating the anode and reflector voltages, so as to keep the ratio 
of their variations constant, it was possible to obtain sinusoidal amplitude 
modulation keeping the frequency practically constant. In a similar way it was 
possible to obtain frequency modulation keeping constant the output power, 
with large frequency swings. Finally, always by means of the diagrams, con- 
ditions of operation characterized by good stability were found. 

The described modulation systems may prove to be of good utility, inas- 
much as they give particular results unobtainable with the signal generators 
available up to date. 
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Zusammenfassung 


Durch eine Reihe von Messungen an einem Reflexklystron wird die Abhangig- 
keit seiner Ausgangsleistung und Frequenz von den angelegten Gleichspannungen 
untersucht und in Diagrammen dargestellt. An Hand dieser Diagramme wird ge- 
zeigt, dass durch gleichzeitige Modulation der Resonatorspannung und der Re- 
flektorspannung entweder reine Amplitudenmodulation oder reine Frequenz- 
modulation des Klystrons erreicht werden kann. Die wichtigste Anwendung 
besteht in der Modulation von Signalgeneratoren fiir Messzwecke. 


(Received: July 14, 1953.) 
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Zur iterativen Auflésung algebraischer Gleichungen 


Von Hans J. MAEHLY, Ziirich?) 


Durch die Entwicklung schneller Rechenautomaten hat das klassische Pro- 
blem der iterativen Berechnung der Nullstellen von Polynomen keineswegs an 
Bedeutung verloren; hingegen haben sich die MaBstabe zur Beurteilung der in 
Frage kommenden Naherungsformeln etwas verschoben: Einerseits k6nnen viel 
kompliziertere Formeln in Kauf genommen werden, andrerseits aber soll die 
Konvergenz in jedem Falle gut sein, auch dann, wenn die Koeffizienten oder der 
Anfangswert vielleicht sehr ungiinstig sind. Bekanntlich erfiillt die Newtonsche 
Formel diese Bedingungen nicht immer, vor allem nicht bei komplexen und bei 
mehrfachen Nullstellen”). Eine weitere Schwierigkeit entsteht oft durch die Ver- 
schleppung von Rundungsfehlern, besonders bei Polynomen hohen Grades, 
worauf wir am Schluss nochmals zuriickkommen. 

Gegeben sei also ein Polynom m-ten Grades: 


TEM ED) :y Cones (1) 
vy=0 
gesucht seine Wurzeln #1, %, ..., ¥,, fiir welche 
n 
Pl) = om [I] (# ~ %) (2) 
a 


Dabei kénnen Koeffizienten und Wurzeln komplex sein. 

Herleitung und Vergleich verschiedener Naherungsformeln zur Berechnung 
der +, werden erleichtert, wenn wir als Hilfsfunktionen die logarithmischen Ablei- 
tungen von P einftihren; wir setzen 


1 d \' ie 1 


eR TL (= Ley Spire ye = Sn a 12 (3) 


1) Institut fiir angewandte Mathematik der ETH. 
*) Vergleiche hierzu zum Beispiel: E. F. Moore, Math. Tables & Aids to Comp. 3, 486 (1949). 


TS a a 


wa © 


a 
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Mit 
Lex! d*® P 
Oh Gage ” 
wird 
S$; = =, , 
hh | 
5 
2S3 = —2 a? + 3a, a,— a, , | ei 
6 Sy = 6 aj — 12 af ag + 3434+ 44, a, — a,. 


Schliesslich fiihren wir — fiir ein fest gewahltes ¥ — zur Abkiirzung ein: 


(4%, — #)* = Uy; (6) 
» somit 
nN 
cn = pa Up, « (7) 
ko 


Wir suchen nun Naherungsformeln, die von jedem x aus méglichst nahe an 
eine Wurzel heranfiihren, und zwar am besten an die nachstgelegene. Wegen (6) 
und (7) bedeutet dies: Aus den ersten s,, soll eines der u;, wenn méglich das absolut 
grosste, mdglichst genau bestimmt werden. 

Ist — von # aus gesehen — eine einzige (p-fache) Nullstelle viel naher als alle 
andern, so gilt bei geeigneter Numerierung: 


Uy=Uzp=- = U,, |U,| > |u,| fiir k>p. (8) 


In erster Naherung ist also: 
: 
Sy= PU, Sg=puy. (9) 


Mit p = 1 ergibt sich die Newtonsche Formel: 


C55 ee Ue Sa =o ; (10) 
allgemeiner bei vorgegebenem Pp: 
ay 
Uy, = — 7 +) 11 
und schliesslich bei unbekanntem Pp: 
Beesies ie 
Uy ra es (12) 


Genauere Formeln, welche allerdings die Auflésung einer quadratischen Glei- 
chung erfordern, erhalt man, wenn man sich die nachstgelegenen Wurzeln nicht 
in einem, sondern in zwei verschiedenen Punkten konzentriert denkt. Ist % eine 
p-fache, ¥ eine q-fache Wurzel, so wird 


Sy = Pum --qo™ mit w= (*#—-*)-1; v= (%¥—-2z)*. (13) 
Sind dabei # und q gegeben, so ergibt sich leicht?) : 
ree Psy Vp q (Pp + q) Sy. — Sq] . (14) 
P(p+ 9) 


1) Diese beiden Formeln finden sich schon bei E. ScuroprER, Math. Ann. 2, 317 (1870). 
2) Vgl. E. Boprwic, Proc. Acad. Sci. Amsterdam 49, 911 (1946), und J. G. VAN DER CorPutT 
Proc. Acad. Sci. Amsterdam 49, 922 (1946). 
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Mit p = 1,qg = n— 1 erhadlt man die Laguerresche Formel!) : 


eo ye ee : (15) 


n 


Ist nur die Summe s, = ~ + ¢ bekannt, nicht aber p und q einzeln, so ergibt sich 
fiir wu die quadratische Gleichung: 


U® (Sp Sp — 8?) — U (Sp Sg — Sy Se) + (Sy Sg — 53) = 0. (16) 
Lasst man auch p + q frei, so erhalt man ganz analog: 
u® (Ss; Sy — $3) — U ($1 Sq — Sq Sg) + (Sp Sq — 53) = 0. (17) 


Die letzten vier Gleichungen fiir « haben jeweils zwei Lésungen; da wir jeweils 
die nachstgelegene Wurzel suchen, ist das Wurzelvorzeichen stets so zu wahlen, 
dass man das absolut grésserve u erhdlt. 

Aus unserer obigen Herleitung ergibt sich, dass die Laguerresche Formel auch 
zur Berechnung komplexer Nulistellen geeignet ist. Bei Polynomen mit reellen 
Koeffizienten hat sie den grossen Vorzug, die reelle Achse von selbst, aber erst 
dann zu verlassen, wenn der Einfluss eines komplexen Wurzelpaars iiberwiegt. 
Die Folge sukzessiver Naherungen konvergiert bei Verwendung der Laguerre- 
schen Formel (auch im Komplexen) kubisch gegen einfache, aber nur linear 
gegen mehrfache Wurzeln?). Analog fiihrt Gleichung (14) nur fiir genau p-fache 
Wurzeln zu kubischer Konvergenz. Dagegen ist die Konvergenz mit (16) stets ku- 
bisch, mit (17) sogay immer von vierter Ordnung, ohne dass die eventuelle Vielfachheit 
dey Wurzel bekannt zu sein braucht. Nach unserer Herleitung ist zu vermuten, 
dass (17) noch etwas bessere Konvergenz ergibt als (16), und in den wenigen 
bisher gerechneten Beispielen hat sich dies bestatigt. Ob und wann sich aber der 
erhebliche Mehraufwand ftir die Berechnung von s, lohnt, kann hier noch nicht 
entschieden werden. 

Nach Bestimmung einer Nullstelle +, pflegt man das Polynom P(*) durch 
(~ — *,) zu dividieren, einerseits um den Grad des Polynoms zu erniedrigen, also 
die folgende Rechnung zu vereinfachen, andrerseits um eine nochmalige Berech- 
nung von #; zu verhindern. Dies fiihrt aber oft zu gefahrlichen Rundungsfehlern. 
Beischnellen Maschinen, wo die Ersparnis an Rechenaufwand weniger ins Gewicht 
fallt, ist daher eher folgendes Verfahren zu empfehlen: 


Kennt man Naherungen %, ..., #, fiir / Nullstellen #,, ..., x, und sucht man 
die restlichen Wurzeln #,,,, ..., ¥,, S80 ersetze man alle s,, durch 
l 
, F- _ 
Sm = Sm — » (4, — *)-™. (18) 
ke] 


Im Gegensatz zum tiblichen Divisionsverfahren haben dann die Ungenauigheiten 
4X, — X_% praktisch keinen Einfluss auf die iterierten Werte und gar keinen auf deren 
Grenzwerte, solange | * — x; | > |¥_— *%,|. Ist diese letztere Bedingung nicht 
erfiillt, so kann die Folge allerdings gegen ein *, oder eventuell gegen ein 4, mit 
k < / konvergieren. 

Theoretische Untersuchungen iiber die hier vorgeschlagenen Iterationsfor- 
meln sowie numerische Versuche sind im Gange, und eine ausfiihrliche Arbeit 
hiertiber ist vorgesehen. 


1) E.N, Lacuerre, Nouv. Ann. Math. [2¢s.1 19 (1880), oder Giuvres de Laguerre, Bd. I, S. 87 
bis 108. 


*) Vgl. E. Bopewic und J. G. VAN DER Corput, a. a. O, 


~ 
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Summary 


The slow convergence of Newton’s method for the iterative approximation 
of the roots of an algebraic equation may be improved by using a formula due to 
LAGUERRE. It is shown that this formula together with certain generalizations 
are useful for the numerical calculation of both real and complex zeros. A method 
is also given for the calculation of a further root of an equation after several 
roots have already been found, such that errors on the determination of the 
previous roots will not affect the accuracy with which the present one may be 
approximated. 


(Eingegangen: 21. Januar 1954.) 


Aerodynamische Krafte auf schlanke Rotationskérper') 


Von Hans R. VoELiMy, Ziirich?) 


Zur Berechnung der Krafte auf schlanke drehsymmetrische K6rper, insbeson- 
dere zur Bestimmung der Querkraft, steht heute eine recht umfangreiche Litera- 
tur zur Verfiigung. Seltener aber findet man Ver6ffentlichungen iiber Messergeb- 
nisse und Vergleiche mit theoretischen Berechnungen. 

Im Verlaufe vieler routinemassig durchgefiihrter Windkanalmessungen wurde 
der Verfasser auf betrachtliche Differenzen der erhaltenen Messresultate mit den 


_vorliegenden potentialtheoretischen Berechnungen aufmerksam. Eine spezielle 


_ Versuchsserie, welche im Institut fiir Aerodynamik der ETH. unter Leitung von 


Herrn Prof. Dr. J. AcKrReEr durchgefiihrt wurde, sollte den Anwendungsbereich 
besonders der einfachen Theorien aufzeigen. 


1. Elementare Theorie 


Die Normalkraft N, welche auf einen zur Stromungsrichtung mit dem kleinen 
Winkel « angestellten schlanken Rotationsk6érper wirkt, wurde seinerzeit schon 
von Munk /[1]’) angegeben und seither von JONES[2] und ACKERET [3] wieder 
aufgenommen. Die hier getrennte Betrachtung von benachbarten Schnittebenen 
senkrecht zur Rotationsachse stellt eine noch starkere Vereinfachung dcr tat- 


' sachlichen Stromungsvorgange dar, als dies mittels der Linearisierung der rdum- 


lichen Potentialstr6mung nach den Methoden von von KARMAN und Moore [4], 
FERRARI [5] und TsIen [6] méglich ist. Dabei zeigt sich, dass nach dieser Quer- 


_ schnittstheorie die Normalkraftverteilung unabhangig von der Mach-Zahl wird und 


nur von der Anderung des Querschnittes F mit der Lange ¥ des Kérpers abhangt: 
aN dF 


Se Te 


Wird iiber die Lange des K6rpers integriert, so erhalt man das erstaunlich ein- 
fache Resultat, dass die Normalkraft zur Bodenflache F, proportional ist: 


N=sin2a-£V?F,. 


1) Vorgetragen am Astronautischen Kongress 1953 in Ziirich, zufolge Abwesenheit des Verfassers 
durch Herrn B. Cuaix, Institut fiir Aerodynamik der ETH. 

2) Werkzeugmaschinenfabrik Oerlikon, Biihrle & Cie., und Institut fiir Aerodynamik der ETH. 

3) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 269. 


264 Kurze Mitteilungen — Brief Reports - Communications bréves PS ZAMP 
Ebenso ergibt sich fiir das Lingsmoment H eine physikalisch anschauliche Grosse : 
M =sin2«- 3 V2(G—F, 5), 


indem G das Volumen des Korpers und s den Abstand des Momentenbezugs- 


punktes von der Bodenflache bedeuten. 
Die Versuche haben die Giiltigkeit dieser einfachen Formeln fiir kleine A nstell- 
winkel und Koérper mit schwach konvergierendem Heckteil bewiesen. 


2. Modelle und Messeinrichtung 


Um den Einfluss der Grésse der Bodenflache auf die Normalkraft zu unter- 
suchen, wurden fiinf Versuchskérper im hinteren Drittel der Lange verschieden 
stark verjiingt (Figur 1). 

Heck= 


tangenten <= oe 


ee 
aes 2,25 
> 450 
_——— 6.75 


Korper Nn 


Cy) 


Se 
eS 


3 SS J 0,396 
— 9,00 0,240 
[ — 
Figur 1 


Meridianschnitt der fiinf Versuchsk6rper mit verschieden grossen Bodenflachen. Die elektrische 
Drei- Komponenten-Waage sitzt im Innern des Modells. 


Die Messung der Krafte geschieht mittels einer tm Modell eingebauten elek- 
tvischen Waage (straingage), welche auf einer Stiitze von veranderlichem Anstell- 
winkel sitzt. An einer Philips-Messbriicke werden die Resultate abgelesen und 
in einem von der Firma Contraves AG., Ziirich, entwickelten elektronischen 
Analogierechner [7] ausgewertet. 


3. Messergebnisse bei grossen Anstellwinkeln 


Wie Figur 2 deutlich zeigt, tibersteigt die gemessene Normalkraft stark die 
aus der elementaren Theorie (Kurve a) oder auch aus der linearisierten Potential- 
theorie (Kurve 6) berechneten Werte bei grossen Anstellwinkeln. Das quadra- 
tische Verhalten dey Normalkraft iiber dem Anstellwinkel kann nur durch die Be- 
trachtung der Winfliisse der Zahigkeit der Luft auf die Stromungsvorgange erklart 
werden, worauf Brrz [8] erstmals hingewiesen hat. Ahnlich der Queranstrémung 
eines Zylinders entsteht durch die Queranstrémung des Rotationskérpers auf der 
Abstromseite, im Querschnitt betrachtet, ein Wirbelpaar. Die gleichzeitige Langs- 
anstromung verhindert aber die Ausbildung einer Karmanschen Wirbelstrasse, 
indem die Wirbel in der dritten Dimension wegschwimmen kénnen. In den beiden 


ae 
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; 
a: 


be 


Korper Nr. 


Figur 2 
Vergleich der gemessenen Normalkraftbeiwerte Cy fiir die fiinf KGrper bei M = 1,38, Rez = 1,2- 108 
mit: a Elementartheorie von MunxK (Q-Theorie); b linearisierter Theorie; c Elementartheorie und 
nichtpotentieller Normalkraft. 


cant 
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Figur 3 


Schlierenaufnahme der Wirbelzopfe auf der Abstrémseite des K6rpers (Pfeile), eine Folge der durch 
die Potentialtheorie nicht erfassten Zahigkeit der Luft. Kérper Nr.3, M = 1,66, Rey = 1,2- 10% 
= 10°. a 10°-Anstellwinkel in der Bildebene; 6 10°-Anstellwinkel senkrecht zur Bildebene. 


in zwei Projektionen aufgenommenen Schlierenbildern zeichnen sich diese Wir- 
belzopfe als starke Helligkeitsveranderungen deutlich ab (Figur 3). 

Die durch die Wirbel hervorgerufene Querkraft pro Einheitslange in Achs- 
richtung ist: 


ipo ae (V sina)? D~a?. 

Mit dem konstanten Widerstandskoeffizienten C, = 1,2 (unterkritisch angebla- 
sener Kreiszylinder) berechnet, ergibt sich der in Figur 4, Kurve c, abgebildete, 
zum Durchmesser D jedes Querschnittes proportionale Verlauf der Normalkraft. 
Die resultierende Kraft wird also durch den Schwerpunkt der Meridianflache 
gehen. Da dieser Schwerpunkt im allgemeinen nahe beim Massenschwerpunkt 
des Ko6rpers liegt, wird das Drehmoment, welches die Stabilitat des K6rpers beein- 
flusst, durch die nichtpotentiellen Krafte haum verdndert. 


aco 


~ 9,002 


Figur 4 
Normalkraftverteilung entlang der Kérperachse fiir Kérper Nr. 2 bei sina = 0,2, nach: a Quer- 
schnittstheorie; b linearisierter Theorie, M = 1,38; c nichtpotentieller Normalkraft aus Quer- 


anstromwiderstand infolge Wirbelbildung, Cw, quer = 1,2. 


3 


“4 


' 
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Aus diesem Grund ist die Abschatzung der Langsstabilitat eines schlanken 

Geschosses mittels der linearisierten Potentialtheorie (Kurve b) oder sogar mit- 
_ tels der einfachen Querschnittstheorie (Kurve a) erstaunlich genau. 

Der Vergleich der von den Kurven eingeschlossenen Flachen hingegen zeigt 
den starken Einfluss der durch die Zahigkeit hervorgerufenen Normalkrdifte, 
welche die von der Potentialtheorie erfassten Krafte libersteigen k6nnen. 

Nach ALLEN [9], VAN Dyke [14] und Frax und Lawrence [11] addiert man 
die nichtpotentielle zur potentiellen Normalkraft und erhalt so die Kurven c in 

' Figur 2, welche nun das quadratische Verhalten der experimentellen Kurven 
besser wiedergeben. 


4. Messergebnisse bei kleinen Anstellwinkeln 


Bei kleinen Anstellwinkeln méchte man erwarten, dass die potentialtheore- 
tisch berechneten Normalkrafte gut mit den Messungen iibereinstimmen, da das 
in « quadratische Glied seine Wirkung verliert. Schon ein Blick auf Figur 2 
belehrt aber, dass besonders bei den sich stark verjiingenden Kérpern die Tangen- 
tenneigungen der Normalkraftkurven im Anstellwinkel-Nullpunkt wenig mit den 


Korper Nr = 4 


Messgenauigkeit\ 


+02 


Feoa ene 


0,24 0,394 0,56 0,766 ] Fax 
; Figur 5 
Vergleich der gemessenen Querkraftgradienten dC y/de bei « = 0° mit: a O-Theorie; 6 linearisierter 
Theorie. 


Neigungen der theoretischen Kurven tibereinstimmen. In Figur 5 sind diese Nor- 
malkraftgradienten iiber der Bodenflache aufgetragen, um die Abweichung der 
Messungen gegentiber der durch eine Gerade dargestellten Querschnittstheorie 
(Kurve a) oder der linearisierten Theorie (Kurve b) besser tiberblicken zu konnen, 

Die Verdrangungsdicke der Grenzschicht wird gegen das K6rperende hin 
um so grésser, je starker sich der Kérper nach hinten verjiingt. Dies zeigt Figur 6 
deutlich, wo die Grenzschicht als weisser Streifen auf der oberén Seite der Kérper 
deutlich sichtbar ist. Die Berechnung der Druckverteilung entlang der verschie- 
denen Kérper zeigt das starke Anwachsen des Druckes gegen das KG6rperende hin, 
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Figur 6 
Wachsende Grenzschichtdicke bei zunehmender Konvergenz des Heckteils. Schlierenaufnahmen 
. ‘ F; Pee 
der K6rper Nr. 0, 2 und 4 bei M = 1,38, Rez = 1,2- 108, a = 0°. 


3 2 1 


Figur 7 


Korper Nr —= 4 


Vergleich der gemessenen Querkraftgradienten mit den bis zum Punkt minimalen Druckes inte- 
grierten Querkraften nach: a Q-Theorie; b linearisierter Theorie. 


. 
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_ woraus mit der rein geometrischen Verkleinerung des Kérperumfanges zusammen 


das starke Anwachsen der Grenzschicht hervorgeht. 
In ganz schematischer Weise wird nun angenommen, dass die Strémung 


_ausserhalb der Grenzschicht vom Querschnitt des tiefsten Druckes an infolge der 
_ stark anwachsenden Grenzschichtdicke nicht weiter gegen die Achse konvergiere. 


An Stelle der Bodenflache werde also der Querschnitt minimalen Druckes in die 


Rechnung gesetzt. Figur 7 zeigt, dass nach dieser einfachen Hypothese die experi- 
mentellen Werte bedeutend besser erreicht werden. 
Zusammenfassend kann gesagt werden, dass es wichtiger ist, die Einfliisse der 


_Ziahigkett der Luft bei der Berechnung der Krafte auf schlanke Rotationskérper 
_abzuschatzen, als die potentialtheoretischen Rechenmethoden weiter zu ver- 


1 al 


feinern. 
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Summary 


Comparison between measurements of the normal force of five slender bodies 
of revolution with different base areas in supersonic flow and the simple airship- 


hull theory of Munk or the more accurate linearised theory of potential flow show 


considerable discrepancy not only at large, but also at small angles of attack. 

To get better agreement with the experimental results, it is shown to be more 
important to estimate the effects of viscosity than to consider higher orders of 
potential flow theory. 


_ (Eingegangen: 5. Februar 1954.) 


270 Kurze Mitteilungen — Brief Reports - Communications bréves Pa ZAMP 


Abklingen eines schwachen Dralls in der Langsstromung zwischen 
koaxialen Kreiszylindern 


Von HERMANN WUNDT, Freiburg i. Br.*) 


L. Cottatz und H. GOrTLER haben neuerdings?) die laminare Strémung mit 
schwachem Drall durch ein ruhendes gerades Kreisrohr rechnerisch behandelt. 
Nach Linearisierung der Bewegungsgleichungen beziiglich der Drallst6rungs- 
grésse wird man auf ein Eigenwertproblem einer gew6hnlichen Differentialglei- 
chung zweiter Ordnung gefiihrt. Die Gréssenordnung der Eigenwerte gibt dabei 
Aufschluss tiber das Abklingen des anfanglichen Dralls in Stromrichtung. 

Diese Betrachtungen sollen hier auf eine Langsstré6mung zwischen koaxialen 
Kreiszylindern verallgemeinert werden. Den Drall mag man sich dabei etwa 
durch ein im ruhenden Rohr angebrachtes geeignetes Leitwerk erzeugt denken; 
der innere Zylinder kann als feste Nabe ausgebildet sein. In den technisch inter- 
essierenden Fallen wird diese Strémung freilich selten laminar sein. Die das Ab- 
klingen des Dralls messenden Eigenwerte werden von der Weite des ringf 6rmigen 
Spalts abhangig sein; man erwartet fiir enge Spalte ein rascheres Abklingen, das 
heisst absolut gréssere Eigenwerte. 

Wegen der Herleitung des Eigenwertproblems sei auf die genannte Arbeit von 
CoLLtatz und GORTLER verwiesen (fortan zitiert als [1]). Mit den dortigen Be- 
zeichnungen seien g = #/7, Re*) die axiale, s=y/ry) die radiale dimensionslose 
Koordinate des der Str6mung angepassten Zylinderkoordinatensystems. s = } 
sei der Radius des inneren Zylinders, s = 1 der Radius des ausseren Zylinders, 
so dass 6 das Radienverhaltnis bedeutet. Betrachtet werden Falle von b= 0 
(kein innerer Zylinder) bis b = 1 (Rohr verstopft). Das s-Intervall ist b= s <= 1. 

u(s) sei die Komponente der Stromungsgeschwindigkeit in g-Richtung, wah- 
rend w(q, s) die tangentielle Geschwindigkeitskomponente, also den Drall, be- 
zeichne. Fiir diesen Drall wird der Separationsansatz 


w(q, s) = e-49 G(s) 


gemacht, welcher auf das untenstehende Eigenwertproblem (1) fiir die Funktion 
G(s) fiihrt. Der Separationsparameter A tritt als Eigenwertparameter in die 
Gleichung fiir G(s) ein. Er ist gegeniiber [1] mit umgekehrtem Vorzeichen defi- 
niert, da dann positive Eigenwerte zu erwarten sind. 

Damit lautet das gleich in selbstadjungierter Form geschriebene Eigenwert- 
problem: 


L{G] = -(sG' + — G = Asis) G (1) 


mit den aus der Haftbedingung am inneren und am Ausseren Zylinder herriih- 
renden Randwerten G(b) = G(1) = 0. 

Der wesentliche Unterschied zum in [1] aufgestellten Eigenwertproblem be- 
steht in der verénderten Funktion f(s), die die axiale Geschwindigkeitsverteilung 


1) Mathematisches Institut der Universitat, Abteilung Angewandte Mathematik. 
*) L. Corzatz und H. Gorter, Rohrstromung mit schwachem Drall, ZAMP 5, 95 (1954). 
3) Uber die Definition der hier Sreaieten Reynoldsschen Zahl Re siehe unten wach Formel (3). 
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_ in weitem Abstand vom Einlauf darstellt. Ohne inneren Zylinder (b = 0) war diese 


u*(s) 
Uo 


= al() = Se 


die bekannte Poiseuille-Strémung (U, die Maximalgeschwindigkeit in der Rohr- 

achse), jetzt hingegen, nach Hinzutreten der weiteren Randbedingung u(b) = 0, 
lautet sie (s) 
u(s 

Sy ree a Sa LS (2) 


mit der vom Radienverhaltnis b abhangigen Konstanten 


1— }? 
AD \ Sa at (3) 
_ Wir wollen hierbei fiir alle Spaltweiten b stets dasselbe konstante Druckgefalle 
— dp/dx zugrunde legen. Damit lasst sich eine von b unabhdngige Bezugsge- 
schwindigkeit 
sil i dpe rt 
. dx 44 


(w ist die konstante Zahigkeit) definieren; sie ist mit der bei [1] im Falle b = 0 
verwendeten Bezugsgeschwindigkeit identisch. Die so zu bildende Reynoldssche 
Zahl 

Wott U%e _ 4 v5 @ 


be dx 4p? 


ist gleichfalls von b unabhangig. 

f(s) ist eine exakte Lésung der vollen Navier-Stokesschen Gleichungen. Die 
zu berechnenden Eigenwerte 4; werden nun natiirlich von 6 abhangen: A; = A,(d), 
und gerade diese Abhangigkeit (etwa bei unverandertem Druckgefalle) besitzt 
als Verallgemeinerung gegenitiber [1] das Hauptinteresse. 

Die Theorie der Eigenwertaufgabe (Existenz von Eigenwerten und Eigen- 
funktionen, Entwicklungssatz) ist fast wortlich aus [1] tibertragbar. Erhdhte 
Vorsicht fordert lediglich der Umstand, dass f(s) nunmehr an beiden Randern 
verschwindet. An die nach Eigenfunktionen G; entwickelbaren Funktionen h(s) 
muss daher ausser den in [1] genannten Bedingungen noch die Forderung 
L{h(s)],-» = 0 gestellt werden (fiir s = 0 war dies von selbst sichergestellt), damit 


pe(g) — LOU 
Vs f(s) 


im Intervall (b,1) quadratisch integrabel ist. Alle iibrigen Aussagen bleiben er- 
halten; insbesondere sind die Eigenwerte A; samtlich positiv. 

Wahrend demnach die Theorie der erweiterten Aufgabe nur eine geringfiigige 
Verscharfung der Voraussetzungen notwendig macht, stellen sich der numerischen 
Berechnung der Eigenwerte und -funktionen betrachtliche Hindernisse in den Weg. 

Um iiber die Abhangigkeit der 2; von b und tiber ihre Gréssenordnung Anhalts- 
punkte zu gewinnen, wird zundchst ein Vergleichsproblem herangezogen. Wir 
ersetzen f(s) durch g(s) = 1. Dann gilt g(s) => f(s) iiberall im Intervall (}, 1). Das 
Vergleichsproblem lautet in selbstadjungierter Form 
~G=xsG, G(b) =G(1) =0. (4) 


(sy + 
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Die allgemeine Lésung hiervon ist Z (Vx s), eine Zylinderfunktion erster Ord- 


nung, die den Forderungen Z,(b x) = Z,(//) = 0 zu geniigen hat. Da fiir die 
Eigenfunktionen ein Faktor frei bleiben kann, so lautet diese Forderung: 


TV*) + 4) NV) =0 und J,0)x) + 4,N,0)%) = 0. (5) 


Zwecks Eliminierung von A, multiplizieren wir die erste Gleichung mit N,(d x), 
die zweite mit N, (Vx), und subtrahieren sie von der ersten. Dann folgt 


SV*) N,(b Vx) a Si(6 Vx) N, (Vx) =0, (6) 


eine viel untersuchte transzendente Gleichung. Man findet ihre Wurzeln fiir ver- 
schiedene Werte von b, etwa bei JAHNKE-EMDE}). 

Mit diesen Werten ergibt sich die folgende Tabelle der Eigenwerte x; in Ab- 
hangigkeit von 6: 


b | Hy | Hy | Hy | 4 | Hs | 6 

0) 14,682 49,219 103,50 ViT,52 271,28 384,78 
1/2 40,883 159,386 359,27 633,15 988,46 1422,72 
2/3 89,927 356,424 800,56 1422,35 2221,79 3198,88 
5/6 356,198 1422,127 3198,65 5 685,79 8883,55 12791,90 

ib Fiir b> 1 wachsen die x, iiber alle Grenzen. 


Wegen g(s) 2 f(s) gilt nun x, S /;. Wir haben damit (grobe) untere Schranken 
fiir die 1,(b). Dass sie sehr grob sind, sieht man schon aus dem Vergleich der 
Zeile b = 0 mit den bei [1] berechneten Werten 4,, A, As. Aus der Theorie der 
Gleichung (6) folgt, dass die x,(b) in 0 < 6 = 1 monoton von einem Anfangswert 
nach + oo steigende Funktionen sind. Dasselbe gilt auch fiir die /,(b), wie in Fuss- 
note 1, S. 274, gezeigt werden wird. Das heisst: je enger der Spalt zwischen den 
koaxialen Zylindern ist, desto rascher klingt der Drall stromabwéarts ab; fiir 
b = 1 (Spalt zu) kann natiirlich nichts mehr fliessen. 

Zur eigentlichen Berechnung der Eigenwerte bieten sich im wesentlichen vier 
Wege: 

1. Differenzenverfahren mit Versuchseigenwerten, die so lange systematisch vari- 
iert werden, bis der zweite Randwert mit gewiinschter Genauigkeit erreicht 
ist; Eigenfunktion und Eigenwert sind dann richtig getroffen; 

2. numerische Integration (etwa Runge-Kutta-Methode) mit demselben Probier- 
verfahren fiir die Eigenwerte; 

3. Reihenentwicklungen um den einen Randpunkt; 

4, Verfahren der schrittweisen Naherungen. 

Die Methoden 1 und 2, die man im allgemeinen nur zur Berechnung der Eigen- 
funktionen, selten zur Bestimmung der Eigenwerte anwendet, scheiden aus, da 
das Problem noch den Parameter b enthalt. Man miisste fiir jeden Parameterwert 
b alle Rechnungen gesondert durchfiihren, was einen unangemessenen Arbeits- 
aufwand bedeuten wiirde. 

Methode 3, die sich fiir den Fall b = 0 sehr bewahrt, versagt auch, weil die 
Funktion f(s) nun einen Logarithmus enthalt. Die iibersichtliche Rekursions- 


1) E, JANKE und F, Empr, Tafeln hoherer Funktionen, 4. Aufl. (Teubner, Leipzig 1948), S. 204. 
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formel fiir die Koeffizienten der Potenzreihe 


ie) 
G(s) = pas a, svt 
v=0 


_ geht verloren, weil man mit der Logarithmusreihe auszumultiplizieren hat. Zu- 
dem ist der giinstige Entwicklungspunkt s = 0 nach ausserhalb des Intervalls 
_ geriickt. 

Es bleibt daher nur Methode 4, das Verfahren der schrittweisen Naherungen, 

welches bekanntlich leider nur den ersten Eigenwert mit ertraglichem Aufwand 
_ liefert. Gerade der kleinste Eigenwert A, ist jedoch in erster Linie interessant, da 

er hinreichend weit stromabwarts allein noch die Drallverteilung bestimmt. 

(Dennoch kann man sich iiber die Gréssenordnung des zweiten Eigenwertes ein 
_ Bild machen; siehe unten). 

Zum Vergleich mit [1] wird der Rechengang grob skizziert: Um nicht fiir 
jede Wahl von 6} getrennte Rechnungen durchfiihren zu miissen, wird von einer 
stetigen Funktion Fy) ausgegangen, die die Randbedingungen erfiillt und den 
Parameter 6 explizit enthalt. Dieser Parameter wird durch alle Operationen ohne 

- Spezialisierung mitgeschleppt. Da der Logarithmus im Verlauf der Rechnung 
ohnedies auftritt [wegen f(s)], bietet sich natiirlich f(s) selbst als erste Vergleichs- 
funktion F, an; es erfiillt beide Randbedingungen. 

| Bezeichnet N[F] die rechte Seite des Eigenwertproblems (1), so berechnet man 
weitere Vergleichsfunktionen F,, F,,... aus den Randwertproblemen 


LF 431) = NB). 
Man erhalt so 


1 ee s6 + a st — = s?+ A? 5? (- 3 In? s + = In s — 3) | (7) 
+A [s* (- ‘s In s + ss) 1 52 (Gms- >) ete ees 

mit 

ae S & At (1 — 63) + 5 A(t BH) + = a—o], = 

ee [op At —d) + Be A (Lb) + Ge (1 —09)] 


Die Formel fiir die Funktion F, wiirde bereits eine sehr betrachtliche Lange an- 


nehmen. 
Daraus werden die Schwarzschen Konstanten 


1 1 
ig 3 }(s) F3(s) ds; a, = [s f(s) Fo(s) Fy(s) ds; ae =/s f(s) Fi(s) ds (9) 
b 5 b 


und die Schwarzschen Quotienten m, = @)/a, und tg = @,/a, berechnet. Sowohl 
die Schwarzschen Konstanten als auch die Quotienten enthalten den Parameter b, 
der im Integranden und in der unteren Integralgrenze steckt. Nun schrumpft 
fiir b > 1 nicht nur die Intervallange auf 0 zusammen, sondern auch der Betrag 


ZAMP V/18 
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des Integranden wird rasch klein. Daher gehen die a; fiir b > 1 sehr stark gegen 0. 
Die y;, deren Folge von oben gegen den ersten Eigenwert 4, konvergiert, nahern 
sich dadurch fiir wachsendes b immer mehr der Form 0/0, worunter die erzielbare 
Rechengenauigkeit stark leidet. 

Um mit Hilfe von y,(b) und ,(b) schliesslich 4,(b) in Grenzen einschliessen zu 
k6énnen, benétigt man noch eine untere Schranke /,(b) von /,(b). Unter der durch 
Uberschlagsrechnungen kontrollierbaren Annahme, dass /,(b) fiir wachsendes 6 
jedenfalls nicht schwacher wachst als ,(b), bildet 


150) = 215) Aal0) 


eine verwendbare, wenn auch vielleicht noch sehr grobe untere Schranke fiir A . 
Die Schranken aus dem Vergleichsproblem von oben geniigen nicht, da sie die 
Forderung /, > mw, nicht erfiillen, das heisst viel zu schlecht sind. Die sogenannte 
Hauptformel 
Sahel ee hii one ig ages 10 
i) (12/42) — 1 SAS bs ( ) 
liefert sodann die in der Tabelle wiedergegebenen Werte fiir A,(b)+). 
Zu Vergleichszwecken werden hier auch noch folgende Gréssen mit angegeben: 
1. Die Strecke, auf der der Drall zu einem vorgebbaren Bruchteil seines anfang- 
lichen Betrages abgeklungen ist, im Verhaltnis zu der Strecke, auf der das 
Abklingen zum selben Bruchteil ohne inneren Zylinder erfolgt. Dieses Ver- 
haltnis d(b) veranschaulicht also, um wieviel schneller das Abklingen des 
Dralls bei abnehmender Spaltweite erfolgt?). 


4,(0) 
A,(6) © 


2. Die Durchflussmenge Q(b) im Verhaltnis zur Durchflussmenge Q(0) bei Fehlen 
des inneren Zylinders. 


Q(b) = $75 Uy (1— b+ Cm): 


d(b) = 


In b 


3. Die noch erreichte Maximalgeschwindigkeit u,,,,(b) im Verhaltnis zur Maxi- 
malgeschwindigkeit w,,¢_(0) = Up, ohne inneren Zylinder. 


Umac() _ ee <7 
ira cea c a In <acege 


4. Die Abszisse sy(b), auf der die Maximalgeschwindigkeit w,,,,(b) angenommen 
wird. 


sy(b) = a 


5. Die Abszisse sp(b), auf der das Maximum der asymptotischen Drallverteilung 


G, max(b) liegt [der graphischen Darstellung der ersten Eigenfunktion G,(s) 
fiir die verschiedenen 6 entnommen]. 


1) Die 4; sind im offenen Intervall <0, 1> stetige, stetig differenzierbare und — wie sich nach- 
rechnen asst — monoton wachsende Funktionen von b. Ihre Folge konvergiert daher in <0, > 
auch gleichmassig gegen A,(b), das demnach gleichfalls monoton wachst. 


*) Die Werte gelten nur asymptotisch fir grosse g, da der Einfluss der héheren Eigenwerte 
nicht beriicksichtigt ist. 


to 
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Es sei betont, dass alle diese Werte auf der Voraussetzung gleichbleibenden 
Druckgefalles, das heisst von konstantem U, bzw. Re fussen. 


21,382 302 

39,6 0,463 64 

60,5 0,54611 
0,614.74 


0,677 03 
0,735.53 
0,79147 
0,84555 
0,89815 
0,949 56 
1 


Alle Daten sind als um etwa eine Einheit der letzten Ziffer unsicher anzusehen. 
Der Ort des asymptotischen Drallmaximums stimmt fiir nicht allzu schlanke 
innere Zylinder praktisch mit dem Ort des Geschwindigkeitsmaximums tiberein. 

Das Maximum der asymptotischen Drallverteilung (wie auch das der maxi- 
malen Durchflussgeschwindigkeit) liegt fiir alle b ndher am inneren als am ausse- 
ren Zylinder, wie aus der Asymmetrie der ersten Eigenfunktion zu ersehen ist. 


Figur 1 
Die erste Eigenfunktion G,(s) fiir verschiedene Spaltweiten 6 [Normierung derart, dass Gy(1) fiir 
alle Kurven gleich ist.] 
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Summary 


The investigations of Corratz and G6RTLER (ZAMP 5, 95 [1954]), are gener- 
alized to a laminar longitudinal flow between coaxial circular cylinders, super- 
posed by a small angular velocity component. It is found that the eigenvalues A 
describing the damping of the angular momentum downstream are rapidly ~ 
increasing with increasing radius b of the inner cylinder. The first eigenvalue 
and eigenfunction are calculated for several values of Db. 


(Eingegangen: 23. November 1953.) 
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Historische Anmerkung zum Aufsatz: Einige Anwendungen funktional- 
analytischer Methoden in der praktischen Analysis‘) 
Von LotHar CoLtiatz, Hamburg 


Einem Vortrage von Herrn A. OSTROWSKI (am 14. Oktober 1953 in Oberwolf- 
ach) und einem Gespraéch mit Herrn OsTRowskI verdanke ich den Hinweis auf 
einen sehr beachtlichen Briefwechsel von R. MEHMKE und R. A. NEKRASSOFF, der 
in Math. Sbornik (Moscow) 76, 437-459 (1892) (drei Briefe von NEKRASSOFF in 
Russisch, zwei Briefe von MEHMKE in Deutsch) ver6ffentlicht ist und in welchem 
bereits verschiedene Konvergenzkriterien fiir das Iterationsverfahren in Einzel- 
schritten bei linearen Gleichungssystemen bewiesen sind. Insbesondere ist das fiir 
praktische Zwecke sehr ntitzliche Kriterium AK <1 mit dem nach Gleichung (28) 
meiner Arbeit berechneten A nicht erst von SASSENFELD eingefiihrt worden, wie 
ich bisher annahm, sondern findet sich auch schon in dem genannten Brief- 
wechsel S. 439, und auf den folgenden Seiten gibt MEHMKE einen Beweis. 

Zu Abschnitt 9 ist noch als Literatur hinzuzufiigen: 

A. Ostrowski, Uber die Determinanten mit tiberwiegender Hauptdiagonale, Com- 
ment. Math. Helv. 70, 69-96 (1937). 


2. Internationale Tagung 
fur Kurzzeitphotographie und Hochfrequenzkinematographie 


Die 2. Internationale Tagung fiir Kurzzeitphotographie und Hochfrequenz- 
kinematographie wird vom 22. bis 28. September 1954 in Paris stattfinden. Das 
Organisationskomitee besteht aus den Herren LuciEN Butt, Direktor des Marey- 
Instituts (Vorsitzender), Ingénieur Général P. Lrpessart, Prof. Dr. H. ScHARDIN, 
Ingénieur en Chef P. Fayorie, J. Vivié, Sekretar der AFITEC (Kassenfiihrer), 
Ingénieur Principal P. Nasiin (Sekretar) sowie auslindischen Delegierten. Der 
Schrifiverkehy, der die Tagung betrifft, ist zu fiihren mit dem Sekretir des 
Organisationskomitees, M. Ingénieur Principal P. Nastrn, Laboratoire Central de 
l’Armement, Fort de Montrouge, Arcueil (Seine), Tél. Alésia 06.46. 

Vortrage werden gehalten nicht nur iiber die optischen, mechanischen, elek- 
trischen und elektronischen Methoden, Verfahren, Gerdte und Zubehovapparate der 
Kurzzeitphotographie und -radiographie und der Hochfrequenzkinematographie 


1) L. Coxtatz, ZAMP 4, 327-357 (1953). 


egy? 


i 
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_ und -radiokinematographie, sondern auch iiber die Anwendungen dieser Methoden 
_ zur Untersuchung schnellaufender Vorgange in den verschiedensten wissenschaft- 


lichen und technischen Gebieten (Physik, Aero- und Gasdynamik, Werkzeug- 
maschinen, Metallurgie, Elektronik, Motoren, Ballistik, Medizin, Chirurgie, Bio- 
logie usw.). Die Hochfrequenzkinematographie bezieht sich im allgemeinen nicht 
auf Aufnahmeserien, die mit Bildfrequenzen unter 100s aufgenommen worden 


sind. 


Die drei Tagungssprachen sind Franzésisch, Englisch und Deutsch. Der Ta- 


_ gungsbeitrag betragt fiir die Vortragenden 1000 Francs, fiir die iibrigen Teilnehmer 


2000 Francs. Jeder Teilnehmer erhalt vor der Tagungseréffnung Zusammen- 


_ fassungen aller Vortrage. In den Grenzen des zur Verfiigung stehenden Raumes 


konnen die Teilnehmer photographische und kinematographische Apparate und 


_ Zubehor, gegebenenfalls kommerziellen Charakters, ausstellen. Fiir die Aus- 
_ stellung jedes Gerdtes wird ein Beitrag von 20000 Francs erhoben. 


Im Auftrage des Organisationskomitees: J. EGGERT, Ziirich 


Generalversammlung der Union Radio-Scientifique Internationale 
vom 23. August bis 2.September 1954 in den Haag 


Die Union Radio-Scientifique Internationale (U.R.S.I.) wird ihre diesjahrige 
Generalversammlung vom 23. August bis 2. September 1954 in den Haag (Holland) 
abhalten. Nahere Auskiinfte erteilen die Prasidenten der Nationalkomitees oder 
das Generalsekretariat in Briissel (42, rue des Minimes). Das Sekretariat des 
Organisationskomitees legt in den Handen von Herrn Dr. C. TH. VAN DER WIJCK 
(Direction générale des P.T.T. des Pays-Bas, Kortenaerkade 12, ’s Gravenhage). 
Anmeldungen zur Teilnahme sind moglichst friihzeitig erbeten. Die Programme 
werden noch erscheinen. F, TANK 


Internationales Symposium tiber experimentelle Meteorologie in Zurich 
; vom 4. bis 6. Oktober 1954 


Die «Eidgenéssische Kommission zum Studium der Hagelbildung und der 


_ Hagelabwehr» wird an der Eidgendssischen Technischen Hochschule in Ziirich 


vom 4. bis 6. Oktober 1954 ein Internationales Symposium tiber experimentelle 
Meteorologie abhalten. Im Rahmen dieses Kolloquiums sollen in erster Linie die 
Fragen der natiirlichen Niederschlagsbildung und ihre kiinstliche, Beeinflussung 
besprochen und im besondern das Problem der Hagelbildung erértert werden, 
Fiir dieses Symposium sind eine beschrankte Zahl von Hauptvortragen ins Pro- 
gramm aufgenommen worden, ferner sind Kurzvortrage vorgesehen; doch soll 
den Diskussionen méglichst reichlich Zeit eingeraumt werden. Als Vortragende 
und Diskussionsteilnehmer haben zugesagt: 

EuGENE Botiay; North American Weather Consultants, Altadena, Calif. ; 
Prof. Dr. Mario Bossorasco, Istituto Geofisico Italiano, Milano; Prof. Dr. 
Horace R. Byers, Department of Meteorology, University of Chicago; HENRI 
DESSENS, Direktor des Observatoire Puy-de-Déme, Clermont-Ferrand; L. Facy, 
Ingénieur en Chef, Météorologie Nationale, Paris; E. M. FoURNIER D’ALBE, Geo- 
physical Observatory, Quetta, Pakistan; Dr. CH. JUNGE, Geophysics Research 
Directorate, Cambridge, Mass.; F. H. LupLtam, Department of Meteorology, Im- 
perial College of Science, London; B. J. Mason, Department of Meteorology, 
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Imperial College of Science, London; R. J. MURGATROYD, Meteorological Research 
Flight, Farnborough; Dr. WALTER Rau, Max-Planck-Institut fiir Physik der 
Stratosphare, Weissenau bei Ravensburg; Dr. VINCENT J. SCHAEFER, Munitalp 
Foundation, Inc., Schenectady, N. Y.; G. SOULAGE, Observatoire Puy-de-Déme, 
Clermont-Ferrand; Prof. Dr. FERD. STEINHAUSER, Zentralanstalt fiir Meteoro- 
logie und Geodynamik, Wien; OTTARIO VITTORI, Osservatorio Monte Cimone, 
Sertola, Italien, und ein Mitglied des U.S. Advisory Committee on Weather Control. 

Anfragen sind zu richten an die Adresse des Prasidenten der Kommission: 
Prof. Dr. R. SANGER, ETH., Postfach, Ziirich 23 (Schweiz). 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Rechenmethoden der Quantentheorie. Von S. FLUGGE und H. MARSCHALL. 
1. Teil: Elementare Quantenmechanik, 2. Auflage (Springer, Berlin 1952). 272 S.; 
DM 29.80, Ganzleinen DM 32.80. 


Dieses schon in der ersten Auflage sehr giinstig aufgenommene Werk legt nun 
in zweiter, umgearbeiteter Auflage vor. Die Anderungen, deren Tendenz vor allem 
war, die physikalischen Fragestellungen starker als bisher hervortreten zu lassen, 
diirften das dem Lernenden wie dem Lehrenden gleich wertvolle Buch noch 
wesentlich bereichert haben. Jedenfalls ist der Einschluss von Aufgaben tiber 
die WKB-Methode, iiber das neue Schwingersche Variationsverfahren fiir Streu- 
probleme und anderes sicherlich ein Vorzug, ebenso die vertiefte Behandlung des 
Diracschen St6rungsverfahrens, das in allen Anwendungen der Quantentheorie 
ja eine zentrale Stellung innehat. Die elementaren Mehrk6rperprobleme und das 
Pauli-Prinzip sind neu aufgenommen worden; dafiir wird die Diracsche Theorie 
des Elektrons auf den zweiten Band verspart, der offenbar der relativistischen 
Quantentheorie und der Theorie der Strahlung gewidmet sein soll. Nach dem 
guten Gelingen der Neufassung des ersten Bandes erwartet man mit doppelter 
Spannung das Erscheinen des zweiten. M. R. Schafroth 


Linear Algebra and Projective Geometry. By REINHOLD BAER (Academic 
Press Inc., Publishers, New York, 1952). 318 pp., 19 Figs; $6.50. 


Das in der von P. A. SMITH und S, EILENBERG herausgegebenen Sammlung 
«Pure and Applied Mathematics» erschienene Buch setzt sich zum Ziel, die struk- 
turelle Identitaét der linearen Algebra und der projektiven Geometrie aufzuzeigen 
durch Hervorhebung der algebraischen Formulierung geometrischer Begriffe. 
Im Vordergrund stehen daher zunachst die linearen Mannigfaltigkeiten, die semi- 
linearen und linearen Transformationen sowie die semilinearen Formen, denen 
entsprechend die projektiven Geometrien, die Projektivitaten und Kollineationen 
sowie die Dualitaten zugeordnet sind. Die Satze werden méglichst allgemein for- 
muliert, also ohne Beschrankung auf die Dimension oder den zugrunde gelegten 
K6rper. Damit wird die von verschiedenen Autoren seit der Jahrhundertwende 
angebahnte Entwicklung zu einem gewissen Abschluss gebracht, und es zeigt sich 
die Tragweite der rein algebraischen Methode. Infolge des gesteckten Zieles fehlen 
die Determinanten und die Matrizen, ebensowenig treten Stetigkeitsbetrachtun- 
gen auf. Aus dem Inhalt: IT. The Basic Properties of a Linear Manitold, III. Pro- 
jectivities, IV. Dualities, V. The Ring of a Linear Manifold, VI. The Groups of 


~ 


we 


- 
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a Linear Manifold, VII. Internal Characterization of the System of Subspaces. 
Zahlreiche Detailfragen, Spezialfalle und Erganzungen sind in Anhangen zu den 
einzelnen Kapiteln untergebracht. Im ganzen handelt es sich um ein sehr ver- 
dienstvolles Werk, das auch manche neue Gesichtspunkte aufweist. 

E. Roth-Desmeules 


Compendium of Meteorology, Prepared under the Direction of the Commit- 
tee on the Compendium of Meteorology H.R. Byers, H.E, LanpsBerc, H. WEx- 
LER, B. Haurwi7z, A. F. SpiruHaus, H.C. WILLET, H. G. Houcuton, Chairman. 
Edited by THomas F.MALone (American Meteorological Society, Boston, Massa- 
chusetts 1951). 1334 pp., Fr. 62.40. 


Dieses grosse Werk der American Meteorological Society in Boston zeigt wohl 


wie kein anderes, welche Wandlungen die wissenschaftliche Meteorologie im 


letzten Jahrzehnt durchgemacht hat. Der Herausgeber, Tuomas F. MALong, meint, 
mit dem Eintritt in die zweite Halfte unseres Jahrhunderts standen wir auf der 
Schwelle einer aufregenden Aera in der Geschichte der Meteorologie, in welcher 
bedeutende Fortschritte zu einem besseren Verstandnis der physikalischen Ge- 
setze, welche in der Atmosphare herrschen, méglich sind. Schon die Wahl der 
Themen und das verschieden grosse Gewicht, das den einzelnen Abhandlungen 
gegeben wird, zeigen diese Wandlung. Wir miissen es uns hier versagen, zu schil- 
dern, in welcher Art die 108 Abhandlungen, an denen 102 Autoren aus der ganzen 
Welt gearbeitet haben, in 25 Abschnitte eingeteilt sind; wir kénnen hier nur einige 
grosse Gruppen hervorheben. Die erste grosse Gruppe befasst sich vornehmlich mit 
der physikalischen Meteorologie, wobei die Physik der Wolken besonders eingehend 
behandelt wird. Im folgenden Abschnitt, der Meteorologie der héheren Atmo- 
sphare, der allein aus zwolf grossen Beitragen besteht, werden — nach einem kurzen 
Abstecher in die kosmische Meteorologie — die Dynamik der Atmosphare, die Zir- 
kulationssysteme und sehr ausfiihrlich die Probleme und Methoden der Wetter- 
vorhersage dargestellt. Dreizehn Beitrage sind der tropischen, der polaren Meteo- 
rologie und der Klimatologie gewidmet. Diesen schliessen sich in mehr oder 
minder logischer Folge Abhandlungen an iiber Hydrometeorologie, Meteorologie 
der Meere, biologische und chemische Meteorologie und Staubmessungen in der 
Atmosphare. In einem besonderen Abschnitt werden die Wolken, der Nebel und 
die Vereisung von Flugzeugen behandelt, wahrend sich die letzten Abschnitte mit 
meteorologischen Messinstrumenten, Untersuchungen im Laboratorium (Modell- 
versuchen), mit Radiometeorologie und Mikroseismik befassen. 

Es handelt sich bei diesem Kompendium keineswegs um eine erschépfende 
Behandlung der meteorologischen Materie; es ist also kein Nachschlagewerk im 
eigentlichen Sinne. Der Hauptzweck besteht vielmehr darin, eine Bestandesauf- 
nahme der im letzten Jahrzehnt gesammelten Erkenntnisse vorzunehmen und 
die vielen noch der Lésung harrenden Probleme zu skizzieren. Die gleiche Materie 
ist oft von zwei Autoren bearbeitet worden. Dies geschah absichtlich, obgleich 
damit eine gewisse Doppelspurigkeit unvermeidbar wurde. Aber es ist ebenso 
interessant wie anregend, zu sehen, wie verschieden das gleiche Problem und 
seine Lésungsméglichkeit von verschiedenen Forschern betrachtet wird. 

Es kann keinem Zweifel unterliegen, dass die meteorologische Forschung durch 
die Bediirfnisse des Zweiten Weltkrieges machtig geférdert worden ist. Die Mittel, 
die namentlich in den Vereinigten Staaten von Amerika hierfiir eingesetzt wurden 
und werden, sind nicht nur ungeheuer gross, es wird auch kéine Miihe gescheut, 
das modernste Instrumentarium aus der gesamten Wissenschaft und Technik fir 
die Meteorologie anzuwenden oder neues zu entwickeln, um die sehr komplexen 


7. 
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Vorginge in der Atmosphare zu erfassen. Die in der Meteorologie so lange herr- 
schend gewesene statistische Methode wird offensichtlich immer mehr von der 
physikalischen Betrachtungsweise abgelést. Nicht nur ist die Atmosphare selber 
zum grossen Experimentierfeld geworden, auch die Modellversuche im Labora- 


arhcer anes 


torium spielen eine immer gréssere Rolle. Die experimentelle Meteorologie hat q 


sich ja schon heute zu einem héchst selbstandigen Wissenszweig entwickelt. 
Dass bei dieser stiirmischen Entwicklung der Forschung ein Werk wie das vor- 
liegende Kompendium in mancher Hinsicht einen recht inhomogenen Charakter 
hat, ist ohne weiteres verstandlich. Es ist aber ein grosses Verdienst, dass die fast 
uniibersehbar gewordenen Einzelpublikationen zu abgerundeten monographischen 
Darstellungen verarbeitet worden sind. Auch dass bei allen Beitragen Bedeutung 


und Anwendbarkeit der Erkenntnisse fiir den praktischen Wetterdienst heraus- 


gearbeitet werden, sei besonders hervorgehoben. Keim Meteorologe wird das Werk 
aus der Hand legen, ohne von ihm angeregt zu sein. J.C. Thams 


Imperfections in Nearly Perfect Crystals. By W. SHOCKLEY, J.H. HOLLo- 
MON, R. Maurer, and F. Seitz (John Wiley & Sons, Inc., New York; Chapman 
& Hall, Ltd., London, 1952). 490 pp.; $7.50. 


Die Eigenschaften von Kristallen hangen weitgehend von der Anordnung der 
sie aufbauenden Teilchen (Atome, Ionen usw.) und der Art der gegenseitigen 
Bindung derselben ab. Aber auch der Grad der Idealitat des Gitterbaues ist be- 
kanntlich oft von entscheidender Bedeutung, insbesondere bei den sogenannten 
stérungsempfindlichen (= strukturempfindlichen) Eigenschaften, wie zum Bei- 
spiel die Festigkeit, die elektrische Leitfahigkeit und anderes mehr. Das vorlie- 
gende Buch (Symposium einer vom 12. bis 14. Oktober 1950 abgehaltenen Kon- 
ferenz des «Committee on Solids» des National Research Council U.S.A.) gibt 
einen guten Einblick in die Natur der St6rungen von Kristallgittern, insbesondere 
von Metallkristallen. Ein einleitendes Kapitel von F. SE1tz macht den Leser mit 
den verschiedenen Typen der Gitterunvollkommenheiten bekannt (Phononen, 
Elektronen und Lécher, Excitonen, Leerstellen und Zwischengitterplatz-Atome, 
Fremdstoffatome, Versetzungen) und beleuchtet deren Wirkung auf die Kristall- 
eigenschaften. Anschliessend werden die Versetzungen geometrisch charakteri- 
siert (READ und SHOCKLEY) und ihre Bedeutung fiir die Verformung von Metallen 
erdrtert (BARRETT, READ, Mott, AVERBACH und andere). In einem weiteren 
Kapitel wird die Diffusion unter anderen vor allem durch BRECKENRIDGE und 
ZENER behandelt. Ferner geht eine Reihe von Arbeiten auf die Eigenschaften und 
Auswirkungen der Korngrenzen ein (insbesondere READ, SHOCKLEY, CHALMERS) 
wie auch auf die Bewegungen in denselben (SMoLUCHOWSKI). Die Natur der inneren 
Korngrenzen (Mosaikgrenzen, «subboundaries» usw.) beleuchtet schliesslich 
Guinier. Das ganze Buch ist gut illustriert und vermag dem Leser einen guten 
Uberblick tiber den Stand der heutigen Kenntnisse vom Realbau der Kristalle 
zu geben. Es kann daher jedem, der sich hieriiber informieren méchte, bestens 
empfohlen werden. W. Epprecht 


C. B. Biezeno, Anniversary Volume on Applied Mechanics (H. Stam, 
Antwerpen, Djakarta 1953). 328 S.; Hfl. 20.-. 

Am 2. Marz 1953 feierte C. B. BrezENo seinen 65, Geburtstag. Ein kleiner 
Kreis von Freunden, Mitarbeitern und Schiilern hat ihm die vorliegende Fest- 
schrift gewidmet, die einen erfolgreichen Lehrer ehrt, aber auch einen der bedeu- 
tendsten Forscher auf dem Gebiete der Elastizitatstheorie, den Mitbegriinder der 


' 
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_ Internationalen Kongresse fiir theoretische und angewandte Mechanik und nicht 


zuletzt eine feine und kultivierte Persénlichkeit. 
Die Festschrift wird durch eine kurze Lebensbeschreibung des Jubilars sowie 
durch ein vollstandiges Verzeichnis seiner Publikationen eingeleitet. Im weiteren 


_ enthalt sie die folgenden Beitrage: 


\ sd 
~ 


W. Boomstra, Triangles équilatéres inscrits dans une conique donnée; H. BRE- 
MEKAMP, Sur la théorie de Sturm-Liouville; TH. Von KARMAN and G. MILLAN, The 
thermal theory of constant pressure deflagration; J. M. BURGERS, Some vemarks on 
detonation and deflagration problems in gases; C. KONING, Some interference pro- 
blems; R. V. SOUTHWELL and G. VatsEy, A problem suggested by Saint-Venant’s 
‘Mémoire sur la torsion des prismes’; R. GRAMMEL, Nichtlineare Schwingungen mit 
unendlich vielen Fretheitsgvaden; R. J. LEGGER, The d’Alembert principle; F. J. 


-PLANTEMA and W. J. VAN ALPHEN, Compressive buckling of sandwich plates having 


various edge conditions; A. D. DE Pater, La stabilité d’un dicéne se déplacant sur 
une voie en alignement droit; A. VAN DER Nevt, The local instability of compression 
members, built up from flat plates; J. A. HaRINGx, Stresses in corrugated dia- 
phragms; A. VAN WIJNGAARDEN, Ut tensio sic vis; D. DRESDEN, Shrink-fit used 
to transmit a torque, W. T. KorterR, On partially plastic thick-walled tubes; C. G. J. 


| VREEDENBURGH and O. STOKMAN, Some new elements in the calculation of flat 


slab floors; J. P. Mazure, Statistical problems in the code of practice for steel win- 
dows; J. J. Kocu, The Laboratory for Applied Mechanics at the Technological 
University of Delft; R. G. Boiren, The design of diaphragms for pressure measur- 
ing devices, based on the use of wire-electrical stvaingauges. H, Ziegler 


Introduction to Measure and Integration. By M. E. Munroe (Addison- 
Wesley Publishing Co., Cambridge, Mass., 1953). 310 pp.; $7.50. 

Cet excellent ouvrage s’adresse a des étudiants avancés, ayant déja suivi un 
cours sur la théorie des fonctions de variables réelles. Ils y trouveront une expo- 
sition didactique de la théorie moderne de la mesure et de l’intégrale. De nom- 
breux exemples et contre-exemples, des problémes allant du facile au difficile 
sont insérés dans chaque chapitre. Plusieurs apergus sur les relations entre la 
notion de probabilité et celle de la mesure, sur les espaces de HILBERT et de 
BANACH, sur la moyenne ergodique et sur la mesure de HAAR dans les groupes, 
donnent une premiére orientation sur des domaines qui sont encore l’objet de 
recherches. 

Table des matiéres: I. Point Set Theory (1-71). II. Measure. General 
Theory (72-114). III. Measure. Specific Examples (115-144). IV. Measurable 
Functions (145-164). V. Integration (165-218). VI. Convergence Theorems 
(219-266). VII. Differentiation (267-302). Bibliography, Index of Postulates, 
Index of Symbols, Index (303-310). M. Plancherel 


Inhalt und Mass. Von Kart MAYRHOFER (Springer-Verlag, Wien 1952). 
269 S., 17 Abb.; broschiert sFr. 37.~, Ganzleinenband sFr. 40.-. 


Sowohl die klassischen als auch die modernen Integrationstheorien fussen 
auf der Theorie der Inhalte und Masse (totaladditiver Inhalt auf einem Sigma- 
mengenkérper). Das Buch von K. MAyRHOFER gibt einen liickenlosen Aufbau 
dieser Theorie unter Wahrung der gréssten Allgemeinheit. Zunachst wird der 
abstrakte Inhalts- und Massbegriff axiomatisch eingefiihrt, um dann, ausgehend 
vom anschaulichen, elementaren Inhaltsbegriff, den Jordanschen Inhalt sowie 
das Borelsche und das Lebesguesche Mass zu definieren. Nach einem Kapitel 
iiber das Verhalten von Inhalt und Mass bei linearen Transformationen wird die 


>. 


=, 


Verallgemeinerung auf Somenfunktionen. 


Das Buch gibt einen trefflichen Uberblick iiber dieses grundlegende Gebiet : 


a 
Jz 


der Mathematik und fiihrt bis zu den Ergebnissen der aktuellen Forschung. Es 


zeichnet sich durch einen strengen, klaren Aufbau aus. E. Roth-Desmeules — 


Introduction to the Differential Equations of Physics. By L. Hopr; 
translated by W. NeF (Dover Publications, New York, 1948). 154 pp.; $1.25. 


Dans cet ouvrage de 154 pages petit format, dont la lecture n’exige que la 
connaissance des éléments du calcul différentiel et intégral, l’auteur expose les 
notions fondamentales de l’analyse vectorielle. Il en donne des applications 4 la 
mécanique et a la physique et il intégre dans quelques cas simples les équations 
différentielles et aux dérivées partielles auxquelles ces applications conduisent. 

Comme il convient dans un livre qui ne veut étre qu’une «introduction», 
Vauteur fait davantage appel a l’intuition géométrique et physique qu’a la rigueur 
du formalisme mathématique (le paragraphe consacré a l’intégrale de FOURIER va 
cependant trop loin dans ce sens en introduisant sans explication l’intégrale 


foe) 


[cose x du 
0 


qui n’a pas de sens). Cette réserve faite, il faut féliciter l’auteur d’avoir réussi 
a donner sous une forme limpide un exposé dont la lecture peut étre recommandée 
au futur ingénieur ou physicien. 

Contenu: I. The Differential Equation as an Expression of a Law of Nature 
(1-4). II. The Ordinary Differential Equations of the Mechanics of Particles 
(6-38). III. The Simplest Partial Differential Expressions (39-63). IV. The Sim- 
plest Partial Differential Equations (64-83). V. Solution by Eigenfunctions 
(84-123). VI. Solution by Change of Variables (124-137). VII. Solution by the 
Use of Singularities (138-152). Index (153-154). M. Plancherel 


Signal, Noise and Resolution in Nuclear Counter Amplifiers. By 
A. B. GILLESPIE (Pergamon Press Ltd., London 1953), 155 pp., 59 figs.; 21s. 


In jedem kernphysikalischen Laboratorium ist man heute auf geeignete Ver- 
starker angewiesen, insbesondere wenn Ionisationskammer, Proportionalzahlrohr 
oder Szintillationszahler zum Nachweis und zur Untersuchung von Kernprozessen 
Verwendung finden. An solche Verstarker sind spezielle Forderungen zu stellen. 
So soll das Verhaltnis von Signal zu Stérpegel, bedingt durch Widerstandsrau- 
schen, Schroteffekt, Funkeleffekt usw., m6dglichst gut sein und gleichzeitig ein 
ausreichendes zeitliches Aufl6sungsvermégen erreicht werden. 

In dem vorliegenden Werke werden zunichst die elektrischen Signale der in 
Frage kommenden Detektoren zum Nachweis von Partikeln diskutiert und an- 
schliessend die auftretenden Stérquellen behandelt. Dies bildet die Grundlage zur 
zweckmassigen Wahl der ersten Réhre und ihrer Arbeitsbedingungen sowie zur 
Wahl des Netzwerkes, welches die Ubertragungsfunktion des Verstarkers be- 
stimmt. 

Das Biichlein bietet eine Ubersicht iiber diese Probleme und ihre Lésung und 
enthalt eine Reihe von experimentellen und theoretischen Ergebnissen, deren 
Kenntnis allen Physikern und Ingenieuren, die sich mit derartigen Fragen zu 
beschaftigen haben, von Nutzen sein werden. E. Baldinger 
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Masstheorie von CARATHEODORY axiomatisch dargelegt und der Zusammenhang 
mit dem abstrakten Mass untersucht. Schliesslich folgt im letzten Kapitel die — 
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ELEKTRONENMIKROSKOPE 
MIKROTOME 
ELEKTRONENDIFFRAKTOGRAPHEN 
HOCHSPANNUNGSOSZILLOGRAPHEN 


TRUB, TAUBER: ZURICH 


Mathematische Neuerscheinungen aus dem Verlag Birkhauser, Basel 


Einfall und Uberlegung in der Mathematik 


von Dr. B. L. van der Waerden, Professor an der Universitat Ziirich. Sonder- 
druck aus «Elemente der Mathematik», 26 Seiten mit 15 Abbildungen, 
Fr. 4.15 (1954). 


Analytische Geometrie 


von Dr. Wilhelm Blaschke, Professor an der Universitat Hamburg. Zweite 
verbesserte Auflage, 190 Seiten mit 67 Figuren, Ganzleinen Fr. 19.60, 
broschiert Fr. 16.65 (1954). — Aus dem Inhalt: Grundbegriffe, Vektoren, 
Matrizen — Kugeln — Stabe — Tragheitsmomente — Quadriken — Konfokale 
Quadriken — Formelsammlung — Namen- und Sachverzeichnis. 


Projektive Geometrie 


von Dr. Wilhelm Blaschke, Professor an der Universitat Hamburg. Dritte 
verbesserte Auflage, 197 Seiten mit 71 Figuren, Ganzleinen Fr. 19.60, 
broschiert Fr. 16.65 (1954). — Aus dem Inhalt: Einleitung: Gegenstand, 
Geschichtlicher Abriss — Schrifttum — Homogene Zeiger, Kollineation — 
Korrelation — Doppelverhaltnis — Staudts Hauptsatz — Kegelschnitte — 
Liniengeometrie — Quadriken — Nichteuklidische Geometrie — Vierflach- 
paare von Mobius — Waben und Gruppen — Namen- und Sachverzeichnis. 


Ubungen zur projektiven Geometrie 


von Dr. Horst Herrmann, Braunschweig. 168 Seiten mit 90 Figuren, 4 zwei- 

farbige Raumbilder und 2 rotgriine Betrachtungsbrillen, Ganzleinen Fr. 17.—, 

broschiert Fr. 14.— (1952). — Aus dem Inhalt: Projektive Ebene — Die Matrix 

als Figur — Die Figur als Matrix — Kollineation — Korrelation — Inversion — 

Projektiver Raum — Vereinigte Lage — Geraden — Kollineation und Korre- 

lation — Liniengeometrie — Quadriken — Kurven dritter Ordnung — Flache 
dritter Ordnung — Konfigurationen. 


«Das eigenartige und reichhaltige Buch stellt eine wertvolle Erganzung zu 
den Werken ,Projektive Geometrie‘ und ,Analytische Geometrie‘ von 
W. Blaschke dar.» (Elemente der Mathematik) 


Sonderprospekte stehen auf Verlangen gern zur Verfiigung 
Zu beziehen durch jede Buchhandlung 


ARCHIVES OF MATHEMATICS 


ARCHIV DER MATHEMATIK 


ARCHIVES MATHEMATIQUES 


Herausgegeben in Verbindung mit dem Mathematischen Forschungsinstitut in 
Oberwolfach von H. BirHarz, H. KNESER und W. Stss : 


VOL. V — 1954 
ERSCHEINT ALS FESTSCHRIFT FUR 


ALEXANDER M. OSTROWSKI 


mit Arbeiten, die Freunde, Kollegen und Schiiler dem verehrten Lehrer und 
Forscher zum 60. Geburtstag am 25. September 1953 widmeten 


Vol. V (1954) erscheint in zwei Teilen 
Teil I (248 Seiten) im Mai 1954, Teil II (etwa 240 Seiten) im August 1954 


Preis von Vols Vi. *. ice ste 4k oe lien Nee 
Subskriptionspreis bis 31. Juli 1954 Fr. 60.— (DM 60.-) 


Auf vorstehend genannte Preise erhalten Mitglieder mathematischer Vereinigungen 
einen Spezialrabatt von 20% 


Inhalt: Beitrage von Baer, R.; Bellman, R., und Hoffman, A.; Bieberbach, L.; 
Bilharz, H., und Schottlander, St.; Blanc, Ch.; Brelot, M.; Carlitz, L.; Collatz, L.; 
Conforto, F.; Deny, J.; Dieudonné, J.; Eichler, M.; Erdés, P., Herzog, F., und 
Piranian, G.; Fan, Ky; Fantappié, L.; Gdértler, H.; Goldberg, K., Newman, M., 
Straus, E. G:, und Swift, J. D.; Green, J. W.; Grébner, W.; Grin, O.; Haack, W., 
und Hellwig, G.; Hadwiger, H.; Hartman, Ph., und Wintner, A.; Hasse, H.; 
Heinhold, J.; Henrici, P.; Hélder, E.; Hofmann, J. E.; Hoheisel, G.; Jager, A.; 
Kaluza, Th. Jr.; Kneser, H.; Krull, W.; Lammel, E.; Levi, F. W.; Lichnerowicz, A.; 
Lorentz, G. G.; Maier, W.; Meixner, J., und Schaefke, F. W.; Nagell, T.; Neumer, W.; 
Norlund, N. E.; Obrechkoff, N.; Radon, J.; Richter, H.; Rogosinski, W. W.; 
Rothstein, W.; Rutishauser, H.; Schmeidler, W.; Schmidt, Herm.; Seibert, P.; 
Specker, E.; Sperner, E.; Stiefel, E.; Steuerwald, R.; Strubecker, K.; Siiss, W.; 
Taussky-Todd, O.; Tautz, G. L.; Todd, J.; Tricomi, F:;\Turan, Ps Varga, O.; 
van der Waerden, B. L.; Wall, H. S.; Walsh, J. L., und Evans, J.; Warschwski, S. E.; 
Weil, A.; Wever, F.; Wielandt, H.; Wittich, H.; Wu, T. T.; Zassenhaus, H. 
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BIEBERBACH, L.: Theorie der geometrischen Konstruktionen . 
BLASCHKE, W.: Projektive Geometrie .. . 
BLAscukE, W.: Analytische Geometric. ..... 
> Burckuarot, J. J.: Die Bewegungsgruppen. der Evibtaliogranhie 
’ CaRATHLODORY, C.: Funktionentheorie, Bd. 1 geb. Fr. 39.50, Bd. 11 
Caratutopory, C.: Mass und Integral und ihre Algebraisierung . 
Doetscu,G.: Handbuch der Laplace-Transformation. Band I: 
Die theoretischen Grundlagen der Laplace-Transformation . . 
Fister, P.; Kurven und Flachen in allgemeinen Raumen . . .- 
FurTer, R.: Analytische Geometrie der Ebene und des Raumes . 
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Werke der reinen und angewandten Mathematik und Physik 
‘aus dem Verlag Birkhauser - Basel und Sere 


geb. Fr. 18.70 
geb. Fr. 19.60 
geb. Fr. 19.60 


geb. Fr. 33.30. 


geb. Fr. 27.05 
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geb. Fr. 83.20 
geb. Fr. 15.40 
geb. Fr. 24.95 


’ Herrmann, H.: Ubungen zur projektiven Geometrie . .. . . . geb. Fr. 17.— 
Hurwit1z,A.: Mathematische Werke, 2 Bande .......- geb. je Fr. 54.10 
Linpver,A.: Statistische Methoden, 2. Auflage 1951. ...... geb. Fr. 31.20 
Locuer-Ernst, L.: Differential- und Integralrechnung .... . geb. Fr. 49.90 


Locuer-Ernst,L.: Einfiihrungin die freie Geometrie ebener Kurven 
MICHAEL, W.: Ortskurvengeometrie in der komplexen Zahlenebene 
Osrrowsk1,A.: Vorlesungen iiber Diff.- u. Integralrechng., Band I 
(Band III erscheint demnachst) Band II 

ScHLAFtt, L.: Gesammelte mathematische Abhandlungen, Band I 
(Band III in Vorbereitung) Band ITI 

StieFEL, E.: Lehrbuch der darstellenden Geometrie. ...... 

' Trost, E.: Primzahlen ..... 
-VoELKER,D., Doerscu, G.: Die. aweiderpucinadle Laplace-Transf. 
ZWINGGI,E.: "Mersicherungsmathematik Se RTH Sy SO ge 


Physik und Astronomie 
Gomsds, P.: Theorie und Lésungsmethoden des Mehrteilchenpro- 
- blems der Wellenmechanik...\.00 7 = SD jee.gecs. pw een a ke 
NIETHAMMER,TH.: Astronomisch-geogr. Ortsbestimmung . .. . 


geb. Fr. 12.50 
geb. Fr. 13.50 
geb. Fr. 36.40 
geb. Fr. 69.70 
geb. Fr, 56.15 
geb. Fr. 56.15 
geb. Fr. 29.65 
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geb. Fr. 47.85 


geb. Fr. 31.20 ~ 


geb. Fr. 33.30 
geb. Fr. 36.40 


Von DER PAHLEN, E.: Einfiihrung in die Dynamik von Sternsystemen geb. Fr. 39.50 


Wa.pmelrr, M.:; Einfiihrung indie Astrophysik ........ geb. Fr. 49.40 
Kongressbericht der Internationalen Konferenz iiber Kernphysik 
und theoretische Physik, Basel 1949... ........ brosch. Fr. 16.65 
' Ingenteurwissenschaften 
Bopea,E.: Giorgis rationales MKS-MaBsystem ....... « brosch. Fr. 28.10 
GoLpDsTEIN, J.: Die Messwandler . . .... 6, | ani, Set a asia Sie geb. Fr. 29.10 
JAEGER,Cu.: Technische Hydraulik. .. ......4.5.-. geb. Fr. 52.— 
Kuximann, K.; Theoretische Elektrotechnik. Band II: Crandeaue ; 
der Theorie der Wechselstréme u. des einphasigen Transformators geb. Fr. 27.55 
Band III: Grundziige der Theorie elektrischer Maschinen. . . geb. Fr. 74.90 
MEIssNeER, E., und Z1zGLER, H.: Mechanik. Band I: Statik der star- 
ren, fliissigen und elastischen Kérper. . . .......4. geb. Fr. 87.45 
Band II: Dynamik der starren Kérper. .......... geb. Fr. 33.80 
Band III: Dynamik der Systeme... ........4..+. geb. Fr. 46.80 
SANGER, R.:-Ballistische Stérungstheorie ..........-. geb. Fr. 15.10 
Srvss1,R.: Vorlesungen tiber Baustatik, BandI ........ geb. Fr. 40.05 
RusTErRHo1z, A.;Elektronenoptik, BandI: Grundziige . . .... geb. Fr. 31.20 
RicuTter, R.: Elektrische Maschinen BandI_ geb. Fr. 49.40 
Band II geb. Fr. 46.80, Band III geb. Fr. 26.— 
VipMar, M.: Die Gestalt der elektrischen Freileitung ...... geb. Fr. 19.75 
Wacn_r, K. W.: Elektromagnetische Wellen ......... geb. Fr. 33.30 


Zu bexiehen durch Ihre Buchhandlung 


